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Abstract 


This study shows an analysis of a ship power plant with waste heat recovery systems — Thermo 
Efficiency System (TES) with main engine exhaust gases utilization and exhaust power gas and steam turbines 
application. The systems were compared each other. The energetic efficiency in chosen waste heat recovery 
systems was compared. The advantages and disadvantages of turbines application in waste heat recovery systems 
were listed in a summary. 


Keywords: turbine, engine, waste heat energy, recovery, ship, power plant, 
1. Introduction 


Ship’s engine room efficiency determines the stage of the heat, which is received by fuel 
oil combustion in the main components of her power plants. Modern solutions of ship’s power 
plants ensure the biggest possible efficiencies of converting fuel oil chemical energy for other 
aspects. One of the methods increasing efficiency of the power plant is recovering the heat 
which is lost with exhaust gases, main engine cooling water and charging air [2, 3, 4, 5]. 
Solution of this problem is possible through maximization of waste heat recovery generated by 
marine diesel engines and applying it in other installations to produce heat, electric and 
mechanical or combination of these energies. As a consequence of increasing main propulsion 
unit power of modern ships (increasing its deadweight), is forming a surplus energy which is 
produced in exhaust gas boilers as heating steam. In these circumstances, applying the 
combined waste heat recovery systems with exhaust gas and steam turbo generators is rational. 
There are a lot of waste heat recovery systems possible to apply on the motor ships. Their 
configuration depends on type of the vessel, her capacities, operating speeds and the output 
values implemented by the main propulsion plant [2]. Selection of the system solution on the 
specific ship should be result of widely comprehended, penetrating technical-economical 
analysis, based on solid mounts of thermodynamic and reliability analysis — methods keeping 
the engine room mobility. 


2. Waste heat energy availability 


Main and auxiliary engines are sources of waste heat energy in ship’s propulsion plant. 
Quantity of these streams are diversified, energy included inside is not utilized fully. Waste 
heat energy losses size depend on main propulsion plant which has been used (diesel engine 


or/and gas turbine). Low speed diesel engines have the highest thermal efficiency which is 
formed in range between 45+55% (for comparison the turbine engines have lower thermal 
efficiency even 20%) [3, 4, 7]. Waste heat energy include heat losses of an exhaust gases, 
scavenging air and cooling water streams. 

Increasing ship’s energetic power plant efficiency is lashed with taking off assist of huge, 
available waste heat energy streams coming from main components of the system. On the 
pictures no. 1, 2 there are adequately heat balance diagrams of marine low speed diesel engines 
without and with recovery of energy, 12RT-flex96C type [9] and 12K98ME/MC type [8] on 
Fig. 3 and 4. Both of the engines evolve the 68640 kW output power which is realized with 94 
rev/min. 

Shaft power Overall efficiency 49.3% 


Exhaust gases 25,4% 


Scavenge air cooling water 14.1% 


Jacket water 6,3% 
Lubricating oil 4,3% 


Radiation 0,6% 


Fig.1 Heat balance diagram of marine low speed diesel engine 12RT-flex96C type [9] 


In WARTSILA solution, both waste heat recovery system uses exhaust gas stream 
generated by main propulsion plant and scavenging air and cooling water streams. The biggest, 
possible waste heat recovery stream is carried by exhaust gases leaving the main engine. The 
lesser are: scavenge air, cooling water and lubricating oil streams. Systems offered by MAN 
B&W utilize only the exhaust gas stream. Overall efficiency of these engines without recovery 
is equal in total 49,3%. 


Shaft power Electrical power 5.9% Overall 
efficiency 54,9% 


Condenser 8,6% 


Exhaust gas 12,6% 


Scavenge air 
cooling water 12.9% 


Jacket water 6.2% 

Lubricating oil 4,2% 

Radiation 0,6% 

Fig. 2 Heat balance diagram of marine low speed diesel engine 12RT-flex96C type cooperating with waste heat 
recovery system [9] 
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Lubricating oil cooler 2.9% 


Jacket water cooler 5,2% 


Exhaust gas 25,5% 


Air cooler 16,5% 


Radiation 0,6% 


Fuel input 100% 


Fig. 3 Heat balance diagram of marine low speed diesel engine 12K98ME/MC type [8] 


Shaft power Gain 9,9% 


Lubricating oil cooler 2,9% 


Jacket water cooler 5,2% 


Exhaust gas 22,9% 


Air cooler 14,2% 


Radiation 0,6% 


Fuel input 100% 


Fig. 4 Heat balance diagram of marine low speed diesel engine 12K98ME/MC type 
cooperating with waste heat recovery system [8] 


3. Waste heat recovery systems 


There are a lot of waste heat recovery systems with exhaust and steam turbines [7, 8, 9]. 
They are representing various energetic efficiency. 

In this study, the standard, referencing solution of waste heat recovery system is system 
with exhaust power gas turbine, working according to the scheme shown on Fig. 8. The basic 
waste heat recovery system consists of a main engine with high efficiency turbochargers, 
exhaust power gas turbine via reduction gear and clutch driven AC generator. The energy of 
exhaust gases, leaving the main engine, is partly recovered in a turbochargers to compress, 
charge air or supply the exhaust power gas turbine driving the generator. About 10% of exhaust 
gas leaving the main engine is used to drive the generator [8]. The exhaust gases are being 
carried off through the cumulative manifold to the funnel in which there dampers installed and 
then to the atmosphere. 

This solution is an alternative source of producing an electric power and now it is being 
introduced to the new building ships. 

IC 
Exhaust 
gas lines 


with dumper wem 
35,33) 


Tch (Tich: Ptch) 


Turbochargers 


Cluth Exhaust 
power gas 
turbine 


Fig. 5 Waste heat recovery system with exhaust power gas turbine [8] 


Engine working waste heat recovery system (Fig. 8) is equipped with an air waste gate, 
stifle valve aim at keeping the normal cylinder pressure without the other restrictions at the 
very low intake air temperatures. At present, the propulsion plant cooperating with waste heat 
recovery systems can work at the intake air temperatures maintaining at the level of —5 to 35 
centigrade. After exceeding the maximal pressure in a cylinder, the flow rate of exhaust gases 
feeding the exhaust power gas turbine, turbine and generator rotor speed would increase. It will 
cause the generator load disturbance. 

Adopting the main engine to work at lower intake air temperatures involves the costs of 
fuel oil consumption, however the rate of recovered energy compensates this disadvantage. 

Exhaust power turbine is working in the range 55+100% of main engine load. Exhaust gas 
outflow is controlled by orifice on the outlet of the exhaust gas receiver, keeping the constant 
value of exhaust gas flow rate feeding the power turbine. In this way, the constant turbine load 
rate is being kept. With the main engine output lower than 55%, the exhaust gases feeding the 
turbine are being cut off. The air flow which is given by turbochargers is too low to generate 
properly large exhaust gas flow rate to ensure stabilized turbine load rate. Expand stage and 
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exhaust gas flow rate value are compared with suitable turbochargers parameters. Turbine 
outlet exhaust gases temperature is close to the temperature after the turbocharger. 

Basic waste heat recovery system (Fig. 5) is often compared to the system with the exhaust 
gas boiler and steam turbine driving the generator. These systems are different in respect to the 
configuration and outputs. On the picture 6 there is waste heat recovery system with steam 
turbine feeding the steam from exhaust gas boiler [9]. 


Fig. 6 Waste heat recovery system with an exhaust gas boiler and steam turbine [9] 


On the Fig. 7 there is waste heat recovery system with an exhaust gas boiler exhaust power 
gas turbine and shaft generator. This system consists of main engine, turbochargers, exhaust 
gas turbine mechanically connected with the generator, exhaust gas boiler, shaft generator and 
self-contained auxiliary engines feeding the common board grid Described system differs from 
above one in relation to waste heat energy division between turbines. 

Main exhaust gas energy stream: carrying by exhaust gases, after energy conversion in 
turbochargers is divided in two: exhaust gas stream heating the exhaust gas boiler and exhaust 
gas stream feeding the power turbine. During main engine working one part of exhaust energy 
stream is used to produce heating steam in exhaust gas boiler. Second part of this exhaust gas 
stream is converted by power turbine set mechanically connected with the generator to electric 
energy for the board grid. Controlling of the exhaust gas amount feeding the power turbine is 
realized in relation of main engine load through an orifice. About 10% of exhaust gases energy 
are used in this waste heat recovery system in power turbine. 

Exhaust gases leaving the exhaust power gas turbine are guided to exhaust gas boiler too. 
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Exhaust 
Generator — power 
gas turbine 
Steam 
Exhaust > for heating 
gas boiler purposes 
Auniliary engines 
1 | | 
p H 
- 7 , | 
Turbochargers la Gi EJ es Da 
A 4 a 
Shaft ; 
generator 
Switchboard 


Fig. 7 Waste heat recovery system with an exhaust gas boiler, power turbine and shaft generator [8] 


The efficiency of whole system is increased by a shaft generator application. This 
generator can work in a two variants. It can taking off the energy to the board grid or taking in 
assist the main propulsion plant. Engine running variant can be used for reduce fuel oil 
consumption amount with the reduced main engine output. The second variant is electric 
energy production for ship requirements, while the waste heat recovery system is stopped. 

Beyond supplying the electric power, this system, allows to save the fuel oil (even to 5% 
per year) and essentially reduces the emission of toxic combustion substances. 

Advanced waste heat recovery system, shown on Fig. 8 consists of an exhaust gas boiler, 
steam turbine, exhaust power turbine, synchronous generator driven by these turbines and shaft 
generator working for the board grid. Exhaust power gas turbine impeller connected with steam 
turbine impeller via power transferring system (reduction gear and clutch) drives the AC 
generator. Steam turbine is fed by superheated steam from the exhaust gas boiler. Shaft 
generator application increases whole efficiency of the system. This solution has been used for 
a first time on m/v Gudrun Maersk containership [9]. 
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Reduction Reduction gear 


ear with clutch 
SĄ Steam | Emergency 
Exhaust Generator \ turbine | generator 
È \ | f CH 
gas boiler | Exhaust l 
| power gas 
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-OH | 
generator in | 


\ Switchboard 


Fig. 8 Advanced waste heat recovery system with an exhaust power turbine, steam turbine and shaft generator [9] 


Combined system application with en exhaust and steam turbine considerably increases its 
efficiency and reliability. Due to increase of an electric energy production at main engine loads 
above 50%, system efficiency grows even to 10% (see Fig. 9). 

One possibility to improve the efficiency of the recovery of more heat without increasing 
the heat exchange surfaces can be applied fluidized bed exhaust boiler [1]. 


4. Chosen aspects of turbine application advisability in waste heat recovery systems 


Economical efficiency of examined waste heat recovery systems with exhaust and steam 
turbines can be measure of its quality solution assessment [6]. Its final measurement will 
always be derivative brief fore designs accepted for the specific ship, conditions of her 
operating and stabilized technical — economical criterion. In this study, it was restricted to some 
aspects of turbine applications in waste heat recovery systems [5, 6]. 

Profitability of use and operate turbine in waste heat recovery system is conditioned by 
vessel state operating (sailing speed and the power involved by the main propulsion plant)and 
the time of being in its, which determines the quantity of energy possible to utilize. Modern 
waste heat recovery systems, with gas and steam turbines or their combination, are different 
due to their configuration, machines property and working parameters. 

With increasing the exhaust power gas turbine output, disposal waste heat recovery 
streams increases in waste heat recovery system. The power turbine output increases and as a 
consequence the amount of the saved fuel combusted by propulsion plant increases too [9]. On 
Fig. 9 there are example ranges of possible output values reached by exhaust power gas turbine 
in relation to main engine output type 12K98ME/MC manufactured by MAN B&W [8]. 
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Exhaust power gas turbine output [kW] 
N 
DA 
8 


2000 
1500 
1000 
500 : 
41,130 kW i 
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Main engine power output [kW] 


Fig. 9 Exhaust power gas turbine output waste heat recovery system in relation to main engine power [8] 


With the main engine power 41130 kW exhaust power gas turbine working in TCS — PTG 
system evolves repayable power 2260 kW. 

As a result of extracting specified outputs, there are fuel oil savings shown for the same engine 
on Fig. 10 in relation to exhaust power gas turbine output and main engine operating hours [8]. 


1650kW, 6000h Pon 
Oszczędność 2000t paliwa ;„+** 
zeta te data zac zat zadka L 


1650kW 


0 500 1000 1500 2000 2500 3000 3500 4000 
Exhaust power gas turbine output [kW] 


Fig. 10 Prediction fuel oil savings obtained as a result of waste heat recovery application in relation to exhaust 
turbine output [8] 
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In a case of waste heat recovery system with exhaust power gas turbine output 1650 kW, after 
6000 working hours, the 2000 tons of fuel oil can be saved in a year. Such rational advantage 
coming from waste heat recovery system application can be gained in case of huge main 
propulsion plants, here exceeding 40000 kW [8]. 


5. Thermodynamic analysis of examined waste heat recovery system 


Exhaust gases ability was analyzed to conversion part of the waste to effective electric 
energy with intermediate aspect reciprocating energy generated by power turbine. On the Fig. 5 
there is basic waste heat recovery system with alluvial points enable the thermodynamic 
identification. Inlet exhaust gas parameters marked as 1, outlet parameters as 2. 

Large flow rate and exhaust gases temperatures relatively keeping at the same level ensure 
working stability of power turbine and constant internal efficiency close to project value. On 
table 1 there are technical data of radial and axial power turbines applying in waste heat 
recovery systems manufactured by MAN B&W [8]. Expand stage turbine 7; is 3,3, exhaust gas 
temperature before turbine is 450°C. 


Tab. 1. Basic technical data of power turbines 


| Radial power turbines | 
Output [kW] Flow rate [kg/s] Speed [1/min] 
(TCS-) PTG16 600 3.0 41,000 
(TCS-) PTG18 900 4.3 34,000 
(TCS-) PTG20 1,300 6.3 28,500 
(TCS-) PTG22 2,300 13.0 21,500 
(TCS-) PTG55 3,400 19.0 17,000 
(TCS-) PTG66 4,800 27.0 14,500 


They can work in system TCS/PTG (Turbo Compound System/Power Turbine Generator) — the 
system taking off main engine assistance/exhaust power generator set. 

For the turbines described in table 1 the disposal enthalpy drops were calculated according 
to the dependence [3, 4, 6]: 


P 
Mio mam (1) 
Moch ` Nir ` Vint 
and heat streams [4, 6]: 


AQ, = UM j Ah, (2) 


The results of calculations are show in table 2. The biggest heat streams convert axial power 
turbines whereas disposal enthalpy drop values render that it can be even single stage turbines. 
Bigger disposal enthalpy drops are converted by radial power turbines, in comparison with 
axial they give less outputs. Higher output values gained by axial turbines are result of bigger 
flow rates exhaust gases conversion. Bigger turbines internal efficiencies are fostering in 
gaining higher axial turbine outputs, which influences on waste heat recovery system 
efficiency. 
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Tab. 2. Calculated dispose enthalpy and heat streams in turbine 


Radial power turbines 
Turbine type Ahsr [kJ/kg] AQ, [kW] 
TCS -PTG 16 266,6 799,8 
TCS -PTG 18 279,1 1200,13 
TCS-PTG 20 275,1 1733,13 
TCS -PTG 22 235,9 3066,7 
Axial Power turbines 
TCS -PTG 55 210,5 3999,5 
TCS — PTG 66 209,1 5645,7 


5. Summary and conclusion 


Heat balance diagrams analysis of marine low speed diesel engines confirms possibility 
and advisability use of waste heat recovery energy in marine energetic power plants. Marine 
engine power output, temperature of medium feeding the waste heat recovery system (exhaust 
gases, cooling water), pressure and flow rate factor determine the quantity of available heat in 
waste heat recovery system. Exhaust gas and steam turbine application takes effect in 
additional electric energy production, fuel oil savings and limitation of emission toxic 
substances to the atmosphere [8]. Exhaust gas and steam turbine cooperating increases systems 
efficiency especially in direction to enlarge of electric energy production. Use character of 
waste heat recovery system with power turbine is conditioned by continuity of steady turbine 
work at main engines load ranges above 50% rated power ensuring adequately large waste heat 
streams carrying by exhaust gases. Quantity of saved energy can be estimated after penetrating 
analysis of each other of waste heat recovery systems. 
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Abstract 
This paper presents proposal of a method for determining operational load parameters of dredge pumps, one of 
the crucial mechanical energy consumers on dredgers. The method based on results of the author's operational 
investigations, deals with two main service states of dredge pumps on dredgers, namely: the state of loading the 
solid into soil hold (of dredger or hopper barge) and the state of pumping ashore. 
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1. Introduction 


Dredge pumps belong to the most important consumers of mechanical energy on 
dredgers. Their function is to hydraulically transport loosened soil from the sea bed into soil 
hold of the dredger or hopper barge (a service state called loading the spoil) as well as from 
the hold (sometimes directly from the sea bed) through long transfer piping to a dump on 
shore (a service state called transferring the spoil ashore). The states occur always on suction 
dredgers (e.g. trailing suction hopper dredgers, cutter suction dredgers, barge unloading 
dredgers), sometimes also on dredgers with mechanical dredging systems (e.g. bucket ladder 
dredgers) [8]. Power demand of dredge pumps depends on their use and design assumptions 
as well as on size of dredger. It is contained in a broad interval ranging from several hundreds 
kW to even a dozen or so thousands kW [7,8]. 

Irrespective of a type of dredger the dredge pumps can operate in two basic service states 
[3,4,6,7,8]: 

— the loading of the spoil into the hopper (soil hold) on the dredger or assisting hopper 
barge; operational conditions of the pump system are characterized by the following 
features: the static lifting height of the system is as a rule greater than the dynamic one 
(H,,2H,,, ), similar values of flow drag of water-soil mixture on suction and 


pressure side of pump (Ah, x Ah, ); 


— the hydraulic emptying of the soil hold or transferring the spoil directly to a dump on 
shore (the pumping the spoil ashore); in this case operational conditions of the pump 
system are characterized by a much greater dynamic lifting height than the static one 
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(H,,<<H,,,,) and much greater values of flow drag on the pump pressure side than on 
its suction side (Ah, << Ah, ). 


The state called loading the spoil into the hold always occurs on trailing suction hopper 
dredgers (it concerns their own holds) , and may also occur on cutter suction dredgers (in this 
case it concerns hopper barge holds). The pumping-away the spoil occurs on trailing suction 
hopper dredgers and cutter suction dredgers, sometimes also on bucket ladder ones. Great 
differences in the parameters which characterize the pump systems operating in the above 
mentioned service states must result in great differences in the loads applied upon dredge 
pumps during loading and pumping-away the spoil. On the trailing suction hopper dredgers 
the using of the same pumps both for the loading and pumping-away the spoil is common. 
Then their driving systems are fitted with multi-speed gear transmission devices. 

This paper presents a proposal of determining the distribution parameters of operational 
loads of dredge pumps installed on various types of dredgers, depending on their service 
states. For the determining of the distribution parameters of pump loads the use of linear 
form of dependence of mean driving loads of dredge pumps on their rated power outputs, is 
proposed [1,2]. Standard deviations are proposed to be determined by means of the data 
concerning the load distribution variability coefficient of dredge pumps, obtained on the basis 
of operational investigations. 


2. Investigations of relations between mean operational loads of dredge pumps and their 
rated power outputs 


Investigations of relations between the mean operational loads of main power consumers 
Ne, , and their rated power outputs (N;/.)””” are important in view of possible making use 
of their results further in preliminary design stages of ship power plants. 


As far as the main consumers are concerned, it can be considered effective energy flow 
(i.e. power output or effective power) associated with them, each case expressed by the 
product of the so called ,, generalized potential” and the ,, generalized flow” [1]. 


The investigations of relations between the mean operational loads of main power 
consumers and their rated power outputs were already performed for fish factory trawlers as 
well as certain main consumers on dredgers [1,2]. 


The performed investigations show that for all main power consumers on the investigated 
ships the following linear relation takes place [1,2]: 


eff _ eff \nom 
Nic =a + b (Nyc) (1) 
where: a, b - constants 


The statement on validity of the linear relation (1) as well as on possibility of 
determining the constants a, b, is very important as the relation can be used for determining 


the power demand of main engines for ships of the considered types in the stage of offer or 
preliminary design [1,2,3]. 


The dredge pump power output N;/, is determined by the relation: 
Nop = pp Qpe (2) 


where: H, - the dredge pump lifting height determined for water , 


Opp- the dredge pump volumetric rate of delivery determined for water. 
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It is important that producers of dredge pumps usually provide nominal (rated) 
parameters and characteristics of the pumps valid for the conditions of water pumping but not 
water-soil mixture pumping. Change of characteristics of dredge pump handling water-soil 
mixture, as well as change of pipeline characteristics is most influenced by density of the 
mixture and soil graining (a.o. mean grain diameter, grain-size distribution, grain shape) [7,8]. 
High variability of the parameters causes that the providing of the parameters of dredge 
pumps for the conditions of water-soil mixture pumping would be unjustifiable. 

In Tab. 1 and 2 are given the rated parameters of the dredge pumps on the investigated 
dredgers as well as the distribution parameters of operational loads of the pumps during 
operation in two basic stages of their service (loading the spoil to hopper and pumping the 
spoil ashore) [6]. The distribution parameters of operational loads of dredge pumps have been 
obtained as a result of long-lasting operational investigations carried out by this author on a 
dozen or so dredgers of various types. The problem of operational loads of dredge pumps, 
which covers measurement methods, measurement system characteristics, distributions of 
pump loads, has been presented more thoroughly in a few publications of this author [3,4,5,6]. 


Tab.1. Rated parameters of dredge pumps and characteristics of their operational loads during loading the spoil 
to hopper, for 8 investigated dredgers 


Rated parameters of dredge pumps Characteristics of load 
piace: Number - = water) — ue of dredge pumps 
pumps App DP Nop Nop Opp Vpp 
kPa m/s kW kW kW - 

Kostera 1 105 0,65 68,2 82,1 6,9 0,084 
Kronos 1 140 0,65 91 95,1 7,8 0,082 
Łęgowski 2 175 1,8 2x315 774,9 62,3 0,08 
Bukowski 2 175 1,8 2x315 786,4 54,5 0,069 
Nautilus 1 210 2,5 525 706,8 26,2 0,037 
Gogland 2 220 3,2 2x704 1787,1 61,2 0,034 
a e 2 265 3,4 2x901 1987,7 71,9 0,036 
ie 1 395 4,6 1817 2482,6 


Tab.2. Rated parameters of dredge pumps and characteristics of their operational loads during pumping the 
spoil ashore, for 13 investigated dredgers 


Rated parameters of dredge pumps Characteristics of load 
Rise Number (pumping water) distributions of dredge pumps 
pumps Ay. Ee NZ, Nop Opp Vpp 

kPa mĉ/s kW kW kW - 
Kostera 1 370 0,55 203,5 216,5 29,6 0,137 
Kronos 1 490 0,5 245 189,4 20,6 0,109 
Łęgowski 2 385 1,6 616 832,9 101,4 0,122 
Bukowski 2 385 1,6 616 835,7 52,8 0,064 
Gogland 2 430 3,0 1290 1815,9 178,4 0,098 
Geopotes 15 2 560 3,0 1680 2211,1 239,9 0,109 
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renee 2 1240 4,1 5084 3861,9 
Trojan 1 600 1,0 600 536,7 122,9 0,229 
Torun 1 530 0,95 503,5 431,5 52,3 0,121 
Scorpio 1 610 2,05 1250,5 1667,5 278,6 0,167 
Rozkolec 2 1160 1,05 2030 1188,1 334,8 0,282 
Raja 1 440 0,7 308 240,5 32,3 0,134 
Małż II 1 410 0,4 164 158,26 30,78 0,194 


Tab.3. Linear regression equations which determine mean loads of dredge pumps during two basic states of 


their service 


States of dredge Postać zależności Statistical evaluation 
pumps o F F, | m 
loading of the spoil to (AI =1,242 - (A) 5 745 0,988 | 81,4 | 276,6 | 5,99 | 8 
hopper 
pumping the spoil ashore | (NE)? =0,755-(N;)” +244,43 | 0,934 | 147,4 | 75,8 | 4,84 | 13 
3000 7 
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2000 + 
> 1500 
z 
5 A 1000 
Z 
500 
0 t U T U 
0 200 400 600 800 1200 1400 1600 1800 2000 


WALIA 


Fig. la. The relations Np, = f(N Pj for dredge pumps in state of loading the spoil to hopper 
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Fig. 1b. The relations Np, = f(N EA) for dredge pump in state of pumping ashore 


If rated parameters of dredge pump (-s) are known it is possible - by making use of the 
relations given in Tab.3 - to predict its (their) mean load during a considered service state. 


3. Standard deviations of operational load distributions of dredge pumps, working time 
fractions of the pumps in a given service state 


Value of the standard deviation cpp can be determined by using data which deal with the 
variability coefficient v,, =o,»/N;p (Tab.1 and 2). Values of the variability coefficients of 
load distributions of dredge pumps during loading the winning are contained within the 
interval of 0,034-+0,084 at the mean value of v,,equal to 0,06 (Tab.1). In the case of the 
service state of transferring the winning values of the variability coefficients of load 
distributions of dredge pumps are contained within the interval of 0,064 + 0,229 at the mean 
value of Vp, equal to 0,146 (Tab.2). 

The analyzed loads of dredge pumps concern duration time of loading the spoil to the 
hold (of the dredger or hopper barge) or hydraulic unloading the spoil from the hold (of the 
dredger or hopper barge), that is determined by values of the usage time coefficients 25, 2%, 


(Tab.1 and 2 given in [6]).Values of the coefficient 15, are contained in the interval of 


0,95+0,98, at its mean value equal to 0,972. Values of the coefficient 2%, are close to the 45, 
values and are contained in the interval of 0,96--0,99 , at its mean value of 0,978. 

If the relation of the duration time of „dredging operations” has to be determined the 
coefficient A” or Ae is to be additionally taken into account. Values of the coefficients 


depend on a type of dredger and its design solution. 

In addition the case of cutter suction dredger with underwater pump should be 
highlighted. The so applied pumps operate both during loading the winning into hopper barge 
(dredging to hopper barges) and during dredging with simultaneous pumping ashore with the 
use of dredge pump (pumps) installed onboard. For both the situations working conditions of 
the underwater pump can be assumed the same and corresponding with the conditions of 
loading the spoil. Of course, the dredge pump installed onboard operates only during 
pumping the spoil ashore. 
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4. Summary 


All the calculation results which concern load characteristics of dredge pumps, presented 
in this paper, reflect operational reality of the pumps on dredgers. The results can be deemed 
representative for the whole population of dredge pumps used on dredgers, in view of the 
large number of the investigated dredgers and wide range of their size. 

The presented results may be useful in predicting operational loads of dredge pumps, 
depending on their service states typical on various types of dredgers. That will make it 
possible — in association with knowledge of loads of other main consumers and efficiency 
characteristics of power transmission systems of particular consumers - to determine 
characteristics of operational loads of main engines on dredgers. It is especially important in 
preliminary design stages of power systems for dredgers. 
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Abstract 


The paper deals with the problem of energy demand for the main propulsion as a function of deadweight 
and speed for general cargo vessels built in the 60-ties, multi-purpose general cargo vessels built In the 80-ties 
as well as recently built container vessels. Changes in power of the main propulsion and trends observed in the 
matter are defined. In the same way the analysies of electric power and boilers capacity are carried out. In the 
summary conclusions and prognosis concerning energetic plants of general cargo and container vessels are 
expressed. 


Keywords: cargo ship, container ship, main propulsion power, electrical power, auxiliary steam delivery, 
Statistics 


1. Introduction 


General cargo vessels built in the 60-ties were designed for carriage of general cargo i.e. 
industrial goods counted in number and packed in boxes, drums, bales, bags or other similar 
packages. Such vessels were usually provided with cargo handling equipment (deck cranes, 
cargo booms) to make possible cargo operation without aid of harbour cargo facilities. Cargo 
space of the vessel was divided by bulkheads and twin decks to optimise the space utilisation 
and to separate different kinds of cargo as well as to separate cargo of different destinations. 
The speed of these vessels was 16 +17 knots. At the turn of the 60-ties and the 70-ties a new 
kind of general cargo vessels fitted for container carrying named universal or multi purpose 
general cargo vessels were built. They achieved service speed of 18+20 knots. In the mid of 
70-ties traditional general cargo as well as multi purpose general cargo vessels were almost 
entirely replaced by containerised cargo and ships named container vessels. 


Container vessels are the ships equipped with special guides to freight vertically loaded 
and unloaded containers. The first ship adapted in 1956 from tanker for container transport 
was Ideal-X. Nowadays, the biggest container vessels can carry above 10,000 TEU and can 
hardly be situated in panamax class dimensions. The number of those ships is bigger than 
panamax class and they are used on routes excluding Panama Canal passing e.g. China — USA 
West Coast. Today the biggest container vessel is MS EMMA MAERSK with carrying 
capacity of 11,500 TEU shown in fig. 1.To ensure quick transport of containers all around the 
world and to minimise transport expenses large container vessels travel overseas calling at a 
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number of big ports called hubs. Containers are delivered from smaller ports to hubs by small 
container vessels (200-500 TEU) named feeders. A number of small container vessels (below 
3000 TEU) is equipped with cargo handling facilities (deck cranes or bridge cranes), so they 
are able to call at ports not equipped with cargo facilities. Bigger container vessels have to use 
harbour equipment. Regarding vessel dimensions and weight of containers a special goliath 
gantry cranes are used. Container vessels are the quickest freighters. They achieve service 
speed of 24-26 knots. 


Fig.1. M/S EMMA MAERSK at sea 


An initial analysis shows as follow: 

- the main propulsion is executed by low speed diesel engines; in case of modern 
container vessels steaming with high speed very large and powerful diesel engines of 
100,000 HP and bigger shaft power are used, 

- three diesel generators create onboard power station on general cargo vessels; 
sometimes an additional shaft generator is used; in case of container vessels the 
onboard power station is considerably bigger due to the necessity of bow thrusters 
and refrigerated containers supply, 

- a boiler room usually consists of two auxiliary steam boilers one fuel fired and the 
second one heated by main engine exhaust gases; capacity of boilers 2000-3000 kg/h 
and considerably higher on big container vessel. 

The aim of this paper is an analysis of the trends in development of main propulsion 
power, electric power and auxiliary boiler capacity on general cargo and container vessels by 
means of statistics. 


2. Analysis of main propulsion plants development 


To determine main propulsion power of general cargo vessels built in the 60-ties the data 
of 287 such ships was taken to considerations [3]. The result of statistic calculations is the 
formula (1) [3] [5] which shows the dependency of main propulsion power on deadweight and 
ship service speed: 


N, = 0,039-D°* .v*** [kW], (1) 


where: Na[kW] — shaft power of the main engine, 
D [tons] — deadweight, 
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v [knots] — ship service speed. 
Statistic researches of container vessels built in the last years gave formula (2) [2]: 


N,= (0,99179 + 0,00003412 - D) v? [kW], (2) 


On the basis of formulas (1) and (2) the dependency of main propulsion power on deadweight 
and different service speed for general cargo vessels is shown in figure 2 and the same for 
container vessels is shown in figure 3. To enable the comparison of main propulsion plants 
development in last years for both cases, the approximation of both formulas was done by 
means of linear functions. An example of the final analysis concerning main propulsion 
power of both types of vessels for given speed v=18 knots is shown in figure 4. 
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Fig. 2. Dependency of general cargo vessels main propulsion power on deadweight for different speed 
Nn = 0,0287 * D°* * y°"8 [kW] DWT [tons], v [knots] 


26knots 


24 
50000 
22 
40000 Za 
20 
30000 
18 
20000 
16 
——— — ia 


0 
0 20000 40000 60000 80000 100000 
DWT[tons] 


Fig. 3. Dependency of general container vessels main propulsion power on deadweight for different speed 
Nn = (0,9179 + 0,00003412 * D) * vê [kW] DWT [tons], v [knots] 
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Fig. 4. Comparison of main propulsion power of general cargo vessels and container vessels 
at service speed 18 knots 


As the largest built general cargo vessels achieved deadweight of about 20,000 tons to 
make it possible to compare with container vessels an extrapolation of power describing 
function was executed as a prognosis. The result of comparison (fig. 4) shows that for vessels 
larger than about 15,000 tons the main propulsion energy demand of general cargo vessels is 
much bigger than the same of contemporary built container vessels. It is probably possible 
due to the improvement in hull construction. General cargo vessels built in the 60-ties and the 
70-ties did not have bow bulb and modern stern construction thus the resistance of ship hull in 
motion was higher and the same for main propulsion power demand. On the other hand, a 
high main propulsion power of container vessels is the result of high service speed because 
container vessels are the fastest merchant vessels. 


3. Analysis of onboard power station development 
The total electric power of general cargo vessels onboard power station can be 


approximately estimated from formula (3) elaborated by Centrum Techniki Okrętowej in 
Gdansk [4]: 


Z Ny= 23 + 0,1088 N, [kW], (3) 
where: © Na[kW]-— total electric power, 
N, [kW] — main propulsion power (main engine shaft power). 


Next, the total electric power of container vessels power station is given in formula (4) 
described in method [1]: 


Z Na= 1077 + 0,1580 N, [kW], (4) 


The dependencies of total electric power on main propulsion shaft power for general cargo 
vessels and container vessels from formulas (3) and (4) are compared in figure 5. It is possible 
to observe much bigger demand of electric power for container vessels than for general cargo 
vessels. However some general cargo vessels were equipped with cargo handling facilities 
which increase electric power demand, container vessels have bigger electric power demand 
due to installation of bow thrusters and the necessity of power supply for a big number of 
refrigerated containers. 
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Fig. 5. Comparison of total electric power installed on general cargo vessels and container vessel 
as a function of main propulsion power 


4. Analysis of auxiliary steam boilers development 


To determine total auxiliary boilers capacity on general cargo vessels formula (5) given in 
[2] was used: 
D; = 0,075045 N, + 1054,5 [kg/h], (5) 
where: Dx [kg/h] — total boilers capacity, 
Nn [kW] — main propulsion power (main engine shaft power). 


In turn total auxiliary boilers capacity on container vessels is given in formula (6) 
described in [1]: 
D; = 0,0657 Na + 2536,6 [kg/h], (6) 


Dependency of total boilers capacity on main propulsion power for general cargo vessels 
and container vessels calculated according to formulas (5) and (6) is compared in figure 6. 
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Fig. 6. Comparison of total boilers capacity on general cargo vessels and container vessel 
as a function of main propulsion power 


27 


Much bigger heat energy consumption on container vessels (fig. 6) is the result of a 
considerably higher heavy fuel oil consumption. Heavy fuel oil system is the biggest heat 
energy consumer on board up to 70% of produced energy. On the other hand powerful 
engines of container vessels produce more exhaust gases which are used in gas heat boilers. 


5. Summary 


Obviously contemporary container vessels belong to ships with the biggest energetic plants 
i.e. main propulsion plant, electric power station and steam boilers. An example of these is the 
biggest container vessel in the world MS EMMA MAERSK built in 2006 with loading 
capacity 11,500 TEU and service speed 25 knots (fig.1). The main propulsion of the vessel is 
low speed diesel engine Wartsila Sulzer 14RT-flex 96C nominal shaft power 80,080 kW. 
Onboard power station consists of 5 diesel generators of total power 20,700 kW and 
additional 8,500 kW steam turbo alternator using steam from main engine exhaust gases heat 
auxiliary boiler. Thanks to high waste heat utilisation the energetic efficiency of engine room 
achieves 70% during sea passage. Economic analysis show that even bigger container vessels 
are expected. However there is a limitation in the maximum power of the main engine. The 
main engine mounted on MS EMMA MAERSK is the biggest diesel engine offered up till 
now by diesel engine producers. Today, the only alternative is double engine propulsion the 
same as used on modern biggest liquefied gas tankers. Other possibilities are turbine 
propulsion, COGES propulsion system and V-type slow speed diesel engines which have not 
been constructed yet. 
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Abstract 


The paper provides a proposal of a quantitative interpretation of operation which (as the operation of 
Hamilton and Maupertius presented in the classical mechanics and the operation after change of the body 
momentum) is considered as a physical quantity with the measurement unit called a joule-second 
[joulexsecond]. An original method for analyzing and estimating the engine operation has been demonstrated in 
energetic aspect for operational needs. Herein it has been shown that the operation of this kind of engines 
considered in the proposed aspect enables obtaining essential information about energetic properties of the 
engines, that completes the information regarding energy conversion in the form of work and heat. Possibilities 
of analyzing the diesel engine operation have been demonstrated in deterministic and probabilistic aspects. 
Basing on deterministic aspect of operation of this kind of engines there have been presented possibilities of 
determining the operational usability of the engines through determining the possible operation and demanded 
operation to perform a task. 


Keywords: operation, energy, technical state, diesel engine, Poisson process, semi-Markov process 


1. Introduction 


During operation of diesel engines we need to identify not only their technical states but 
also their energetic properties [5, 7, 8, 12, 13, 15]. The properties characterize the medium 
torque (Mo) and rotational speed of a crankshaft (n) of this type of engines. The torque M, and 
rotational speed n (as measurable values) enable to define a useful power (Ne) [12, 13, 15]. 
The useful power (No) is a quantity that characterizes the stream of energy converted in a form 
of useful work (Le) at a defined time (4). From this reason the work Le can be interpreted as 


the output of the delivered power N, = L, at time ¢ and therefore expressed with the formula: 
Loo N of (1) 


From the formula (1) results that the power M, is a quantity containing information how 
quickly the work L, has been (or can be) performed by a combustion engine. 

However in practice also a quantity is significant that provides information how long 
the work Le must be delivered by engine to a receiver (screw propeller of a ship, generator, 
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compressor) in order the given task could be performed. This quantity can be called operation 
[5, 6, 7]. Because each type of work being performed by engine (ex. useful work, 
compression, expansion, etc.) is a form of energy conversion, thus just understood operation 
(D) is a quantity expressing the energy (E) released over the time (£), and that is the reason it 
can be defined (when E=idem can be accepted) with the formula: 


D =Et. (2) 


The operation (2) determines thus the energy released over the time during which this 
energy has been consumed. When the engine wear is considered, the operation equals to the 
energy drop (decrease) at time at which it proceeded [4, 8, 11]. The energy can reveal only 
when converted into form of work or heat [2, 6, 8, 13, 14]. 

In case of any diesel engine the useful energy (Ee) generated by the engine with a 
defined useful power (Ne), in strictly determined conditions, can be considered as a measure 
of its ability to perform the work Le at a defined time ¢. Therefore the work as a form of 
energy conversion, generated by the engine, can be defined from the formula [13] 


L, =2mM_t, (3) 
in the case when: M, = idem and n = idem. 


When M, + idem and n + idem the work performed in the time interval [ż1, f] can be 
presented in the form of dependences: 


L= 2n| n(r)M (r)jdr. (4) 


ti 


In the operating practice of diesel engines (main engines) being applied to marine 
propulsion systems it is extremely important how long the work Le can be released for the 
needs of the propulsion system of the given ship. This refers especially to the ships of which 
propulsion systems are equipped with such engines. In the case when due to the wear, the 
main engine cannot be loaded with the demanded useful power (Ne) at time ¢, it is not able to 
perform in this time the demanded work Le needed to ensure generation of the demanded 
pressure force (T) by the screw propeller of the ship. In consequence the ship is not able to 
perform the transportation task. Moreover, when the cruise runs in storm conditions, it can 
lead to a catastrophe [9]. 

From the above considerations follows that it is reasonable to analyze not only the 
power N, released in diesel engine's workspaces, and simultaneously the work Le, but also the 
operation (D) of this type of engines, understood in this case as energy conversion in these 
workspaces that leads to obtaining the demanded useful work (Le) at a defined time (4). This 
will enable to fix whether the possible engine operation (Dm) for the given conditions is at 
least equal to the demanded operation (Dy) being indispensable to perform a defined task Z. 


2. Diesel engine operation as energy conversion in the form of heat and 
work 


Operation of engines consists in converting and transferring the supplied energy. In case 
of diesel engines, first the chemical energy contained in fuel-air mixture, generated in 
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workspaces, is converted into thermal energy and then the thermal energy — into mechanical 
energy [2, 5, 12, 13, 14]. 

It is obvious that the energy conversion in the form of heat in workspaces of each diesel 
engine can proceed at a different time. In practice it is essential to make the performance of 
the work as greatest as possible or as quickly as possible at a defined time. If it is not possible 
to obtain such energy conversion which is favorable the engine is considered to work 
incorrectly and to be in the state of partial usability [7, 12, 14]. 

In case of diesel engines, conversion of chemical energy into thermal energy and then 
into mechanical energy, enables creation of a torque (M,) of a crankshaft at a defined 
rotational speed (n) of each engine [12, 13]. Thus, the operation of engine, interpreted as 
energy conversion in form of useful work Le expressed with the formula (4) can be defined by 
the equation as follows 


D, = ft, (z)dr =2z [n(r)M, (rade (5) 


ti ti 


Engine operation connected with energy conversion in the form of work like 
compression of fresh charge, expansion of combustion gases in a cylinder, etc. can be 
considered in a similar way. 

Determination of engine operation consisting in conversion of chemical energy (Een) 
contained in fuel-air mixture generated in engine combustion chambers into thermal energy 
(E.) is equally important. Such operation (Fig. 1) when conversion of this kind of energy 
proceeds in the form of heat (Q) can be defined by the formula: 


D, = O()dr (6) 


Because the operation of this kind of engines consists in converting the energy £ in the 
form of work and heat, can be generally interpreted as follows 


D= | E(r)dr (7) 
where: 
D —engine operation, E — converted (obtained) energy enabling realization of a task Z, 
t —time of E energy conversion (consumption). 


Usability of particular combustion engines can be inferred after making value 
calculations of their operations (7) which are. in the interpretation proposed herein, equaled to 
physical quantities with the measurement unit: „joule-second”. Apparently, the functional 
dependence of energy from time, so E = ft) must be known in order to determine the field of 
operation (D). Because D = AE, t) the operation of machines can be presented in the 
coordinate system ,, D—E-?” [5, 6, 7]. 

Such understood operation defined by the formula (7) can be presented in the coordinate 
system ,,E—f” so in the form of graph which I propose to call operation graph. An example of 
such an operation graph for the range of the energy transformation from Ej into Ez for any 
selected time f = 0 and ¢ is presented in Fig. 1. 
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From the formulas (2) and (7) follows that the functional dependence of energy (£) 
from time (f) must be known in order to determine the operation field (D). Because D = KE, 
t), the machine operation can be displayed in the co-ordinate system „D, E, t”. 


operation field (D) 


Emax 


Fig. 1. An exemplary graph of engine operation: E[T /energy, Emax — maximum energy, 
E, — energy in the moment t, t [ fime 


Diesel engine operation can be and sometimes must be considered as a stochastic 
process [1, 3. 6, 7]. Such operation can be then displayed in a form of stochastic process 
realization as the dependence {E(f); t > 0], where energy £ is a random value. The process is 
characterized by the expected value E[E(f)] and the standard deviation o|E(t)] of energy £E. 
Such approach follows from that the analysis and the resulted estimation of the operation of 
diesel engines can be presented in a probabilistic aspect by applying the theory of the 
stochastic processes. An exemplary graph of such engine operation is show in Fig. 2. 


EIE + ol EO] 


E[E(9)] 


Fig. 2. An example of a stochastic process showing the dependence E(t), where E is a random value: E[[] 
energy, E, — energy assigned to time tı, E, — energy assigned to time t, t [ fime being a parameter of the 
process, E[E(t)] — expected value E, ofE(t)] — standard deviation of E 


A stochastic process is a random function of which the parameter is the time ¢. The time 
is not a random variable [1, 3]. This approach towards the issue of expressing the diesel 
engine operation as a value, results from the necessity of getting information what the 
operation can be in the interval defined by two arbitrary moments, ex. in the interval [fo, tn]. In 
this case, analyzing operation of each combustion engine, each time ¢ from the considered 
time interval [fo, ¢,)] can be assigned by a state called momentary state of the process, which 
is a random variable X, with the excepted value E(X;) and variation D(X) dependent from the 
t value. For the considerations the variable can be energy (£) or forms of its conversion, so 
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work (Le) or heat (Q). Thus the stochastic process (a random function) is a set of random 
variables X, for t € [fo, tn], so for to < t < tn. The function’s expected value E[X(A)] and 
variation D”[X(1)] are defined by the sets of expected values E(X) and variations D(X) for to 
< t < tn. It should be pointed here that the expected value E(X;) and the variation D(X) of the 
random function (X(t): £ > 0} depend on time ¢ because the values E(X) and D”(X)) can be 
different for different £ values. They are not, however, random functions X(t) because E(X) 
and D”(X,) are not random variables but the constants for the given £ value and the given set of 
realizations of the random variable X, [3]. 

An example of dependences of E(X) and D”(X,) from time ¢ is shown in Fig. 2, at the 
assumption that the random variable X is the energy E supplied by a combustion engine to a 
receiver of the energy. In this Fig. the o[£(£)] quantity is a standard deviation of the random 
variable Æ. This quantity is a square root of the variation D”(E)). 

In this case to define the operation D from the formula (7) the integral calculus can be 
applied because the integral defined by the formula is a definite Riemann integral with the 
integration range [0, £] and the integrand E(7). Because the function E(7) is continuous for the 
examined exemplary range [0, ¢], it can be stated in compliance with the second fundamental 
theorem of integral calculus (Newton-Leibniz Theorem) that [16] 


| EG)dr =D); (8) 


whereas: 
D(t) 


, = D(t)-D(0) 


Estimation of the expected value E(E,) for each value of time ¢ requires application of 
statistical inference, so point or interval estimation [1, 3]. 

From the presented interpretation of combustion engine operation follows that the 
operation consists in converting and transferring the released energy in the form of work (L) 
and heat (Q), whereas both of the forms of energy conversion can be presented as the fields 
[13, 14]: 

e in the Clapeyron diagram (diagram of work) when analyzing the work L (Fig. 3 and 
4), 
e inthe Belpaire diagram (diagram of heat) when analyzing the heat Q (Fig. 5). 


For instance, in case of a piston engine the absolute work (L4) of exhaust decompression 
(i.e. work determined in relation to the ambient pressure pọ = 0), of which the field is 
displayed in Fig. 3a in the p-V co-ordinate system, called the Clapeyron diagram, can be 
calculated from the formula 


L = Í pV)dV (9) 
where: p — pressure, V — volume. 
That means that the integral (9) can be determined when the functional dependence p = 


KV) is known. The technical work of exhaust decompression (L), of which the field is 
displayed in Fig. 3b can be determined from the formula 


L =-[V wap (10) 


Pi 
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Also in this case the integral (10) can be determined when the functional dependence V 
= fp) is known. The technical work of exhaust decompression can be however determined 
also by employing the definition of absolute work. Then the following formula is of 
application (Fig. 3) 


P2 
L, = pV, — pV,- |V(p)dp (11) 


Pi 


Area of absolute work (La) 


Area of technical work (Lì) 


Fig. 3. Graphical examples of work: a) absolute, b) technical: p — pressure, V — volume 


Taking into account that the work of exhaust decompression in the piston engine space 
is performed till the moment of opening the exhaust valve, when the exhaust gases are 
removed outside to the environment with the pressure p, = p, the useful work can be then 
considered (Fig. 4). 

The mentioned useful work (Z,,) can be defined (Fig. 4) from the formula 


L, = | pV)aV - p,(V,-V) (12) 


U 


where: p — pressure, / — volume, p, — barometric pressure (ambient pressure). 


Fig. 4. An exemplary graph of useful work: p — pressure, V — volume, p, — barometric pressure 
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A field of any other work (e.g. air compression in a cylinder, indicated work, useful 
work, etc.) can be presented in a similar way. 

The carried away heat can be shown (Fig. 5) in the form of a graph by using the ,,7— S” 
co-ordinate system (in the Belpaire diagram). The mentioned heat (Q) can be determined (Fig. 
5) from the formula: 


Sy 
Q= | T(S)dS (13) 
Sı 
where: T — absolute temperature, S — entropy. 


The following restrictions must be taken into account for the formula (13): 


T(t) 


Field of heat (Q) 


Fig. 5. An exemplary graph of heat: T — absolute temperature, S — entropy 


In the paper’s introduction it has been signaled that the operation in the interpretation 
presented herein will enable to determine if the possible engine operation (Dy) for the given 
conditions is at least equal to the demanded operation (Dy) being indispensable to perform the 
task (Z). That means the operation in the presented interpretation is of essential practical 
significance. 


3. Practical significance of engine operation with value interpretation 


Task for which a combustion engine has been designed and manufactured can be 
performed only when the following inequality is satisfied 


Dy Dy, (14) 
so when: 
tw =tw, when at the same time Ey > Ey. 


where: £y — possible operating time, ty — demanded operating time, Ey — energy that can be 


converted by engine, Ew — demanded (desired) energy to perform the task Z (energy that must 
be converted to enable performance of the task). 
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That means that when analyzing the energetic properties of combustion engines (not only 
diesel ones) ability of this type of engines (as well as other energetic systems) to work can be 
considered in the following alternatives: 

tm =tw, when simultameously Ey= Ey, 

ty =tw, when simultameously Ey> Ey, 

tm > tw, when simultameously Ey= Ey, 


tm > tw, when simultameously Ey> Ew. (15) 
In case when the inequality emerges: 
Dy< DW, (16) 


the engine is damaged and is not able to perform the task Z. 

The inequalities (14), (15) and (16) are also true when dissipation of energy converted 
in the form of heat is considered. In such a case the heat is carried away from the engine and 
therefore (in accordance with the interpretation used by thermodynamics) gets a negative 
value. 

Fig. 6. presents a case when the possible operation of engine Dy = E1 (t — tı) which 
while working can supply the energy E, = idem being indispensable to perform a defined task 
at time /,. In this case however, in order to perform the task the demanded operation is 
required to be greater than the possible one, so Dm < Dw, where Dy = E (t — tı). That means 
that the possible operation (Dy) of engine will not ensure the task performance and that is 
why before starting the task realization the engine should be submitted to reconditioning 
through performing the adequate preventive service work. 


| — Possible Operation Field (Dy) 
a + a — Demanded Operation Field (Dy) 


ti b b 


Fig. 6. An exemplary graph of possible operation (Dy) and demanded operation (Dy): 
E — energy, E; — energy assigned to time tı, E, — energy assigned to time t2, t time 


When no preventive service is carried out and any refurbishment on the machine is not 
done the performance of its operation can be presented in the form of a field provided in Fig. 
7. 


— Possible Operation Field (Dm) 
B + [] — Demanded Operation Field (Dy) 


ti to b 


Fig. 7. An exemplary graph of possible operation (Dy) and demanded operation (Dy): 
E — energy, E; — energy assigned to time tı, E — energy assigned to time t, t —time 
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When diesel engine energy transmitted to a receiver is constant (E = idem) in the time 
interval [¢1, 4], the operation follows in accordance with the formula (7) 


D= | E(p)dt = E| dt = Elp = E(t, -1) (17) 


Taking into account the ways of energy E conversion (7) in energetic machines, work 
(L) and heat (Q), their operation can be determined by using the formulas (5) and (6), as 
follows: 


D, = Lat = Lit, —4,); Do = |0at =Q(t, -t,) (18) 


ti 


Dependences (17) and (18) own interesting cognitive attributes, but may be also of 
utility significance. Applying them it is very easy to determine the demanded operation (Dw) 
as well as the possible operation (Dm) for each energetic machine and to obtain preliminary 
information on its usefulness (operational) to perform a defined task. Moreover, the formulas 
enable making simple graphs of demanded operation (Dw) and possible operation (Dm), which 
are presented in Fig. 6. 

The characterized operation of diesel engines and the examples of its demonstration in 
the form of operation fields refer to the case when important is to determine energy or its 
possible conversions in forms of work and/or heat being indispensable to ensure performance 
of the given task. However, each operation of an energetic system at determined time is 
followed by energy dissipation in accordance with the second law of thermodynamics. In case 
of diesel engines, a part of the produced energy is used to overcome their mechanical 
resistances. From this reason for the piston diesel engines we distinguish indicated work (L,) 
and connected with it indicated power (N;) and useful work (Le) and connected with it 
effective power (Ne). The indicated work of engine is a power being produced by the engine 
in its working spaces (in cylinders), without regard to its own mechanical resistances. The 
effective power, however, is a power which can be delivered to a power receiver in any 
conditions of engine operation. The difference between the two powers is the power (Nm) lost 
for overcoming the mechanical resistances of the engine. Taking this power into account, the 
energy or forms of its conversion, which are work and heat, can be considered from the 
conventional zero level. Then, the operation of a diesel engine can be presented with the same 
formulas but in the form of a field determined by the function E(f) or L(t) or Q(t) and the time 
axis (1), what is displayed in Fig. 8. for the case when change in energy at time ¢ is considered. 


Operation Field (D) 


Fig. 8. An exemplary graph of operation for the case when E = f(t): E — energy, E; — energy 
assigned to time t;, E, — energy assigned to time t, t —time 
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In case, when the energy average value is considered, so when it can be accepted that E 
= idem in the time interval [/, t2], the operation field can be presented as in Fig. 9. 


E Operation field (D) 


Fig. 9. An exemplary graph of operation for the case when E = idem: E — energy, E; — energy assigned to time 
tı, E;— energy assigned to time tz, t —time 


Dependences (17) and (18) follow from that the integral expressed with the formula (7) 
or (5) or (6) is the Riemman integral with the interval [ż, £] in this case and the sub-interval 
function E(t) = E. This function is integrable in terms of Riemman in the mentioned time 
interval according to the formula: 


| Edt = Bt, -1)=(E, -Ea h) (19) 


q 


However, in practice it is not always possible to accept that energy supplied by an 
energetic machine to a receiver is constant [8, 12, 17, 18, 19]. Then we need to define the 
functional dependence of energy (E) from the time (f) of machine operation, so the 
dependence E = f(t). In case when this is possible because the function E = f(t) is continuous 
in the considered exemplary interval [t1, t2], following the second fundamental theorem of 
calculus (theorem by Newton and Leibniz) we can write in accordance with the formula (8) 
that (Fig. 8). 


eo = D(t,) - D(t,) (20) 


fi 


Application of the theorem by Newton and Leibniz is here necessary because it enables 
effective determination of the definite integral of the continuous functions if determination of 
any primitive for the functions is possible. 

Generally, the functional dependence E = ff) is composite. In case, when for such a 
function the internal function derivative is a constant function, the definite integral of the 
function f(t) can be determined by application of integration consisting in substitution [22]. 

Not always however the elementary function describing the dependence of energy from 
time is possible to be defined with elementary functions. Then determination of the definite 
integral from the Newton-Leibniz formula is troublesome, and sometimes even impossible. 
The trouble is that determination of a primitive is connected then with necessity of making 
difficult transformations. In such cases, just like when sub-integral function is determined in a 
table form, we can calculate an approximate value of the operation of an energetic machine, 
being the definite integral value, applying the method of trapezoids or the Simpson method. 
The attention must be paid that at a fixed n the Simpson method enables obtaining more 
accurate results in integration than the method of trapezoids. 
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9. Conclusion 


Operation of combustion engines is understood as generating the energy E by them at 
a defined time t. It has been equated to a physical quantity which can be expressed with a 
numerical value and the measurement unit called joule-second [joulexsecond]. Such 
understood operation gets worse with the growing wear of this type of engines. This means 
that the operation value at a defined time decreases in the result of decreasing energy 
generated by the engines. It has been signaled herein that in case of application of the theory 
of stochastic processes to the analysis of changes of such understood operation the integral 
calculus can be applied. Two stochastic models of decreasing useful energy generated by the 
engines have been proposed for defining the range of worsening operation. The first model 
has been presented in the form of a homogenous Poisson process and the second — in the form 
of a discrete-state, continuous-time semi-Markov process. 

Operation in such interpretation depends on the technical state of the engines and is 
characterized simultaneously by the energy converted by the engines and the energy 
generation time. 

The advantage of the engine operation in the presented interpretation is that it can be 
tested through doing precise measurements and then expressed in the form of: 

— a number with the measurement unit called joule-second [joulexsecond] (formulas 5, 
6 and 7); 

— a graph, as a field of operation (Fig. 2, 3 and 4). 

Operation in the presented interpretation, although formulated for diesel engines, refers 
also to spark-ignition engines. Similar interpretation of operation can be provided for turbine 
combustion engines and other energetic machines. 
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Abstract 


In this article the vertical mixed flow pumps and problems accompanying during working of the devices vere 
presented. At the begining few issues which are consider in this type of machines were mentioned. The results of 
measurements a velocity of vibrations on the pump were presented in next point of the paper and compared with norm. 
Important issue in the article is modeling, therefore 3D model, discreet model of the rotor, numerical analysis of the 
rotor, four-masses discreet model and equations of the free motion of the rotor were presented. For finding stiffness 
and suppression parameters of bearings the special research stand were built and shown in the article. Range of 
future works was also described. 


Keywords: vibrations, dynamic analysis, impeller pump, mixed flow pump, structural model 
1. Introduction 


The object considered in this work is mixed flow pump with two impellers, which works in 
company from Włocławek. Large dimensions, complicated structure and mechanical processes 
which accompany during working of the device provide too difficult dynamic analysis of the 
machine. Additionaly is not possible to measure vibrations on part of the device which is plunged 
in the water. Rotary motion of the mojor shaft with two impellers is a basic motion in the pump. In 
these devices, a special attention is directed on: 

- stability of rotor and critical rotations, 

- unbalancing of blades and the rotor, 

-defining the level of dynamic loading of bearings and the support structure, 
- identyfication loads generated during working of the machine [2]. 


2. Specifity of research object 


High power - 1250kW and efficiency - 5000 m°/h are characteristic for the device. 
These parameters have influence on the way of supporting the pump, the device is spupported on 
two foundations. Basic parts of the machine: electrical motor which is based on upper foundation 
(fig. 1.) and pump — supported by bottom foundation (fig. 2.). 
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Fig.1. Electrical motor of mixed flow Fig.2. Part of machine to escape the 
pump water 


In one hall were installed 8 same pumps, because demand of water is different. The measure of 
vibrations is difficult when few devices are working together. Ratory speed of the pump is 
established on 740 rpm. 


3. Results of measures 


Measures of vibrations were based on norm PN — ISO 10816-1:1998, which describe border 
levels of intensity of vibrations for classified devices (tab. 1). The considered pump is classified as 
II class, but intesity of vibrations in B zone (for devices admitted to long exploitation). 


Tab. 1. Border values of intensity of vibrations according to PN-ISO 10816-1:1998 


Velocity of vibrations | Class | Class | Class | Class 
RMS in mm/s I II III IV 
0,28 
0,45 A 
0,71 à A A 
1,12 
B 
1,8 
2,8 j 
45 E B 
; [e B 
7,1 
[e 
11,2 c 
18 D D 
28 D D 
45 


Measure points of vibrations are shown on figure 3, results of measures in table 2. 


42 


Tab. 2. Results of measures 


Direction of measures vibrations 
Measure x y 

point Veloc. | Accel. | Veloc. | Accel. | Veloc. | Accel. 
[mm/s] | [m/s] | [mm/s] | [m/s”] | [mm/s] | [m/s] 

1 4,0108 | 1,2332 | 5,5884 | 1,8884 - - 

2 5,0412 | 1,3568 | 5,158 | 1,6092 - - 

3 1,394 | 1,2896 | 1,5112 | 1,7412 - - 

4 1,276 0,61 1,1312 | 0,7532 - - 
5 š Š - - 2,2176 | 0,7012 

6 0,6392 | 0,7288 | 1,3912 | 0,6996 - - 

7 1,1968 | 0,4456 | 0,7972 | 0,7984 - - 
8 - - - - 0,4552 | 0,5628 


Fig. 3. Measure pionts of vibrations 


4. Modelling 


For creating spatial model (fig. 4.) of the pump, documentation 2D were used. The model was 
useful to numerical analysis and wizualization how the device works. Additionally, the model is 
used to fast estimate physical properties of elements like for example: moments of inetria or 


centers of masses. 


Fig.4. Cutting view of vertical mixed flow pump 


Basic part of the pump is consisting from three rotors, which are connected by stiff clutches and 
two impellers. Whole part was numerical analyzed in MSC.ADAMS program. Parameters of 
stiffness and suppression were taken from literature: k=1-10° N/mm, c=INs/mm. Results of the 


analysis are shown on figure 5. 
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Rys. 5. Numerical analysis in MSC.ADAMS program 


Results of the numerical analysis and energetical Rayleigh’s method were used to creation four- 
masses discreet model of rotor stiffness supported (fig. 6.). 


Fig. 6. Discreet model of rotor 


For the model equations of the free motion were written (1): 


ô mM, Yı + ÔM, Y, + ÔM V3 + ÔM, y+ y =0 


01M, Yı + ÔM, Y+ 33M; Y3 + Oy4M, Y4 +y, =0 (1) 


31M, Yı + ÓW, Y, + Oy, Y3 + 04M, Y4* y; =0 


Ó4M, Yı+ pM, Y,+Ó GM, Y;+0yM, Yat Yy, =O 
where: 


Ojj — Maxwell’s coefficients, 
m; — discreet masses, 
yi — statical deformation of rotor, 


and solved. At present the model is werryfied. 


The discreet model is simply in relation with real object. Therefore, in next step the parameters of 
stiffness and suppression of supports will be regarded. Parameters of k and c were taken from 
literature. Verification of the parameters will be done in experiments. Price of real bearings in high, 
therefore research stand for bearing with different diameter, clearance between toe and bearing, 
pressure of reinforced water, rotary speed of rotor was built (fig. 7.). 
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Fig. 6. Research stand 


For different diameters of bearings, it is planed to find formula which will be use for finding k and 
c parameters. 


5. Range of future works 


Next works will concern the following topisc: 

- creation of mathematical model which will nclude stiffness and suppression parameters, 

- identyfication of bearings parameters, 

- identyfication loads generated during working of the machine, 

- reserching of stability of rotor and critical rotations, 

- reserching of inluences between different places of supports and durability of plugging and 
bearing, 

- reserching of inluences kinematical random impacts from trunk and rotor’s vibrations, 

- describing guidelines for construction changing. 


6. Conclusions 


1. Problems with veryfication of ceated model exist, because is not possible to measure vibrations 
on all parts of the object. 

2. Measures of velocity of vibrations showed that them (these, their?) values are higher than those 
described in norm. 

3. Created discreet model will be developed by consider parameters based on reserches. 

4. Beginning calculations showed that torsional vibrations can be skipped. 
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Abstract 


The fundamental problem with wide application of the Wankel engine is to assure oil film 
continuity on trochoidal cylinder bearing surface. For the sake of considerable difference 
between curvature radius of the apex seal sliding surface and the trochoidal cylinder the oil film 
can be generated on the short fragment of the apex seal. However there is certain field of 
maneuver in the apex seal shape determination which allows to approach close radiuses of 
curvature of both elements in the areas where the highest gas forces load is occurred. In the paper 
simulation of the revised apex seal oil film parameters are presented. Author described also the 
test stand which is going to be used for experimental verification of the simulation results. 
Conclusions refer to possibility of replacement the constant radius shape of the apex seal sliding 
surface with shape that consist of two different curvatures. 


Keywords: the Wankel engine, oil film parameters 


1. Introduction 


Problems with the apex seal of the Wankel engine has been known since the first engine of that 
type was designed. The first apex seal designs where just to make start of the engine possible, but 
in the later engine development apex seal was limiting the proper working time of the engine. In 
early prototypes usually after just few hour of working the blow-by between adjacent working 
chambers were enough to stop the engine. Later rotary engine constructions which were applied in 
serial produced NSU Ro80 cars or even the Mazda Renesis engine apex seal was not the main 
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disadvantage but still it is less effective, less durable and generates more friction losses than 
modern piston rings in conventional reciprocating engines. Previously described difficulties can be 
easily explained if the apex seal angle of attack is taken into consideration. The angle between the 
apex seal axis and the normal to the trochoidal cylinder liner varies within limit of one radian 
while in conventional piston engine analogous angle do not exceed 0,001 rad. As a consequence in 
conventional engine the maximum gap thickness between piston ring and cylinder liner do not cut 
across few micrometers whereas in the Wankel engine gap reach size of several hundred 
micrometers. There is a question if the apex seal sections with circular or parabolic profile, that are 
being used nowadays, provide the best results both of the oil film thickness and friction losses 
caused by apex seal sliding on the trochoidal cylinder surface. As it is pointed out in present 
elaboration it is possible to propose the apex seal sliding surface shapes which guarantee 
acceptable oil film parameters when the apex seal is situated perpendicularly to the trochoid 
surface. What is more the oil film parameters can not get extremely worsen during inclination of 
the apex seal, mainly in areas where the long and short axis of the trochoid are evenly distant. 
When looking for analogy in functions of the apex seal and the piston ring it can be assumed that 
in the reciprocating engine the piston ring could cooperate with cylinder liner in similarly hard 
conditions as apex seal of the Wankel engine if the cylinder liner would change its diameter by 
few millimeter. The generating line of this cylinder liner would be sinusoidal shaped with 
amplitude of few millimeters and period comparable to piston stroke. This kind of model was 
engaged in computer simulations during the new profile of the apex seal examinations. 


2. Parameters of traditional cylinder liner that creates conditions of the apex seal [1] 
On the figure 1 changes course of the angle between tangent to the piston ring and tangent to 


the cylinder liner which is deformed in such manner that generating line of the cylinder liner has 
sinusoidal shape with amplitude of 1,6mm and period of 160mm. 


Fig. 1. The angle course changes of tangent to the cylinder liner generating line, the amplitude is 1,6mm and the 
period is 160mm 


The tangent line to the sliding surface of the classical piston ring forms with piston axis lesser 
angle than angles shown on fig. 1. As a consequence the piston ring contacts with cylinder liner in 
predominant part of the working cycle alternately with its upper and lower edge and there is no 
possibility to generate continuous oil film between the piston ring and the cylinder liner. This 
situation is very similar to that one which occurs in the Wankel engine, where the apex seal also at 
considerable part of its path works with its edge. In apex seal that have been used so far the sliding 
surface was a parabolic shape. With the purpose of improving the oil film parameters the apex 
seal shape that consist of two parabolas was considered. 50% of the profile keep normal shape and 
the other part, closer to the edges, is multiplied. As a result the apex seal shape looks as it is shown 
on upper right corner of fig. 2. 
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When the corrected sliding surface was used oil film parameters presented on fig. 2 it the part 
signed as “Pier. Zg.” were obtained. If the next apex seal would not be corrected, so its profile is 
formed by one parabola, the oil film on the most part of apex seal path does not exist. It can be 
observed on fig. 2 part that is signed as “ Piersc. 2”. The not corrected profile shape is also 
presented in the middle part of fig. 2 above the diagrams. 


WANKEL 0.200mm Pierscen 2 0.200mm Pierse. zędr. 
Grubose warstwy i filmu olejowego 
do piers. 2 ( = ) 
i dla pierse. zg. ( ) 
lodzi tule’ d 
na glodzi tulei cylindrowej im 280 


um 
22 


O Warstwa 180 Pier. zg. 180 Piersc. 2 180 Film 180 
Predk. obr om — 92 [rad/s] Prom. korby ra- 47 [mm] Kod — 1 / m 
Lepkosc ol. Eta — 139 [mPas] Nacisk pier. Ps2 — 0.20 [MPa] Pso — 0.20 [MPa] 


Fig. 2. Comparison of the oil layer thickness course that is left on the cylinder liner surface by two profiles of the 
apex seal — first part of the diagram; oil film thickness and layer thickness for double parabolic profile — second part 
of the diagram; oil film thickness course and layer thickness course for the parabolic profile with very high 
curvature radius — third part of the diagram; comparison of the oil film thickness that was generated by the one and 
double parabolic profiles — fourth part of the diagram 


The sliding surface correction of the apex seal has to be done very carefully because 
improving oil film parameters in the key areas leads to deterioration of this parameters in other 
areas. On fig. 3 synthetic course of the oil film parameters generated by two following apex seals 
of the Wankel engine. 


WANKEL 0.200mm Pierscien 2 0.200mm Pierse. zgar, 
Grubose warstwy i filmu olejowego 


© Warstwa 180 Pier. zg. 180 Piersc. 2 180 Film 180 
Predk. obr. om = 92 [rad/s] Prom. korby ra — 47 [mm] Kod -5/m 
Lepkosc ol. Eta = 139 [mPas Naclsk pler. Ps2 — 0.20 [MPa] Pso — 0.20 [MPa] 


Fig. 3. Comparison of oil film thickness course which is lest on the cylinder surface by two apex seal profiles, for 
conventional apex seals of the Wankel engine — first part of the diagram; oil film and layer thicknesses for one 
parabolic profile — second part of the diagram; oil film and layer thicknesses for parabolic with very high curvature 
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radius — third part of the diagram; comparison of oil film thickness generated by two one parabolic apex seals — 
fourth part of the diagram 


The results obtained from the simulations are classical course of the oil film parameters for types 
of apex seal that have been used so far. It can be noticed on the first part of the fig. 3 diagram that 
in the 90 degrees to 180 degrees area the classical one parabolic apex seal leaves no oil layer on 
the cylinder surface, which means that there is no oil film and intense wear of the trochoidal 
cylinder surface occurs. However comparing the maximum oil film thicknesses it can be observed 
that for classical apex seal this parameter is better because it reaches value of 33 um fig. 3, while 
for the corrected two parabolic profiles only 25 um — fig. 2. 

The maximum oil film thickness areas for the classical apex seal provide to low friction forces 
values. It is possible that mechanical efficiency of the Wankel engine with regular apex seals is 
higher than with corrected two parabolic solutions. Simulations that have been used do not allow 
to unequivocally confirm this assumption. The partially lack of oil film on the trochoid surface 
needs other ways of friction losses estimation. The best way to estimate friction losses is to use 
one of the experimental method what is in realization phase with authors cooperation. In 
continuous oil film case the computer simulations allow to specify friction forces and example of 
its results is presented on fig. 4. 


WANKEL — Grubose filmu, warstwy i tarcie ZZM BA a 


Kod — 5 m om — 92 [rad/s] ps3 — 0.20 [MPa] 
eta — 139.00 [mPas] ms — 1.00 [mrad] ra — 47.00 [mm] 
po — 0.10 [MPa] pd — 0.10 [MPa] pw — 0.11 [MPa] 
eps — 18.50 T = TAS fos 110 
hr — 2.00 [mm] du — 82.00 [mm] bet — 1.10 i 
Nr — 46.23 [W] Amp — —0.80 [mm]Tok — 160.00 [mm] ew 
[mikrom] [N 
25.00 20.0 
18.75 15.0 
12.50 10.0 
6.25 5.0 
0.00 0.0 
ie) 45 90 .. 135 [owk] 
fim — = warstwa — tarcie — — 


Fig. 4. The oil film and layer thicknesses for one parabolic profile and friction force course for the same apex seal 
profile 


Very often next to basic oil film parameters which are the film thickness, layer thickness and 
friction force also piston profile filling by the oil film is determined. For analyzed cases profile 
filling by the oil film is illustrated on fig. 5 and fig. 6. In two parabolic corrected profile case the 
area where the oil film is generating by the parabola near the edge of slat is evident. It is area of 
90 to 180 degrees range — fig. 5. For the typical apex seal with one parabolic profile filling by the 
oil film in 100 to 150 degrees range is lower so it is probable that the friction force is also lower — 
fig. 6. 

In so far considerations because of its comparative nature the course of angle between tangent 
to the apex seal and tangent to the trochoid, which is presented on fig. 1, was taken into account. In 
actual conceptions of the Wankel engine the course of attack angle can be different and it depends 
on the Z parameter. On figure 7 the angle of attack course for greater value of the Z parameter. 

In order to determine the initial oil film parameters which are obtained for matching the 
trochoid with two parabolic profile and for inclination angles shown on fig. 7 the suitable 
computer simulations were realized and it results are put together on fig. 8 and fig. 9. 
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WANKEL — Pokrycie piersc. zgarn. olejem o.20mm __Pler. zgi 


Kod — 1 m om — 82 [rad/s] ps3 — 0.20 [MPa] 
eta — 139.00 [mPas] ms — 1.00 [mrad] ra — 47.00 [mm] 
po — 0.10 [MPa] pd — 0.10 [MPa] pw — 4.11 [MPa] 
eps — 18.50 fi - 110 ro - 1.10 20.00 
hr — 2.00 [mm] du — 82.00 [mm] bet — 1.10 um 
Nr — 75.42 [W] Amp — -1.60 [mm] Tok — 160.00 [mm] 


ZW 


[mm] 


2.00 


1.50 


1.00 


0 45 90 35 [owk] 
Fig. 5. Course of two parabolic profile filling by the oil film 


WANKEL — Pokrycie piersc. zgarn. olejem o.20mm_Pler. zg. 


Kod - 5m om — 92 [rad/s] ps3 — 0.20 [MPa] c: 
eta — 139.00 [mPas] ms — 1.00 [mrad] ra — 47.00 [mm] J 
po — 0.10 [MPa] pd — 0.10 [MPa] pw — 0.11 [MPa] 
eps — 18.50 fi — 110 ra — 1.10 20.00 
hr — 2.00 [mm] du — 82.00 [mm] bet — 1.10 um = e 
Nr — 46.23 [W] Amp — —0.80 [mm]Tok — 160.00 [mm]?2 zw 
[mm] 
2.00 = 
| h 
1250. 
H AL 
1.00 ka 
+4 tH 
0.50 || | | 
o.oo III fl LINII! W 
0 45 90 35 [Owk] 


Fig. 6. Course of parabolic profile filling by the oil film, a high curvature radius profile 


[mrad] Kod = 5 / m 


0 80 360 540 [Owk] 
Fig. 7. The course of angle changes of tangent to the cylinder liner generating line, the cylinder liner deformation is 
sinusoidal with the amplitude of 0,8mm and the period of 160mm 
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WANKEL 0.200mm Pierscien 2 0.200mmm Pierze. zgdr. 


Grubose warstwy I filmu olejowego 

dla piese. 2 ( ) 

1 dla plerse, ( ) 

na gladzi tulei cyindrowej sc z." 
um = um 
zz 


O Warstwa 180 Pier. zg. 180 Piersc. 2 180 Film 180 
Predk. obr. om 92 [rad/s] Prom. korby ra 47 [mm] Kod - 2 / m 
Lepkose ol. Eta — 139 [mPas Nacisk pier. Ps2 — 0.20 [MPo] Pso = 0.20 [MPa] 


Fig. 8. Comparison of the oil layer thickness course that is left on the cylinder liner surface by two profiles of the 
Wankel engine apex seal with high value of the Z parameter — first part of the diagram; oil film thickness and layer 
thickness courses for two parabolic profile — second part of the diagram; oil film thickness course and layer 
thickness course for the parabolic profile with very high curvature radius — third part of the diagram; comparison of 
the oil film thickness that was generated by the one and two parabolic profiles — fourth part of the diagram 


WANKEL 0.200mm Plersclen 2 0.200mm Plersc. zgar. 
Grubosc warstwy i filmu olejowega $ í j 
dla piers. 2 ( —_—< ) 
1 dla piersc. zg. ( ) 
na gladzi tulei cyindrowej 


25,04 
um 


Zw 


0.0 


O Warstwa 180 Pier. zg. 180 Piersc. 2 180 Film 180 
Predk. obr. om = 92 [rad/s] Prom. korby ra — 47 [mm] Kod = 3 / m 
Lepkose ol. Eta — 139 [mPas Nacisk pier. Ps2 — 0.20 [MPa] Pso — 0.20 [MPa] 


Fig. 9. Comparison of the oil layer thickness course that is left on the cylinder liner surface by two profiles of the 
Wankel engine apex seal with high value of the Z parameter — first part of the diagram; oil film thickness and layer 
thickness courses for two parabolic profile — second part of the diagram; oil film thickness course and layer thickness 
course for the parabolic profile with high curvature radius — third part of the diagram; comparison of the oil film 
thickness that was generated by the one parabolic profile and two parabolic profile which curvature radius is four 
times higher than on fig. 8 — fourth part of the diagram 


The two parabolic profile was accepted for calculations, smaller curvature radius — fig. 8 and four 
times higher radius of curvature — fig. 9. It is easily to notice that in second case the results are 
better. The oil film in second part of the diagram demonstrated on fig. 9 is thicker and what is 
more important there is no break of the oil film what can be noticed on fig. 8 between 150 and 180 
degrees. 
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3. The test stand used to verification of the computer simulations results 


To the investigations authors have adopted the test stand which was formerly used for oil film 
investigations of piston rings of the marine two stroke engines. Scheme of the test stand is 
presented on fig. 10 [2]. 


Fig. 10. Schema of construction of the research stand (description in text) [2] 


On presented stand reciprocating movement is executed by cart with settled plate simulating 
cylinder sliding surface 1. This cart is driven by hydraulic servo-motor 2. Pressures in working 
chambers of servo-motor are regulated by electromagnetic valve 3 which uses the difference of 
voltage between signals from steering apparatus 4 and sliding resistor which is connected with 
piston rod. The element that simulates piston ring 5 is seated in piston block model 6 which is 
settled in the saddle 8 with force transducer 7. 

Because in the apex seal motion there are no dead centers and the plate has finite length, the to- 
and-fro motion had to be substituted by a linear translational motion in the plate length limit. The 
next step was to reflect the apex seal kinematics was to replace analog steering method by its 
numerical equivalent. This allowed to specify any plate path and velocity profile, also the apex 
seal velocity profile. The additional advantage of numerical way of plate motion programming is 
that we can choose any part of apex seal track on trochoid cylinder liner to be simulated. The plate 
path determination consists in specifying its position for every of 4096 points. In case of 
reciprocating engine 4096 point mean full crankshaft turn and for the apex seal motion simulation 
this number is discrete division of any rotor rotation angle that has been chosen to be investigated. 
It was said that there are no death centers in the apex seal motion so it was important to make the 
acceleration and deceleration safe. Additional procedures were added to the main program which 
are responsible for initial and final phase realization of the plat motion. Initially authors assumed 
that with total plate length of 800mm only 500mm will be used to the apex seal motion simulation 
while 100mm will serve to acceleration and deceleration. This results from the necessity of 
minimization the inertia force that is generated during plate acceleration and deceleration. Actual 
view of the test stand is presented on fig. 11. 
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Fig. 11. View of the test stand, upper left corner — hydraulic drive, upper right corner — saddle with piston model 
attached, lower left corner — piston model with attached slat, lower right corner — the plate that simulates the cylinder 
liner surface 


4. Conclusions 


1. One of possible ways to improve the oil film parameters between the apex seal and trochoid 
sliding surface is to shape the apex seal profile as a two parabolic curve. 

2. In most of the apex seal path, where the apex seal axis is perpendicular or almost perpendicular 
to the trochoid sliding surface, the oil film is generated by high curvature radius parabola while 
in the vicinity of shorter axis of trochoid the oil film is sustained by parabolic profile with lower 
radius of curvature. 

3. The low radius of curvature profile part should be formed in such way that the convergent gap 
would not be less than 25% of total gap length between the apex seal and the trochoidal sliding 
surface. 
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Abstract 


The research in presented paper has shown that the physical aspects of interfacial phenomena, described by the 
total value of surface free energy and the values of its components, make it possible to select more suitable materials 
for sliding pairs. The total value of surface free energy depends on the molecular structure and the bonds characteris- 
tic of a given material and determines its hardness. In order to reduce friction losses in a sliding pair which is being 
designed, it is proposed to match such materials for the pair that the surface of one of them has a high sum of surface 
free energy components originating from van der Waals interactions while the other material 's surface has a possibly 
low value of the sum. Furthermore, proper values of the components of surface free energy ensure proper wettability 
with lubricating oil. Pursuing the practical goal of this research, a new piston packing ring/combustion engine PRC 
unit cylinder liner sliding pair was designed and made. 


Keywords: Surface Energy, Wear, Coatings, Nitriding, Diesel Engines 
1. Introduction 


Interfacial forces produced by intermolecular interactions occur on the surfaces of sliding pair 
components in both rotary motion and to-and-fro motion. In solid/liquid and liquid/gas interfacial 
regions atoms belonging to each of the phases are subject to a different system of forces than at- 
oms located within the phases [1]. Atoms at an interface are attracted by both their own phase at- 
oms and the adjoining phase atoms whereby they are located in an asymmetric field of forces. 
When the forces of attraction towards one of the phases are greater, the atoms migrate into the 
phase until they reach an equilibrium through changes in the distances between the atoms located 
at the interface [5]. The differentiation of forces in the interfacial region is the cause of many phe- 
nomena such as: adsorption, wetting, adhesion and so on. 

The state of phase boundaries can be described on the basis of thermodynamics [6]. Various 
thermodynamic functions such as: internal energy (U), free energy (F) and free enthalpy (G) (also 
called the Gibbs function) are used for this purpose. In diphase systems one can distinguish a re- 
gion of pure phases, referred to as a and , and a region of transient surface phase o between the 
boundary phases [7]. Therefore free enthalpy G (at a constant temperature and pressure) or free 
energy F (at a constant temperature and volume) is used describe surface properties. Hence inter- 
facial tension y is defined as: 
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rE) = (E a) 
dA |, dA ),, 


were: 


y — surface free energy, 
F — free energy, 

4 — surface, 

G — free enthalpy. 


Interfacial tension y is also referred to as interfacial free energy or in short, interfacial energy 
(especially in the case of solids). When air is one of the phases or the considered phase adjoins a 
vacuum, then instead of the attribute ‘interfacial’, the attribute ‘surface’ is used. In his latest com- 
prehensive study [18] H. Lyklema recommends to use the term interfacial (surface) tension for 
both liquids and solids. Sometimes instead of y, symbol F’ , referred to as excess surface free en- 


ergy, is used. In the present work, surface free energy y and surface tension are used for respec- 
tively solids and liquids. 


2. Calculation of surface free energy from surface tension of molten solid body 


It is difficult to determine surface free energy for solids. According to many researchers, it is 
possible to calculate surface free energy from measurements of the free surface of melted solids 
[18]. Such problems feature prominently in metallurgy where the surface tension of fused metals is 
of major importance [3]. Allen [2] compared surface tension values for metals at their melting 
points and found a clear relationship between fusion temperature and surface tension for different 
metals. [5]. The relationship between surface tension and temperature is expressed by the Gibbs- 
Helmholtz equation: 


d 
Ys =¥m+—(T-T,), (2) 


where: 


ys — surface tension of molten metal at temperature 7, 
Ym — Surface tension of molten metal at fusion temperature Tn, 


T,„ — the fusion temperature of the metal, 
d . . . ; 
os — Temperature coefficient, whose value is always negative for metals, was introduced [1]. 


3. Calculation of surface free energy from Young’s modulus and other parameters 


R.J. Good [10] proposed a relationship between surface free energy and Young’s modulus, as- 
suming that the value of Young’s modulus EF is a consequence of intermolecular force F»: 


dF. 
E= at dr | (3) 


where: 
E — Young’s modulus, 


ro — a distance between the centres of molecules being in equilibrium between attractive and 
repulsive forces, 
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F, — the force of intermolecular interactions at distance r. 


Expression (3) holds for intermolecular distance r equal to ro and for constant temperature. 
From R.J. Good’s considerations [10], based on the Lennard-Jones equation [17] for the energy of 
a system of molecules versus the distance between them, the following expression interrelating 
surface free energy with Young’s modulus was derived [10]: 


327, Er, 
=—— > =—, 4 
a Ys 32 (4) 


where: Young’s modulus E and distance rọ are in respectively pascals (Pa) and meters (m) and 
calculated surface free energy ys of a solid is in joules per square meter (J/m*) or newtons 
per meter (N/m). 


Expression (4) allows one to estimate surface free energy on the basis of Young’s modulus of 
a solid body and the intermolecular distance at which intermolecular interactions are at equilib- 
rium. Since the material property resulting from intermolecular interactions is taken into account 
through expression (4) it is possible to calculate the total surface free energy if one knows the kind 
and number of bonds in a given solid. Since the surface free energy calculated on the basis of 
Young’s modulus is for a material devoid of lattice defects, it may be considerably overestimated. 

It follows from the above considerations that surface free energy is a result of interactions 
characteristic of the molecular structure of a given phase. Another material property is hardness 
and similarly as surface free energy it depends on a given material’s structure. Material hardness is 
resistance to the forcing of elastic and plastic deformations. Therefore it is associated with volume 
elasticity modulus K which is a measurable and, as opposed to hardness [23], well defined quan- 
tity. The hardness of selected materials versus their volume elasticity modulus K is shown in fig. 1. 
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Fig.l. Hardness HV versus K for hard materials 


By approximating the relation shown in fig. 1 with a linear function the following expres- 
sion was obtained: 
HV =a. K =10.236-K, 16) 
where: 


HV — Vickers hardness, dimensionless quantity, 
K — the modulus of volume elasticity of a perfect crystal [GPa], 
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a — a constant equal to 10.236 [1/GPa]. 


According to A. Witek [22], the value of the volume elasticity modulus is consistent with the 
Marvin Cohen (Berkeley University) theory proposed in 1985 and expanded in 1987. The theory 
allows one to calculate (with an error below 2%) the modulus of volume elasticity from this equa- 
tion: 


K=1761:r;**, (6) 
where: 


K — the modulus of volume elasticity of a perfect crystal [GPa], 
ro — an equilibrium distance between atoms or ions, i.e. the length of a chemical bond in the crystal 
lattice [angstroms]. 


According to the Cohen theory, the modulus of volume elasticity can be calculated from 
the bond length. Volume elasticity coefficient K is related to Young’s modulus E by the following 
formula [24]: 


K =——_~ (7) 
where: 


K — a modulus of volume elasticity [Pa], 
E — Young's modulus [Pa], 
v — a Poisson fraction of 0-0.5. 


Young’s modulus in turn is related through relations defined by R.J. Good [10] (described by 
formula 15) to surface free energy. Hence total surface free energy and hardness are proportionally 
correlated through the modulus of volume elasticity and the Poisson fraction as follows: 


Er, 1-2v 
= = 3Kr, ; 8 
Ys 32 USZY) (8) 


where: 


the meaning and the units of all the variables are consistent with the denotations in equations (4) 
and (7). 


Using equations (5), (6) and (8) the following expression interrelating surface free energy ys, 
hardness HV, intermolecular distance rp and Poisson fraction v was obtained: 


1-2v 
Vs SE a S (9) 


where: 
ys — surface free energy [mN/m], 
HV — Vickers hardness, dimensionless quantity, 


rg — a mean intermolecular distance [nm], 
v — a Poisson ratio of 0-0.5. 
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Equations (8) and (9) were formulated for a perfect crystal free of lattice defects. Their appli- 
cation is, however, limited since they yield, similarly as equation (23), overestimated surface free 
energy values. 


3. Calculation of surface free energy from its components 


F.M. Fowkes [8] proposed to divide surface free energy into the following components: 


yay +y" ty +7, (10) 


where: 


Y —a dispersion component, 

7” — a metallic bonds component, 
h 

Y — a hydrogen bonds component, 


7 — other components, e.g. an ionic component, a covalent bonds component, etc. 


In a simplified formulation, the surface free energy of liquids consists of a dispersion part and a 
polar part. For example, the surface free energy of water is 72.8 mN/m and it consists of dispersion 
part y” = 21.8 mN/m and polar part y” = 51 mN/m. Water has neither a metallic component nor a 
covalent component. 

There are several methods of determining surface free energy components, based on solid body 
wetting phenomena. They employ two or three calibrating liquids with known surface free energy 
components and assume that surface free energy has only two components: a dispersion compo- 
nent and a polar component [13]. The latter is composed of Keesom (permanent dipoles) and De- 
bye interactions (inductive dipoles) of the dipole-induced dipole type, of confirmative and multi- 
pole interactions and hydrogen bonds: 


yL=yi +y? and y,=75 +78, (11) 
where: 


yı — surface free energy as a sum of dispersion and polar interactions for the liquid, 
ys — surface free energy as a sum of dispersion and polar interactions for the solid, 
yi. - the dispersion part of the liquid’s surface free energy, 


y£ —the polar part of the liquid’s surface free energy, 
73 — the dispersion part of the solid’s surface free energy, 


y£ — the polar part of the solid’s surface free energy. 


The above measuring methods are: the Fowkes method [8], the Owens-Wendt method [12] and 
the Wu method [24]. The Owens-Wendt and Wu methods are based on measurements by two cali- 
brating liquids. The accuracy with which the methods identify dispersion and polar interactions 
considerably improves when a third calibrating liquid is employed. Currently the Owens-Wend 
and Wu methods are being replaced with the van Oss-Good acid-base method with three calibrat- 
ing liquids, which is based on the Lewis theory of acids and bases [11] and Fowkes’s concept of 
components [9]. The method is still under development and requires further theoretical- 
experimental research [24]. 


59 


4. Taking interfacial phenomena described by surface free energy into account when 
matching materials for sliding pairs 


Considering the above, it is proposed to match materials for friction pairs also on the basis of a 
quantitative description of interfacial phenomena, using calculations and measurements of the sur- 
face free energy of solid bodies [16]. In order to reduce friction losses in a sliding pair which is 
being designed, it is proposed to match such materials for the pair that the surface of one of them 
has a high sum of surface free energy components originating from van der Waals interactions 
while the other material’s surface has a possibly low value of the sum (fig. 3). 

The two surfaces should have high hardness, i.e. a high value of the total surface free energy. 
Thanks to the other surface’s low sum of van der Waals interactions frictional resistance will be 
reduced due to a reduction in the work of adhesion to this surface, especially in mixed friction 
conditions. 


FIRST ELEMENT 
WITH LARGE SURFACE AREA, 
1. high sum of surface free energy components 
yi and P or y" and y*, 
2. high total surface free energy and so high hardness, 
KAP 


minimization 
of friction 
losses and 
SECOND ELEMENT wear 
WITH SMALL SURFACE AREA, 


1. low sum of surface free energy components 
y? and P or y" and °°, 
2. high total surface free energy and so high hardness 
(higher than that of first element), 
KAM 


Fig.3. Scheme of matching sliding pair elements [16] 


Moreover, the large difference between the sums of surface free energy components is due to 
the materials” insolubility in the solid state, which is pointed out by Bowle in his numerous papers. 
Thanks to the use of surface free energy and its components for matching sliding pair elements this 
process can be described quantitatively [15]. Especially in the case of the second element it is pos- 
tulated that, except for the generally low value of the surface free energy component originating 
from dispersion forces, the dispersion component be dominant relative to the other components. If 
the elements are so matched, friction losses will be reduced. In mixed friction conditions, for 
which, in the author's opinion, the proposed solutions will be most beneficial, the features of the 
second element’s surface will contribute to a reduction in friction losses thanks to their reduction 
in the boundary layer formed on the surface of this element. 


5. Calculations of surface free energy of exemplary materials for sliding pairs 
5.1. Calculations of surface free energy from Young’s modulus and other parameters 


Calculations based on Young’s modulus give a full picture of surface free energy. As men- 
tioned above, one can use equation (4) for this purpose. The calculated values of surface free en- 
ergy y for the particular materials, Young’s modulus £, hardness HV, Poisson fraction v and other 
are shown in table 1. 
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Tab. 1. Surface free energy values calculated from Young’s modulus (yy), 
hardness and Poisson fraction (yyy) and surface tension at melting point (yr) 


Material Young’s | Hardness | Poisson Melting Temperature Vy Yuy Yri 
modulus HV fraction point coefficient mN/m mN/m | mN/m 
E, Y T, dy/dT, 
GPa K mN/Km 
Cast iron 100 [20] 250** 0.25 1500 [3] -0.49 [3] 875 342 2460 
EN-GJL-200 [20] 
Nitrided cast | 200 [20] 800** 0.25 1500 [3] -0.49 [3] 1250 1093 2460 
iron EN- [20] 
GJL-200 
TiN 590 [4] | 2100 [4] 0.23 3223 [4] -0.25 [3] 3503 2067 2734 
(1600)* [20] (1574)* 


* values for coating, ** measurement values 


Similarly, surface free energy values were calculated from relation (9) for surface free energy 
yny as a function of hardness HV and Poisson fraction v. The values are presented in table 6 which 
also shows surface free energy values calculated on the basis of molten solid surface tension yrs 
using relation (2). Although the presented methods of estimating surface free energy yield differ- 
ent results, they are still useful considering that there are no other methods of estimating total sur- 
face free energy. In the author’s opinion, the method based on hardness and the Poisson fraction is 
the most adequate, mainly because of the fact that hardness is resistance to the forcing of elastic 
and plastic deformations. The method’s applicability is limited to very hard (above 200 HV) mate- 
rials. 


5.2. Measurements of surface free energy components from angle of wetting 


The measurements were carried using a D2 Kri goniometer with the DSA (Drop Shape Sof t- 
ware). Three calibrating liquids with known surface free energy polar component y” and dispersion 
component y” for the Owens-Wendt, Fowkes and Kaeble methods and with known component y“” 
and component y* for the van Oss-Good acidic-basic method were used for the measurements. 
The calibrating liquids were: water, formamide and diodomethane. The measurement results 
shown in table 3 confirm that after nitriding the dispersion component increases, which is accom- 
panied by a slight increase in the sum of the components generated by van der Waals forces. Table 
4 shows the results of surface free energy dispersion and polar component measurements for dif- 
ferent coatings deposited using the PAPVD method. 


Tab. 3. Measured surface free energy components of cast iron EN-GJL-200 before and after nitriding [16] 


Material Wetting angle © Owen-Wendt Van Oss-Good 
diodomethane | formamide | water method acidic-basic method 
Sum of components: Sum of components 
EN-GJL- 54.18 57.53 69.87 y =39.20 y=36.07 
200 Dispersion component Component LW 
129.06 y” =31.92 
Polar component Component AB y** =4.15 
y=10.14 Acidic component=0.28 
Basic component=15.51 
Sum of components Sum of components 
Nitrided 47.64 54.32 75.00 y =40.41 y=39.51 
EN-GJL- Dispersion component Component LW 
200 [14] y*34.01 yo" =35.58 
Polar component Component AB T =3.93 
y7=6.40 Acidic component=0.47 
Basic component=8.26 
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Tab. 4. Sum of surface free energy components for TiN 


Material Wetting angle © Surface free energy Calculation method 
components 
diodomethane | formamide | Water Fowkes | Kaelbe 
TiN 59.27 69.04 78.80 | Sum of components y 29.44 29.45 
[32, 33] Dispersion component yf 24.33 23.75 
Polar component y? 5.11 5.69 


6. Experimental investigation of sliding pair and surface free energy of its elements 


The influence of surface free energy and its components on tribological phenomena in a sliding 
pair was verified in a roll-block tester. Clinging to the counterface the specimen forms a distrib- 
uted contact of 100 mm? with it. In accordance with the considerations presented in section 5, the 
counterspecimen was made from cast iron EN-GJL-200. Ten counterspecimens, of which five 
were subjected to vacuum nitriding, were prepared [19]. The specimens were made from cast iron 
designated in accordance with PN-EN 1560 as EN-GJL-350. As specified in section 5, their work- 
ing surfaces were coated with titanium nitride, using the plasma-assisted physical vapour deposi- 
tion (PAPVD) technique described in [21]. A linear speed of 1.25 m/s and a pressure of 5 MPa 
were used to create extremely difficult operating conditions for the sliding pairs, similar to the 
conditions prevailing in the first piston ring’s inner dead centre area at the beginning of the inter- 
nal-combustion engine’s power stroke. As a result, mixed friction was obtained, as evidenced by 
the fact that the coefficient of friction was much above 0.01 — a value typical for fluid friction [15]. 
The range of tests carried out in the TO5 microprocessor roller-block tribotester included meas- 
urements of friction forces and the wear and mass temperature of the specimen, which as time 
characteristics were archived in a computer system. After statistical handling the time characteris- 
tics were used to calculate the average friction coefficients for the sliding pairs (table 5). TiN 
coated specimen/cast iron counterspecimen pairs sliding in the presence of a lubricant — ELF 
SYNTHESE 5W50 synthetic oil (pairs: 12, 13, 15, 18). TiN coated specimen/nitrided cast iron 
counterspecimen pairs sliding in the presence of a lubricant - ELF SYNTHESE 5W5S0 synthetic 
oil ELF SYNTHESE 5W50 (pairs: 11, 14, 16, 17, 19). 


Tab. 5. Plan of tests in tribotester TOS 


Specimen number | Counterspecimen Load | Speed | Type of oil 
MPa m/s 
11 nitrided cast iron 5 1.25 synthetic 
12 cast iron synthetic 
13 cast iron synthetic 
14 nitrided cast iron synthetic 
15 cast iron synthetic 
16 nitrided cast iron synthetic 
17 nitrided cast iron synthetic 
18 cast iron synthetic 
19 nitrided cast iron synthetic 


The friction distance for one test run was set at 1350 m. The total friction distance was set at 
20 runs. Once the sliding pair cooled down after each run, wear, being the sum of the wear of the 
specimen and that of the counterspecimen, was measured. The tests were conducted according to a 
static, determinate and complete test plan. 
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Tab. 6. Friction forces and friction and wear coefficients of sliding pairs tested for distance of 27000 m 


Average Average friction Average speci- 
SLIDING PAIR tangent force coefficient men temperature 
N °C 


TiN coated specimen 
Cast iron counterspecimen 28.1 0.044 54.5 
Synthetic oil Elf Synthese 5W50 
TiN coated specimen 

Nitrided cast iron counterspecimen 8.0 0.015 40.5 
Synthetic oil Elf Synthese 5W50 


Friction coefficients as a function of friction distance were calculated from friction force meas- 
urement results, as a quotient of the friction force and the pressure multiplied by the contact sur- 
face area. The graphs of the friction coefficients are shown in figs 5 and 6. Through regression 
analysis regression functions (in the polynomial form) were fitted to all the parameter values in the 
figures showing friction coefficient graphs. Through comparative studies the coefficients of fric- 
tion of cast iron and nitrided cast iron counterspecimens sliding against TiN in the presence of 
synthetic oil were determined (tab. 6). 
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Fig.5. Coefficient of friction versus friction distance for cast iron counterface 
sliding against TiN coated specimen lubricated with synthetic oil 
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Fig. 6. Coefficient of friction versus friction distance for nitrided cast iron counterface 
sliding against TiN coated specimen lubricated with synthetic oil 
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The differences in surface free energy values between cast iron EN GJL 200 and nitrided cast 
iron EN GJL 200 resulted in considerable differences in the friction and wear coefficient values of 
the tested sliding pairs. Figures 9 show the average friction coefficients for the tested sliding pairs. 
The average surface free energy of the cast iron (y=1225 mN/m) and the nitrided cast iron (y=1600 
mN/m) from which the tested counterspecimens were made, calculated from the data shown in 
table 4, are marked on the Y-axis. Nitriding of the cast iron resulted in an increase in total surface 
free energy (table 1) and consequently, in an increase in hardness. Also the dispersion component 
value (tables 2 and 3) of this energy increased as indicated by calculations using wetting angle 
measurements. High values of the surface free energy dispersion component are conducive to the 
formation of a uniform film of oil, which is well bounded with the base thanks to good wetting. 

The causes of the lower friction coefficients are sought in the greatly reduced intensity of adhe- 
sion wear during mixed friction when microcontacts between surface irregularities occur as in the 
case of coatings with a lower sum of the dispersion component and the polar component of surface 
free energy, including the TiN coatings on the tested specimens (table 4). Also significant is the 
low value of the surface free energy polar component which for the TiN coating is several times 
lower than for the other coatings with a comparably low sum of the surface free energy dispersion 
and polar components (table 4). This is associated with a reduction in Keesom interactions at ele- 
vated temperatures [16]. 
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y=1225; yd=29,06 y=1600; yd=34,01 
surface free energy of counterface; dispersion component determined from wetting angle 
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Fig. 9. Friction coefficient of sliding pair versus surface free energy of counterface 


An additional factor which played a part in preventing adhesion wear was the use of synthetic 
oil as the lubricant. The low friction coefficients were obtained owing to the properties of the syn- 
thetic oil (mainly its ability to form an oil film already at the startup) combined with the property 
of the specimen’s surface: low intensity of adhesion wear and the property of the counterspecimen: 
good durability of the lubricating film. 


7. Conclusion 


e Wear resistance (often identified with hardness) and wettability by a lubricant (usually lubri- 
cating oil) come to the fore among the well known material properties which should character- 
ize sliding pair elements. 

e Surface free energy estimates based on Young’s modulus, hardness and surface tension at the 
fusion temperature allow one to assess general strength properties. For this only the main ma- 
terial properties are needed. 

e By calculating surface free energy dispersion component y{ from the Hamaker constant one 


can assess the wettability of a given body and consequently, the energy of adhesion to oil. 
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Then one can draw conclusions about the behaviour of the two interacting elements, especially 
in mix friction conditions. 


By calculating dispersion component y$ or Lifshitz-van der Waals component 75” and polar 
component y? or acidic-basic component 72” of surface free energy on the basis of the angle 


of wetting by calibrating liquids one can gain insight into the components resulting from inter- 
facial interactions. Using these methods one cannot calculate the total surface free energy of 
solid bodies characterized by metallic and covalent bonds. Nevertheless, measurements by 
calibrating liquids are important for the identification of the components playing a role in inter- 


facial interactions, i.e. dispersion component yi or Lifshitz-van der Waals component y” 


. . . AB 
and polar component y? or acidic-basic component 75 . 


The estimations and calculations of surface free energy and its components done for the com- 
bustion engine's PRC unit made it possible to well match the materials for the sliding pair 
elements. The goodness of the match was confirmed by the results of tribological and combus- 
tion engine tests. 
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Abstract 


The paper shows; the results of studies on the influence of an external normal harmonic force upon reduction of a 
friction force, in the system of two bodies in a planar contact. It has been pointed out that, the main reason for the 
friction force reduction is due to dynamical effect; kind of stick — slip monition. Nature of the phenomenon has been 
described and explained in the article. The computational studies have been validated with experimental results 
available in literature. Moreover a nondimensional relation, between forces acting on a slider and average velocity of 
the slider, has been found. 


Keywords: stick-slip, friction, reduction, contact, nonlinear vibration 
1. Introduction 


Interaction between the friction force and the normal forces takes place in many dynamical 
systems, e.g. in. vibrating plate compactors or vibrating feeders. The phenomenon was described 
in many publications [1-16], in context of friction reduction due to normal vibrations or forces. 
The reduction of the friction force, under normal vibrations and a variable normal force, is usually 
explained in the literature thorough: 

— decreased value of friction coefficient [15]; 

— decreased average value of contact deflection [8, 16]; 
decreased “efective” value of the reaction force [16]; 
decreased average value of the real area of contact [3, 5, 7, 14]; and 

— stick-slip motion caused by a variable normal reaction [4, 9, 10, 12, 13]. 
Although the matter has been known for a long time, there is not one widely accepted opinion for 
the reason of friction force reduction. The main goal of this paper is to prove that, the main reason 
of the friction force reduction due to an external normal harmonic force, is the stick-slip motion. 

In the considered work; the friction coefficient has a constant value, contact deflection has 
been neglected, and the mean value of the reaction equals the earthpull (equilibrium of momentum 
and impulse). Moreover according to Bowden and Tabor [1]; the value of the true contact area is 
directly proportional to the value of the normal reaction, and of the friction force is directly 
proportional to value of true contact area. Therefore the factors are not considered as a reason of 
the friction force reduction. 
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2. Model 


The aforementioned system of the two bodies in a planar contact, shown in Fig. 1, is 
investigated in the section. The system consists of a slider and a slideway. The rough surfaces 
undergoing friction create an elastoplastic interface. Nevertheless in this model, the contact 
flexibility is neglected, because the force of inertia is very slight, during vibrations excited by low 
frequency external force [9]. Thus the reaction force is the sum of earthpull and external force, 


R=Q - P, sin(2 qf?) , (1) 


where: 
R- denotes reaction force, 
Q- earthpull, 
P,- amplitude external force, 
f- frequency, 
t- time. 
As evident from Eq. (1) the average value of reaction equals earthpull (R = Q). 
A constant driving force F4 is applied to the slider (Fig.1), and tends to introduce a sliding 
motion. While the friction force Fy acting in the contact region, tends to stop the motion. The 
friction force is expressed by the Coulomb’s formula 


If x=0 then Fy= -F4 else F;=-uR sgnx , 2 


where: 
4- is a constant friction coefficient, 
x- displacement, 
Fp- friction force, 
Fa- driving force. 
The directly proportional relation between the reaction force and the friction force, during slip (Eq, 
(2)) means, that the decrease of the coeficient of friction is not considered at all. Moreover while 
smooth sliding takes place, the average friction force is equal to the product friction coefficient 
and the earthpull (F, = UR = uQ). 

The sliding motion of the slider of mass m, which is assumed to be a rigid block, is described 
by the ordinary differential equation: 


£ = m '(Fa-F), (3) 


where: 
m- denotes mass of the slider. 
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Fig. 1. Model of a simple dynamical system with external force P and Coulomb 's friction 
3. Investigation of friction force and sliding motion 


In the case, when amplitude of external force equals zero P,=0; then driving force must be 
higher or equal to friction force, to cause or maintain sliding motion in the considered system; 


FEF w= LQ = umg, (4) 


where: 
Fp- denotes the nominal static force in the absence of contact vibrations, and 
g- is gravitation acceleration. 

The sliding motion of the slider may occur at lower values of the driving force then shown 
previously (Eq. (4)), when the normal external force is present. The conditions for motion to occur 
may be written down as follows 


Fa> u(O-P„) = Frmin > (5) 


where: 

Ffmin - denotes the minimum of the friction force. 

The formula delimits areas of motion and stillness (Fig. 3.). The equation of slider motion was 
studied in detail. Firstly the time history was simulated. Diagrams (Fig. 2.) show clearly the time 
histories of the motion and the mechanism of reduction friction force. The slider begins a short- 
time slip, at the point (4), in which the potential friction force wR begins to be lower than the 
tangential force F4,. The reduction of the friction force is coupled with unloading of the contact by 
the external force. The slider accelerates to the time instance (B), in which the friction force Fy 
begins to be higher than that of the tangential force F4. At he point (B) the velocity of the slider is 
the highest. Next the friction force Fy is higher than that tangential force F4, in turn the velocity 
drops to zero at point (C), because the slider’s kinetic energy has been dissipated. Finally is the 
period of stick (immobility of the slider). The period last to the time instance, when the potential 
friction force wR is again lower than that of the tangential force F4. Therefore the stick-slip motion 
is periodical and interrupted. 

In short; temporary and periodical drops of the friction force are used for slider displacement. 
The slider moves when conditions are the best and “waits” when conditions are not good enough 
to move. The periodical motion can be seen as continuos one, when frequency of the motion is 
high. The feature may leads sometimes to misinterpretation of the phenomena. 
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0 0,01 b $ 0,02 


displacement 


0,0E+0 l 
0 0,01 4 s 0,02 


Fig. 2. Results of computer simulations illustrating time histories of forces and 
their influence on kinetic values X, X,x characterising motion of the slider 


(adopted data: P,=5N, Fg=0,7N, f=50Hz, u=0,1, g 10m/s”, m=lkg) 


In general, the occurrence of sliding motion of the slider, and its behaviour depend on: the 
value of the tangential force F4 (which is lower than Fw=umg); amplitude P, of the external force 
(which is lower than Q=mg) and the frequency f of the external force. Velocity of the slider 
depends on time histories of the tangential force F4 and the friction force F; in other words on 
impulse of the forces (lined areas in Fig. 2a). Thus by integration of the areas velocity is 
calculated. Velocity of the slider is increasing, when the lined area is rising. An analysis of factors, 
which influence on impulse of the forces, and consequently on velocity of the slider, leads to 
conclusion that; a general nondimensional solution can be received (Fig. 3). Nondimensional 
velocity of the slider U may be written down as follows 


U=, (6) 


where: 
V [m/s] - denotes average velocity of the slider. 
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Influence of P,/Q and F/Fyy , on the nondimensional velocity U of the sliding motion, has been 
presented on Fig. 3. In the figure, the area of variations of controllable data (F/Fyy and P,/Q) is 
noticeable, in which the sliding motion does not occur (U=0). This takes place for, appropriately 
small amplitudes of external force and values of the tangential force. The area, in which the sliding 
motion occurs, is noticeable too. Nondimensional velocity U increases with the rise of the 
amplitude of the external force P, and the tangential force Fy. These forces can be used to steer the 
slider. 


A R A 
OZNA | 
ZDAN 
1 


SSS? SP PS SS 


Fig. 3. Influence of P/Q and F4Fp on the nondimensional velocity U of the slider 


The model of the reduction friction force is verified. The results of simulation were validated 
with the experimental results (Fig. 4) presented by Lomakin [12]. The differences between the 
results are slight, which gives the opportunity to suppose, that the main reason of the friction force 
reduction has been presented in the article. 


71 


Experimental results 
0,8 4 ° P ,/Q=0,30 
^ P ,/Q=0,44 
Results of simulation 
——— P./Q0=0,30 


P ,/Q=0,44 


0 — p" = v "p" KE ZE) Tr —_T TEŻ p" 


0 0,4 0,8 x', mm/s 1,2 


Fig. 4. Comparison of the results of simulations against the experimental results presented by Lomakin (1955) 
(adopted data: u=0,156, m=lkg, g=10m/s’ , f=100H2) 


4. Conclusions 


The adopted model as well as the simulations and experimental results indicate, the mechanism 
of the reduction of the friction force, for planar contact of solids under the normal external force. It 
has been shown: the influence of a normal external force on the reduction of the friction force, and 
that reduction of: 

— mean value of the friction coefficient, 

— mean value of the normal reaction, 

— mean value of the real contact area, 
are not necessary in order to reduce the friction force in a sliding system. The reduction can be 
explained, on the ground of nonlinear dynamics, as a result of stick-slip motion. Moreover the 
nondimensional relation, between the forces acting on the slider and the average velocity of the 
slider, has been found as well. 
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Abstract 


In the paper the results of the thermodynamic and economic analysis of the supercritical coal fired power plant 
integrated with the carbon dioxide capture installation was shown. The paper presents the algorithm for determining 
the power of power plant and its efficiency losses due to the membrane separation of CO, from the flue gases and CO» 
compression. For the purpose of separating CO; a membrane technology was applied. Calculations concerning the 
membrane separation of CO, were carried out with the program Aspen. For the assessment of the separation process 
two indices were applied: the mole fraction of CO; in the permeate and the recovery ratio of CO,. The decision 
variables in the calculation were the pressure on the feed side and on the permeate side. The pressure on the permeate 
side is generated by a vacuum pump and on the feed side by a compressor. The power rating of require components 
determines the energy consumption of the separation and compression processes. The way of determining the 
minimum losses of the power rating and efficiency of the power plant in membrane separation process and 
compression CO; using calculated indices were shown. The power rating losses and efficiency of the power plant were 
determined for both processes. The economic analysis was calculated for power unit, taking into account investments 
and costs connected to the CO; capture installation. For the conducted analysis it is essential to determine the unit 
sale price of electricity as well as the cost of avoided emission CO). 


Keywords: supercritical power plant, capture CO2, membrane separation 
1. Introduction 


One of the priority operations in the UE is environmental protection. Poland as a members 
country has the task of reducing of the anthropogenic emissions of greenhouse gases, mainly CO2 
connected with the energy production from fossil fuels. The reduction of CO2 emission prevents 
global climate change. Also, Improved energy efficiency, effective burning of fossil fuels and the 
use of alternative energy sources will help to reduce this emission. A technique that could make it 
possible to faster achieve large reduction of greenhouse gases emission is CO»capture from power 
generation processes. 

There are many possibilities of reducing the emission of CO» in energy processes [7]. The 
main methods are pre-combustion capture, post-combustion capture and oxy-combustion. 

In the presented a paper use of a post-combustion process, including the membrane separation 
of carbon dioxide from combustion gases is described [8]. 
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2. Thermodynamic analysis and the results of calculations 


High efficiency of the electricity production has a connection to these technologies, which that 
influence more effective and efficient use of fuels. A modern power engineering rely on 
supercritical coal fired plant. Various structures and parameters with electric power of 
Ne.rer = 600 MW are presented in the paper [2]. These are reference systems for the Polish power 
engineering after the year 2010. The model of the supercritical coal power plants integrated with 
CCS installation was applied for the analysis. The diagram of such an installation is presented in 
fig. 1. 

Characteristic parameters of this block are [6]: 

— Temperature and pressure of live steam — 600°C/28.5 MPa, 

— Temperature of resuperheated steam — 620°C, 

— Pressure in the steam condenser — 5 kPa, 

—  Isentropic efficiency of all turbine stage — 0.9, 

— Mechanical efficiency and efficiency of the generator — 0.99 i 0.988, 

— Efficiency of the boiler — 0.95, 

— Efficiency of the electricity production — 0.4878, 

— Lower Heating Value — 24500 kJ/kg. 

Such a systems emits 2.086 Mg/h of flue gases of the mole fraction of CO; amounts to 
14.05%, which gases CO, emission equal to 121.3 kg/s [6]. Assuming that the block annual work 
for 8000h, the CO, emission would amount to 3 493 440 Mg. 


fuel (coal) 


Fig.1 Reference systems of supercritical carbon blocks (B- boiler, ST — steam turbine, G — generator, CND — 
condenser, DA — deaerator, CL — cooling steam, RH — reheater, h, i, | — stage of turbine (high, intermediate, low) 


For these reference systems of supercritical carbon blocks is characterizing by high efficiency 
of electricity generation, the separation of CO from flue gases was applied in order to lower its 
emission. Diagram of the installation of membrane separation and compression of CO from flue 
gases are presented in fig.2. This system consists of membrane module (M), flue gases compressor 
(C), vacuum pump (VP), CO2 compressor (C1) and heat exchangers HE1 and HE2. The aim of 
membrane module is to separate CO» flow from the flue gases (so-called permeate). The permeate 
is this part of the solution flue gases which penetrates the membrane. The driving force of the 
process is the difference of partial pressures on both sides of the membrane. 

There are two essential quantities used for the assessment of the process of separating CO2 
from flue gases: the mole fraction of CO, in the permeate (Yco2)p (purity of permeate) and carbon 
dioxide recovery ratio R, determining how much of CO, in the flue gases leaving the boiler is 
contained in the permeate flow; second one is expressed by the equation: 


76 


R= ; (1) 
Mia Ba ). 
where: 


n-flux of the flue gas, kmol/s, R-carbon dioxide recovery ratio, X,Y-mole fraction. 


Fig.2 Diagram of the system of membrane separation and compression CO; from flue gases emitted by the power 
plant presented in Fig.1 (HE — heat exchanger, C — compressor,C1 — CO, compressor, VP — vacuum pump, 
B — boiler, M — membrane module). 


Calculations concerning the membrane separation of CO» were carried out using Aspen 
software. The model of this process was developed basing on the following assumptions. The 
analysis was carried out in the membrane module consisting of cross-flow capillaries. The gases 
are considered as semi-ideals. Furthermore, constant temperature of the process at 40°C was 
assumed (hence, it is necessary to use the heat exchanger HE1 which is shown in fig.2). It 
permitted to use the constant coefficients of permeability for the gases. It was also assumed that 
the flow of 900kmol/h was passing through the membrane module with a surface area of 18000 
m”. The decision variables in the calculation were the pressure of the flue gases and of the 
permeate [4]. 

For the analysis a hybrid membrane which permits the following properties i.e. the CO2 
permeability equals to 20 m*(STP)/(m?-h-bar) and the ideal coefficient of selectivity equals to 200 
was used. The pressure of the flue gases was changed in the range 1 bar > pa > 1.35 bar and the 
permeate pressure was changed in the range 0.005 bar > ps. > 0.08 bar. The pressure on the 
permeate side was generated by a vacuum pump and on the flue gases side by a compressor. The 
power rating of require components determines the energy consumption of the membrane 
separation process. 

The power station internal load rate, including the power demand devices connected to the CO2 
separation from flue gases of reference to the power of the power plant, is mark by the symbol 6, 


and it is expressed by the equation: 


n 
6, =—_—__" r (Mc,) |1+ 
i Mel REF Nem f -LHV i ed Nic 


where: 

LHV - Lower Heating Value, kJ/kg, (Mcy) - specific molar heat capacity, kJ/kmol-K, p- pressure, bar, 
T - temperature, K, z- number of stages in the compressor or vacuum pump, K- isentropic process 
exponent, q - efficiency, %, 6 - power station internal load rate. 

Index: la, 2a - characteristic points in the system of membrane separation (fig.2), C - compressor, 
VP - vacuum pump, f - fuel, el - electrical, em - electromechanical, REF - reference system 


The equation (2) permits to determine what part of the power rating of the power plant is 


consumed during the membrane separation of CO» from the flue gases in relation to the 
assessment indices of this process e.g. carbon dioxide recovery ratio R and purity of permeate 
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(Yco2)p. The first part in the curly bracket is connected to the flue gases compression process, the 
second is connected to the vacuum pump. 

The results of the calculations from the Aspen software are presented in figure 2. The figure 
presents the line of constant values of the mole fraction of CO» in the permeate and the line of 
constant values of the recovery ratio in the plane p»a-psa (fig.3). This plot includes the line of 
constant values of the rate 6, which is calculated from the equation (2). For the calculation all of 
these isolines the algorithm of artificial networks was applied. 

This figure allows to sort and to select such conditions of the process, for which value reaches 6, 


minimum and in the same time the mole fraction of CO, in the permeate and the carbon dioxide 
recovery ratio are very high. 

The area in which the mole fraction of CO; in the permeate is greater than or equal to 0.8 and in 
the same time the recovery ratio is greater than or equal to 0.9 is marked in the figure. 

It can be found in the literature [3,5], that for the sake of the costs of compression, transport 
and deposition of carbon dioxide its permeate purity [(Y co2)pliimited Should to be contained within 
the range from 0.8 to 0.95, and the carbon dioxide recovery ratio (marked as Riimitea) ought to be 
larger than 0.8 or even 0.9. 


1,35 

P>, bar 
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Fig.3 The line of constant values of the mole fraction CO» in the permeate (Yco2)p = const, the recovery ratio CO» 
R = const, and the power station internal load rate 6; = const in the plane pr - pp for a =200 with area mark 
(Yco2)p2 0.8 iR>0.9 


In this area the minimum power station internal load rate 6, is searched for the range of value 
R > Riimited and (Yco2)p = [CY co2)Pliimited. 

If assume that (Yco2)p= 0.8 and R > 0.9 and use the characteristic from fig.3 we obtain values: 
(Yco2)p= 0.848, R= 0.9, pr= lbar and pp=0.028bar. For such parameters the process power 
station internal load rate amounts to 7.51%. 

It was assumed that separated carbon dioxide was subjected to atmospheric pressure and it had 
a high temperature. In order to transport separated CO», first it has to be cooled and next 
compressed to high pressure. Scheme of these processes is presented in fig.2. Power station 
internal load rate 6, already includes the power needed for carbon dioxide compression process. 

According to the literature [4,10,12] liquefaction of CO; for transport needs the compression 
of CO; to pressure around 100 — 120 bar. It depends on the permeate purity. In this analysis carbon 
dioxide was compressed to 100 bar. 

The electric power needed for the driving motor of the CO» compression equals to: 
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el 
Tya (Mc, je X co, R 14 Pra 1b> (3) 
Mem (Yeo, ). ic 


elC1 


where: N - power, MW 
The efficiency of electricity production in the power plant with the electric power of the motor 
compressors driving and vacuum pump taken into consideration cab be written as 


N REF = Nic t N avp + Naci) 


Merces = m, -LHV = Nearer (| — Ó)> (4) 


where: 
Nac + N ayp + N acı A 
N REF 
and it make up power station internal load rate CCS installation. 
The results of thermodynamic analysis of the influence of the CCS installation for the 
supercritical coal fired power plant are presented in table 1. 


(5) 


Ó, = 


Tab. 1 The parameters of the devices connected to the CCS installation 


Characteristic P Power Plant with CO, 
Unit 
parameters capture 
Narc MW 0 
Neve MW 45,07 
Narci MW 74,07 
ô - 0,0751 
Ó, r 0,1981 
Nel, REF s 0,4878 
Ner, ccs - 0,3912 


3. Economical analysis and the results of the calculations 


In the conducted analysis of the economical effectiveness the net present value (NPV) method 
was used. NPV is one of the fundamental and most frequently applied economic coefficients for 
the assessment of the economical effectiveness of the investments. It can be described by the 
equation: 

N 
NPV => SE : (6) 
T=l (1 +r y 

The value of the discount rate was assumed at 6,2%. In order to determine the cash flow CF, 
the total investment costs (J), profits from sales (Ser), overall costs of production (Kp), the income 
tax (Pa), changes of the working capital (Kor), amortization charges (A), interest (F) and clearance 


value of the designed installation (L) (L, L, =0 when0<t<N- 1)had to be known. 
Thus, the equation describing the cash flow CF, can be written as: 


CF, =[-J +S —(K py +P, +Ko,)+AtF +L], (7) 
The investment costs in signified analysis are given by equation: 
J =iy Na tices Nes (8) 


where: iy — unit investment cost for the power installation, €/kW, iccs — unit investment costs for the CCS 
installation, €/kW. 
The essential component of the equation (7) is the profit from sales, expressed as: 
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Sy = [Nu Cy dt. 2 (9) 
0 


where: Cei— average price of electricity, € /MWh, t.;— the annual time of operation, h. 
Total costs of production were determined as a sum of the following costs: 

K pp =K +K, +K + Keg th, +K. +4 CATS (10) 
where: Kp — cost of fuel, K, — cost of services, Kps — cost of other raw materials, Kg — other operating costs, 
Ks — environmental costs, K, — costs of maintenance and exploitation, Ax — excise duty, A — amortization 
costs, F —interest. 


The economic analysis were made for the reference system and also for similar power unit, 
taking into account the investments and costs connected to the CO2 capture installation. For the 
system with CO» separation previously calculated values permeate purity (Yco2)p= 0.848 and 
carbon dioxide recovery ratio R = 0.9 were assumed. 

Exemplary literature data [1,4,5,9,11] for different supercritical power stations taking into 
account the cost for CO2 separation were served to determine the investment costs. 

For the economic analysis the investment costs were assumed at 1100 €/kW for the reference 
system and at 1450 €/kW for the power plant with CO» capture. In the calculations, operating costs 
of capture installation were assumed at 3 €/MWh. The annual time of operation was 8000 hours. 
Amortization rate was given at 6.67 %. The constant of repairs was determined at the level of 
0.5% of the capital costs for the first ten years of operation and 1% for the following ten years. It 
has been assumed that the investment cost are spread into three years of the construction: 15% in 
the first year, 30% in the second year and 55% in the third year. It has been also assumed that the 
conduction of the power generating plant is financed in 15% by own means and in 75% by the 
commercial credit at the interest of 6%. The credit is assumed to be repaid in equal payments in 
the course of ten years. The income tax rate was assumed at 19%. Excise duty was assumed at 
5 €/MWh and the average cost of fuel was 55 €/Mg. The residual clearance value of the designed 
system is assumed at 20%-J , and the working capital as equal to zero. 


In the economic analysis the limit sale price of electricity C*% was calculated. For the 


reference systems without CO; capture it was described by the index REF and for the power plant 
with carbon capture installation it was described by the index CCS. 
For both systems the limit sale price of electricity is a quantity determined by the condition 
NPV(C*)=0, (11) 


An important economic rate for the power system with CO» capture is also the CO2 emission 
avoided (Eav) and its cost. The cost of CO» emission avoidance (Cay) is calculated according the 
equation: 

gr gr 
Ce = ‘a Crer , (12) 
AV 
where: 
Egy = Epep — Eccs> (13) 
E —CO; emission, index: REF — reference system, CCS — plant with CO, capture. 

The results of the economic analysis are presented in table 2. The analysis of susceptibility 
were also performed. The influence of the investment costs and the fuel cost were tested. The 
results of the susceptibility analysis are presented in figures (fig.4 i fig.5). 


Tab. 2 Results of the economic analysis calculation 


Unit Value 
Annual operating time h 8000 
Unit sale price of electricity C A F €/MWh 37.86 


80 


Emission CO) (Eger) 


Unit sale price of electricity after CO, separation C 5 


CO; emission after separation (Eccs) 
CO; emission avoided (EAy) 
Costs of CO» avoiding emission Cay 
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Fig.4 The influence of the investment costs of the CCS installation on the limit sale price of the electricity and on the 
cost of CO, avoided emission 
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Fig.5 The influence of the fuel costs of the limit sale price of electricity and on the cost of CO, avoided emission for 
the investment costs at 1450 €/kW 


The susceptibility analysis ware carried out. The analysis showed that the costs of CO2 
avoiding emissions decrease with decreasing the investment cost. The cost of fuel has an important 
influence on the price of electricity. 

In the susceptibility analysis costs of fuel ware tested. The price of fuel was changed by 20% 
from considered value. A change of fuel price has an influence the change of the limit sale price of 
electricity by around 4 € for and above 6 € on the costs of CO» avoided emissions. 


4. Conclusions 


The paper presents the relations between the energy consumption of the membrane separation 
of CO; from the flue gases and the recovery ratio and mole fraction of CO; in the permeate as well 
as the pressure of the flue gases and the pressure of the permeate. The power rating losses and the 
efficiency of the power plant were determined for both processes — CO, separation and 
compression. In order to decreases power and efficiency of electricity we integrated CCS 
installation with supercritical coal power plant. 

The power rating of installations connected to the CO; separation equals to 7.51% of the power 
station power. As a result the efficiency of the power plant decreases from 48.78% to 45.02%. 
Take into consideration the influence of CO» compression the power station internal load rate is 
equal to 19.81 % and the efficiency of electricity production decreases to 39.12 %. 

In order to decrease heat losses the heat coming from the flue gases and the permeate 
cooling may be used in the reference system of the steam-water system steam turbine. It makes 
possible to partly eliminate the steam bleeding. Elimination of the steam bleeding causes an 
increase of the steam flow to the steam turbine and increases the power. The efficiency of the 
electricity production is growing in the reference system. 

The economic analysis was calculated for the reference system and also for power unit, taking 
into account the investments and the costs connected to the CO» capture installation. In this 
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analysis the limit sale price of electricity and costs of CO, avoided emissions was calculated. This 
investment is profitable when the limit sale price of electricity equals to 56.91 €/MWh or purchase 
price of CO permission exceed value of 29.18 €/MgCO>. 

In the susceptibility analysis of the investment costs, price of fuel and annual operating time 
was tested. The cost of electricity generation decreases with longer annual operating time and also 
when the investment costs decrease. When the CCS installation investment costs increase to 
650 €/kW the price of electricity increase about 4,5 — 5 €. However, the costs of CO» avoided 
emissions increase about 7 —8€. The fuel price has an influence of about 4€ on the price of 
electricity and of above 6 € on costs of CO» avoided emissions. 
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Abstract 


As far as marine generating sets driven by diesel engines are concerned, it is assumed that the optimal active 
power load ranges between 70 — 90% [9]. At the same time there is no unequivocal way of assuming the value of the 
auxiliary engine excess power factor in relation to the generator’s rated active power regarded as the rated power of 
the set [3, 9, 10]. According to the outcome of the authors’ research carried out on contemporary transport vessels, 
the factor ranges within 1,05 — 1,62. In operational conditions, even at low values of excess power factor, the 
contribution of generating set working time at load 70 — 90% appears relatively short, which has been discussed in 
[10]. Eventually, the process of deterioration of auxiliary engines technical condition due to their long lasting 
operation at low loads gets accelerated and the operation turns out economically unprofitable due to the increase in 
specific fuel consumption. The paper deals with broader spectrum of individual generating sets’ active power load on 
contemporary transport vessels based on long — standing identification tests of marine electric power system loads. 


Keywords: marine generating set, auxiliary engine, active power load 


1. Introduction 


Electric power system load identification tests have been carried out on various types of 
contemporary transport vessels, e.g. container and semi - container ships, bulk carriers and general 
cargo vessels. Mainly, the value of active power peak loads at fixed time intervals, allowing for 
setting the boundary conditions of the systems’ operation, has been focused on. The tests’ 
methodology has been discussed in [11, 12]. The obtained results enable the assessment of the 
generated power of the marine electric power system and individual generating sets in operational 
conditions. According to the literature [2, 9], considering the optimal operational conditions of 
auxiliary engines, the assumed active power load of generating sets driven by diesel engines 
should range between 70 — 90% of their rated power. The tests proves that in practice such load 
value of generating sets occur relatively short when taking into account their total operational time. 
Therefore, in order to adjust the number of power generating sets in the marine power station and 
their excess power factor, it appears relevant to estimate the peak load distributions of 
contemporary operated power generation sets. Due to the paper space limit, mainly the character 
of generating sets’ (GS) active power load of contemporary container vessels has been focused 
upon. They frequently turn out to be equipped with thrusters and capable of carrying vast numbers 
of reefer containers, where the demanded power may be compared to the power of a single 
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generating set [13]. The tests outcome concerning other types of transport ships have been broadly 
discussed in the following works [2, 5, 6, 7]. 

The subject of the tests have been six container ships of various capacity: 7500 TEU, 7500 
TEU, 5500 TEU, 3050 TEU, 2200 TEU and 1100 TEU equipped with certain fixed number of 
sockets to connect reefer containers. To ease the analysis the container ships have been marked 
with Roman numbers from I to VI. The marine power stations of I, II, IV and V consist of four 
generating sets, III contains three generating sets and a shaft generator and VI is equipped with 
two generating sets and a shaft generator. The vessels I, II, III, IV and V are supplied with 
thrusters, whereas VI is equipped with emergency electric drive. 


2. Active power load of marine generating sets 


For the analysis of the active power load of the marine generating sets the data on peak loads of 
the sets at the consecutive hours of the operation have been used. In case of IV, a 3050 TEU 
container vessel, the data on 24 hour peak loads have been made use of. In fig. 1 box-and-whisker 
plots of the obtained empirical distributions of the generating set loads have been shown; 
presenting the measures of location, dispersion and asymmetry of the distributions [14, 15]. The 
symbol of the generating set in fig. 1 consists of the ship’s mark on which it was installed and its 
number it has been assigned on the vessel. 


the number of observations n=16550 


Load of GS 


0,0 


M 13 14 138 4 M3 N2 M4 v2 V4 WE ee 

O” o 

1214 02 14 W2 Ma N3 VA V3 ViN Tene 
Symbol of GS 


Fig. 1. Box-and-whisker plots of the generating sets' peak loads distributions from the tested container vessels 


The location of the empirical distribution medians (in fig. 1) clearly shows that for over 50% of 
the time of operation most of the generating sets of the container vessels, the subject of the tests, 
(except of GS III 2 and III 3), their peak load turns out lower than 70%, which means that it is 
lower than the lower limit of the optimal load range. Only for 8 out of 21 analyzed generating sets 
the values of the recorded peak loads exceeded 80%. 

Differences between peak loads distributions of particular generating sets installed on a vessel 
(if they are identical) come first of all from operational strategy accepted by the engine crew or 
ship owner technical services. Therefore, the data on peak loads in case of identical sets on 
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individual vessels have been aggregated for the sake of a more general approach to the analysis 
outcome. In this way there have been constructed box-and-whisker plots of marine generating sets’ 
peak loads distributions, shown in fig. 2. In case of vessels equipped with not identical power 
generating sets, they have been marked respectively with capital letters in the alphabetical order. 


Load of GS 


0 5 Median 
: ; : ę 7 ; O 25%-75% 
IA IB IIA IIB III IM V VI T Min-Max 
Symbol of GS 
Fig. 2. Box-and-whisker plots of peak loads distributions of generating sets of the tested container ships for the 
aggregated data 


The location of the medians of the obtained peak loads distributions for the generating sets I A, 
II A and III for the aggregated data (fig. 2) proves that for half of the operational time their hourly 
peak loads exceeded 60% of the set’s rated active power. The median highest value has been 
recorded for generating sets installed on 5500 TEU container ship (III). As a rule, the peak loads of 
values higher than 70% do not occur more frequently than for 25% of the generating sets’ 
operational time. An exception to that appear the data obtained for the generating sets of the 
container ship III, where the peak loads of values beyond 70% occurred for almost half of the time 
of operation. 


3. The characteristic of the obtained empirical distributions of the generating sets’ peak 
loads 


The obtained peak loads empirical distributions of the container vessel generating sets are 
characterized for their asymmetry and they differ from the normal distribution, which is confirmed 
by the results of the statistical tests made by means of STATISTICA 8.0 (there have been carried 
out the Shapiro-Wilk test if the sample number n < 2000, the Lilliefors and chi-square tests if the 
sample number 7 > 2000 as well as the normal probability plots have been computed). Some of the 
obtained distributions are known for their multi-modality. The example of the obtained results for 
the generating sets of 7500 TEU and 2200 TEU (I and V) container ships have been shown in fig. 
3 and 4. 
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number of observations n=1634 


the Shapiro-Wilk test: W: 


Ordered response values 


0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 
Load of GS 


Fig. 3. Normal probability plot and the results of the Shapiro-Wilk test for the generating sets’ peak loads of 7500 
TEU container ship (I A) 
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Fig. 4. Histogram of generating sets’ peak loads of 2200 TEU container vessel (V) and the plotted probability density 
function of normal distribution (solid line) with the tests of goodness of fit results inserted 
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The character of the obtained empirical data on peak loads of the container vessels’ generating 
sets does not allow for the estimation of the probability density functions of empirical distributions 
due to such tools like the program STATISTICA 8.0 by means of well known functions of 
relatively simple theoretical distributions like normal, gamma, log-normal or exponential. The 
differences between the assumed models and the empirical distributions were statistically 
significant. The problem of estimating the probability density functions of empirical distributions 
of container ships’ generating sets peak loads may be solved by the application of the so called 
distribution compositions (mixtures of distributions); however, it appears a dull, time consuming 
method and would go beyond the frames of the paper. 


4. Final remarks 


The presented analysis of the six container ships’ generating sets peak loads (fig. 2) shows that 
only three (I A, II A, III) out of eight generating sets installed, for half of the time show hourly 
peak loads higher than 60%. The position of the lower quartiles (fig. 2) in case of the six types of 
generating sets (I A, II A, II B, III, IV, V) occurs above 40% of the sets’ rated active power, thus, 
for 75% of the operational time the peak loads appear higher than this value. The least economic 
turns out the load of the generating sets of the 1100 TEU container vessel (VI). 

The obtained empirical load distributions of the generating sets appear characteristic for their 
asymmetry and multi-modality. They considerably differ from the theoretical normal distribution. 
The applied test of goodness of fit for models in the form of other simple theoretical distributions 
like gamma, log-normal or exponential have not turned out the right ones to describe empirical 
distributions. 

The collected statistical data on the peak loads values of the generating sets originate from 
container ships of various technical parameters; this refers also to the marine electric power 
systems and generating sets. Therefore, the hypothesis, that they belong to the same population, or 
in other words, that they can be treated as the realization of the same random sample, needs to be 
verified. Only such verification, on the basis of all the collected statistical material, enables to 
describe the general peak loads’ value distribution of the considered container vessels’ generating 
sets. For that purpose the Kruskal-Wallis test was applied because of the fact that the peak loads 
values distributions of the particular types of generating sets do not comply with the normal 
distribution [14, 15]. The test’s detailed structure has been dealt with in [8]. The hypothesis 
verification has been carried out by means of the statistical package STATISTICA 8.0. The results 
have been presented in the table of the program STATISTICA 8.0; shown in fig. 5. 


The Kruskal-Wallis test: H (7, N=16550) = 11439,69; p =0,000 
Symbol of GS N - total number of observation 


12902221 
31219749 


Fig. 5. The Kruskal-Wallis test results presented in the table of STATISTICA 8.0; N — total number of observations, 7 
— number of degrees of freedom of the asymptotic distribution x of the statistic H, H — the value of the Kruskal-Wallis 
test statistic, p — p-value 
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The computed p-value proves that there are no bases for accepting the hypothesis of the data on 
peak loads values of particular types of the container ships’ generating sets , in question, coming 
from the same general population at the assumed significance level a=0,05. 
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Abstract 


The paper presents a concept of improving the injection and spraying processes with the use of the oxygen 
dissolved in the diesel fuel and the results of experimental investigations carried out in order to verify it. 
The combustion process of the direct injection compression-ignition engine being normally performed at a high 
excess of combustion air factor for the whole engine operation range is affected by local deficiencies of oxygen inside 
the fuel sprays in a combustion chamber. This fact is one of the main reasons for forming the harmful compounds in 
exhaust gases such as NO, and PM. 

The aim of the presented concept is to improve the fuel spray atomization by the release of the oxygen previously 
dissolved in the fuel. The influence of dissolving the oxygen in the diesel fuel on the run of the combustion process and 
concentration toxic compounds in exhaust gas has been presented in the paper. 


Keywords: spraying processes, combustion process, emission, fuel solution 
1. Introduction 


The specificity of the combustion process being realized in the direct injection compression 
ignition engine (CI) consist in a fact that the liquid fuel in a form of fuel sprays is supplied to the 
engine combustion chamber right before the piston top dead centre (TDC). Thus a complete 
process of preparing the mixture for being burnt, i.e. a disintegration of the fuel spray into drops, 
their evaporation and mixing with air needs to be performed in a very short time. For such a way 
of fuelling the engine some local and very significant differences of the excess air factor values A 
are met (fig. 1). The values met in the combustion chamber range from the infinite high one in the 
area which is not covered by the fuel spray, throughout A ~ 0.8 + 1.5 at the edges of the fuel spray 
initializing the ignition to A = 0 being recorded in the spray core. The local deficiencies of oxygen 
occur despite the high value of the global excess of air the value of which changes with the engine 
load and ranges from A = 11 at the engine idle running to A ~ 1.4 + 1.3 at the engine full load, i.e. 
for the operating conditions corresponding to the outer engine operating characteristic. 

In case of the combustion process the local deficiencies of oxygen are one of the most 
important reasons for forming carbon oxides, hydrocarbons, and partly of forming the particulate 
matter, whereas the formation of nitrogen oxides is mainly connected with the kinetics of 
developing the flame which generates the value of the produced heat as the heat delivery speed 
determines the level of temperature in a combustion chamber. All of it causes that the fuel 
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spraying, in addition to the air swirl, is of the decisive importance for preparing the mixture. The 
quality of spraying is determined by two basic physical factors: the pressure existing in the nozzle 
area right before the nozzle whole and the pressure in the combustion chamber where the fuel 
spray is directed to. An increase in the injection pressure values improves the fuel spraying and it 
is a current preferred tendency of the developments of the injection systems for the CI engines. 
The changes which are done in the Common Rail (CR) system confirm the above statement. Every 
next generation of this system is characterized by the injection pressure values that are higher that 
ones of the previously generation. 


Fig. 1. The combustion process in the combustion chamber of a direct injection compression ignition engine 


The improvement of the fuel spraying can be achieved not only by increasing the injection 
pressure but also by modifications to the mechanism causing the disintegration of the fuel spray. 
The velocity of the fuel outflow from the nozzle is a single physical stimulus which causes the 
disintegration of fuel spray in the currently used mechanism of spraying. In order to achieve the 
improvement of spraying it is proposed to use in the discussed mechanism an additional physical 
stimulus resulting from the physical properties of the gas-in liquid solution [1, 2, 5, 6]. The 
amount of gas that can be dissolved in a liquid significantly depends on the pressure. A 
spontaneous release of gas at the non-equilibrium state caused by the pressure fuel is very 
characteristic for such a solution. The process of releasing the gas from a liquid is of a volumetric 
nature, i.e. the gas is being released simultaneously from the whole liquid volume. The energetic 
effects which accompany that process depend on the speed of the stimulus modifications and the 
gas which is released always presents a tendency to break the bonds of liquid molecules. Under 
such conditions the state of liquid is similar to that state of boiling. The presented properties of a 
liquid are very desirable in the injection system of the diesel engine. Thus a concept of using the 
effect accompanying the process of releasing the gas from a liquid for improving the existing 
mechanism of fuel oil spraying has been developed. 

This concept consists in adding the appropriate amount of air to the fuel, its dissolving under 
high pressure conditions (in a high pressure pump) and keeping it in a form of solution in a high- 
pressure section of the engine supply system (up to the nozzle) until the moment of injection, 
occurs as shown in fig. 2. In this case the assumed injection pressure determines the energetic 
solution level at which the equilibrium state is achieved. 
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Fig. 2. The illustration of a concept of the spraying improvement by gas dissolved in fuel [1] 


The use of the effect accompanying the process of releasing the air dissolved in the fuel oil 
occurs in the combustion chamber area when, due to a sudden pressure drop a serious disturbance 
of an equilibrium state occurs in a solution. The pressure drop from the injection pressure level to 
the pressure level existing in the surroundings of the fuel stream occurs at the moment of the 
injection beginning. A detailed description of this concept is presented in [1]. 


2. Test stand 


The tests were carried out on the engine test stand equipped with the direct injection 
compression-ignition test engine AVL 5804. This is an one-cylinder engine equipped with a four- 
valve cylinder head and two camshafts. The injector is situated in the cylinder head centrally in the 
cylinder axis. The engine was equipped with a conceptual supply system of a Common Rail type 
controlled by a SesubCR system — i.e. a system specially developed for electronic controlling the 
Common Rail unit. The test stand was equipped with a brake provided for realization of the set 
value of the engine crankshaft speed regardless of the engine load. The test stand also included the 
lubricating oil and cooling fluid temperature stabilization systems. 
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Fig. 3. The test stand 


The fuel was supplied to the engine by an accumulator supply system including a supply pump 
with an independent drive unit. The fuel compressed in the pump was pumped (forced) to an 
accumulator (a container) called a pressure accumulator from which it was delivered to the 
BOSCH 0445 110 131 injector. Its work was controlled by the SesubCR system. The relevant 
injection parameters i.e. the injection commencement angle and injection duration time were set 
by a computer using a program for the SesubCR system. In order to perform a correct analysis of 
the changes in the engine operation parameters it was necessary to complete the information on the 
cylinder pressure characteristics, the cylinder pressure characteristics just before the injector and 
also on the characteristics of the pulses controlling the injector operation as a function of the crank 
angle. For those reasons the engine was equipped with a piezoelectric sensor of the indicated 
pressure situated in the engine cylinder head just before the injector directly on the injection pipe 
connecting the injector with the accumulator. The voltage pulses generated by those sensors, after 
their amplification, and the signal from pulses controlling the work of the injector as well, were 
sent to the Indiset 620 system provided for recording the quick changing engine processes given in 
the voltage form. For the comparative nature of the performed examinations it was necessary to 
use two supply pumps, the conventional one, which compresses the diesel oil only, and the second 
one which makes the gas-in-fuel dissolving possible while it is being pumped (forced) (fig. 4). 
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During the examinations the tonnage oxygen taken out from the high pressure oxygen cylinder 
with use of the pressure reducing valve was used as a gas to be dissolved in fuel. Oxygen at the 
pressure value of 1 bar was supplied to the forcing section area of the pump through a non-return 
valve during the piston moving downwards (fig. 4a). As soon as the piston reveals the lower 
passage the gas supplying valve closes and the diesel oil supplying to the fuel forcing section starts 
(fig. 4b). As soon as the piston moving upwards closes the passage supplying the fuel the 
compression of the oxygen and diesel oil in the forcing section starts, during which a gas-in-oil 
dissolving process occurs. The liquid solution obtained in this way is pumped (forced) to the fuel 
accumulator (fig. 4c). From the fuel accumulator it is supplied to the injector, and next to the 
combustion chamber of the DI engine. Considering the specificity of examinations the relevant 
methodology of their performing had to be developed. 
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Fig. 4. The realization of the gas-in-diesel oil dissolving process 


3. Examination procedure 


The examination was of a comparative nature so it was carried out in two stages. In the first 
stage the engine was fuelled with diesel oil and the injection parameters were set as follows: 
- the simulated engine speed: 2000 rpm; 

- the fuel pressure in rail: 40 MPa; 
- the duration time of electric pulse opening the injector: 0.55; 0.65; 0.75; 0.87; 0.98 ms. 

In the second research stage the engine was fuelled with diesel oil -and-dissolved oxygen at the 
same settings of the above injection parameters and the oxygen pressure value in pump: | bar. 
During the examinations some measurements of the toxic exhaust gas compounds were performed. 
In each of the set engine operating points the cylinder indicated pressure values, fuel in the 
injection pipe before the injector, and characteristics of the intensity of the injector opening 
current values were additionally measured. Those measurements allowed to estimate the similarity 
of the parameters of the fuel injection realized in each stage of the performed examinations. The 
obtained parameters were subject to the mutual comparative analysis. 

4. Test results and discussion 


The evaluation of the discussed conception will be made on the basis of the comparative 
analysis performed for two groups of values. The first group includes the cylinder pressure value 
and the rate of the cylinder pressure rise. These quantities are closely connected with the 
combustion process kinetics. The second group includes the basic components of exhaust gas, the 
emission of which is subject to the limitation. They are: nitrogen oxides NOx, and exhaust gas 
smoke. 
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Having realized their research program the authors have gathered some very extensive 
comparative material which one cannot help fully presenting here. In this paper only some 
exemplary results are presented on a basis of which the tendencies of the influence of oxygen 
content in the fuel on the engine operation, that is observed in the whole research range, can be 
shown. These results are presented in a graphic form. 

The nature of the influence of the release of oxygen dissolved in a fuel on the cylinder pressure 
Pc characteristic during the combustion process is shown on an example of a single point of work 
of the engine. This influence is presented directly on the pressure characteristic as shown in fig. 5. 
In this figure the characteristic of fuel pressure in the high pressure rail Pr and the characteristic of 
a signal t controlling the opening are additionally plotted. At introducing the fuel solution and 
increasing the mass of the oxygen dissolved in it a clear tendency to shorten compression-ignition 
delay appears. The character of changes values, as one emphasized in fig. 6, recurred for all 
examined engine operating conditions. 
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Fig. 5. Cylinder pressure of engine supply with diesel fuel and diesel fuel with dissolved oxygen; 
n = 2000 rpm, Pw = 35 MPa, t = 1,21 ms 
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Fig. 6. Pressure rise rate as function of crank angle; 
n = 2000 rpm, Pw = 35 MPa, t= 1,21 ms 


The discussed character of changes is closely connected with the compression-ignition delay 
angle. After adding the oxygen to the fuel a compression ignition delay angle is subject to a 
significant shortening in comparison to one noted at supplying with fuel without oxygen. This 
shortening depends on the amount of the oxygen dissolved in a fuel and, what should be 
considered as obvious, on the location of the engine operation point on the general engine 
characteristic. Increasing the amount of the dissolved oxygen shortens the time of the compression 
ignition delay. 
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The decrease in the pressure rise rate which is distinctly seen in the combustion chamber 

should be also explained by a shortened compression ignition delay angle. The tendency of such 
association is also unambiguous: an increase in the amount of the dissolved oxygen results in the 
lower rate of the pressure rise. In case of delivering the oxygen under the pressure of 1 bar the 
maximum rate of the pressure rise was noticeable decreased. 
The results of the toxic exhaust gas compound emission measurements obtained for the individual 
engine operation points during the fuelling the engine with diesel oil containing the dissolved 
oxygen were compared with the emission values obtained for the engine fuelled with a 
conventional diesel oil (table 1). 

The relative changes in the concentration values obtained for the individual compounds are 
given in figure 7. The concentration values for the individual compounds in each operating point 
for the engine fuelled with the conventional diesel oil were assumed as a 100% concentration 
value for a given compound of the exhaust gas emission. The concentration values obtained for the 
engine fuelled with the diesel oil and dissolved oxygen were referred to that value. 


Table 1. The concentration values for the selected exhaust gas components (engine speed 2000 rpm) 


Inj ection Nitro gen Opacity 
Fuel time oxides [FSN] 
[ms] | [ppm] 

1 2 3 4 
_ 0.55 413 2.8 
5 0.65 500 4.1 
3 0.75 516 5.5 
A 0.87 531 1.2 
0.98 538 9.0 
+ 0.55 546 2.0 
a 5. 0.65 561 3.2 
Pi 0.75 598 4.2 
20 0.87 604 5.7 
A 0.98 559 8.2 


Comparison of the obtained results allows to evaluate the effect of the oxygen dissolved in diesel 
oil on the changes in concentrations of the exhaust gas components. Some changes in the 
concentrations of nitrogen oxides were noticed. The increase in concentration of approx. 30% 
occurred at low load of the engine. Next, as the engine load increased that concentration decreased 
and it was of 5-17 %. However, it still was higher than the concentration measured for the engine 
fuelled with diesel oil. The opacity in exhaust gas decreased by 5—25 %. 
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Fig. 7. Comparison of relative changes of the nitrogen oxides and opacity concentrations in exhaust gas 


5, Conclusion 


The observed changes in concentrations of the individual compounds in exhaust gas from the 
engine fuelled with diesel oil with dissolved oxygen can be undoubtedly referred to the changes 
occurring in the combustion process. The increase in the concentration of nitrogen oxides can 
prove that there is an increase in the temperature value in the flame front area and it also affects 
(results in) the increase in the combustion chamber temperature. This can be confirmed by the 
results of torque measurements the values of which increased by 5—20 % for the individual engine 
operating points. The obtained results show the positive changes of mutual correlation between the 
nitrogen oxides emission and the exhaust gas opacity especially in range of high engine load. The 
authors have gathered a rich experience resulting from their previous examinations carried out 
with the use of air and exhaust gas dissolved in diesel oil [3, 4]. Linking the present results with 
the previous ones they dare to say that fuelling the engine with diesel oil with dissolved air and 
exhaust gas gives better effects in reducing the concentrations of the toxic compounds in exhaust 
gas from the compression-ignition engines. Comparison of these examination results shows that 
the effect of releasing the gas dissolved in diesel oil, but not the oxygen concentration in the gas 
being dissolved, is a main factor affecting the improvement of the combustion process. 
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Abstract 


The paper presents a proposal of the application of the entropy analysis for the theoretical research of the 
processes occurring during the thermodynamic transformations taking place in the waste energy recovery systems in 
marine Diesel power plants. In view of the low exergy of waste energy carriers in these systems it becomes significant 
to make a proper selection of parameters in the individual points characterising the thermodynamic transformations 
of the working media. The method as presented herein consists an adaptation of the methods applied in the 
arrangements of shore power plants to the needs of the marine power plants. This article presents basic relations 
allowing to determine increments in entropy in the heating steam generation system, as well as the other basic 
elements forming the waste energy recovery systems in marine power plants. They constitute the basis to evaluate the 
performance of the individual processes and to disclose the places of occurrence of the major losses and to determine 
the manners to minimise same. The proposed method may contribute to the reduction of labour consumption of 
universally applied, traditional methods of searching the effective systems of waste energy recovery in marine Diesel 
power plants. 


Keywords: entropy analysis, waste energy recovery, marine power plants 
1. Introduction 


The specific feature of the marine systems of waste energy recovery in Diesel power plants is 
the relatively low exergy of the energy carriers. This results in the fact that the degree of waste 
energy recovery is small, and the generally applied arrangements are mainly restricted to the 
recovery of the main engine exhaust gas heat to generate steam in the amount that satisfies wholly 
or partly the ship’s heating needs. Only in few cases the available amount of waste energy allows 
additionally to generate mechanical power in gas or steam turbine on in both at the same time, 
which may be utilised to drive the generator or for the propeller drive [7]. However, in each case it 
is significant to properly select the parameter values in the individual circulation points, 
characterising the working media thermodynamic transformations. In order to make this choice in 
the appropriate manner it is necessary to conduct the relevant analyses whose results will enable 
the evaluation of the quality of the course of the individual processes and indicate the places of the 
major losses and the ways to minimise them. 

The basis for the calculations of the waste energy recovery systems in marine Diesel power 
plants consists the mass and energy balance equation systems [9]. While designing the discussed 
systems enthalpy analysis has become generally applied. Exergetic analysis [6, 8, 9] can be used 
for the proper evaluation of, in particular the sources of waste energy, whereas for the evaluation 
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of the performance quality of the processes taking place during system operation entropy analysis 
may be considered as the appropriate one. Such analysis is presently used not only at the design 
stage of new systems of shore power plants, but also to conduct the heat-flow diagnostics [5] or 
providing the guidelines for the modernisation directions of these power plants and evaluating the 
effects of their modernisation [2, 4]. 

The entropy analysis has not been found so far so generally useful for the designing of marine 
power plants. Nevertheless, it should be expected that as the solutions improving widely 
understood effectiveness of the marine power plants are sought, its application will turn also useful 
in this field. Thus this article consists an attempt to present the application of the basis of entropy 
analysis to the evaluation of the processes occurring in the systems and the elements forming the 
waste energy recovery systems in marine Diesel power plants. 

The inspiration to present in this article the entropy method as a supporting tool in the complex 
analysis of waste energy recovery in marine Diesel power plants has been provided inter alia by 
the studies [2, 3, 4, 5]. The general methods of entropy increment minimisation during the 
processes of heat exchange and in power plants have also been referred to in [1]. Its elements have 
also been presented earlier by the author in [6, 8]. 


2. Entropy Analysis of Steam and Mechanical Energy Generation System 


Considering a significant similarity of marine systems of waste energy recovery, where water 
and steam are working media, to the shore systems of steam turbine power plant systems, the 
general model consisting the basis for the entropy analysis of engine exhaust gas heat recovery in 
marine Diesel power plants can also be presented by use of the specification contained in [3]. The 
figure 1 shows general diagram of the analysed part of the system. 


Gal 


tess, m 


Exhaust Gas 


LA) Qoar Tour) 


Fig 1. General diagram of the system of heat recovery of marine Diesel power plant engine exhaust gases 


This part under discussion has been separated by means of cover sheet Kı, marked in the figure 
with the interrupted line. The diagram shows: EGB — exhaust gas boiler; H — heaters; CC — 
condensate cooler behind the heaters; T — waste heat turbine; C — condenser of turbine outlet 
steam; FT — feed tank (hotwell); CP — condensate pump; FP — feed pump. 
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Outside the cover sheet there are: waste heat boiler, condensate cooler and steam condenser. 
Processes taking place in these elements can be a subject of separate analysis which basis have 
been referred to in the further part of this paper. 

By use of denominations applied in the figure 1 the energy balance within the cover sheet can 
be presented by means of the equation: 


Qian + Qir = Q oun = Q our -Q, = Qu T pee =P,’ = Qu +P, (1) 


where: 


Q,,1 - heat flux supplied to the marine heating system within cover sheet, 

Qr - heat flux supplied to the turbocharger within cover sheet, 

Qu - heat flux carried out of ship’s heating system by the condensate cooler, 

Qt - heat flux carried out from the turbine by the condenser, 

Q, - heat flux exchanged by the system with the environment (heat losses to the environment, 
Q, - usable heat flux transferred to shipboard heaters, 


P® - (gross) power of generator driving turbine, 
PQ. - own needs of the processes taking place within cover sheet, 


P, - (net) power of turbocharger. 


The sum of entropy generated within cover sheet can be determined as: 


fia (Se, — Soy) + MorhS w — Sor )— Mu (Sm — Sj | Myr (Sir = Sur) + S, = Sen (2) 


where: 


Moy:M,y - Working medium fluxes at the heater inlet/outlet, respectively, 

Mgr My - working medium fluxes at the turbine inlet/outlet, respectively, 

Sẹ - Specific entropy of boiler supply water, 

Squ>Soq ~ Specific entropy of working medium at inlet/outlet of heaters and turbine, respectively, 
Skq Spr ~ Specific entropy of working medium at the outlet of heaters and turbine, 


Sip» Sjy - Specific entropy of working medium at the outlet of heating steam condensate cooler and 
condenser, respectively, 


$, - increase of entropy fluxes related with the Q, loss, 


S.en ~ the sum of entropy generated in the processes taking place within cover sheet. 


It should be noted that: 
Qim = gy (hou -hę ), (3); Qir = nor (hor -h;,). (4) 
Qu = Mon (kę hin | (5), Qor = Mor (lig = hir). (6) 


The average entropy temperatures for the heat flux flows in the individual system elements are 
defined as: 
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Ta =, (7) igs" —*, (8) 
Sou T Stw Sor T Sty 

= h,,, -h = h,, -h 

Tot = at" (9) Tout = at, (10) 
Sin 7 Sin Skr T Sit 


where: 


Hoyo; - specific enthalpy values of working media at the inlet to heaters and turbine, 
respectively, 

h,„.h,;, - specific enthalpy values of working media at the outlet of the heaters and turbine, 
respectively, 

h,„,,h, - specific enthalpy values of working media at the outlet from heating steam condensate 
cooler and turbine steam condenser, respectively, 

h,, - specific enthalpy of boiler supply water. 


By the application of the equations (2)-(10) the sum of the generated entropy fluxes in the 
processes taking place within the cover sheet can be determined as: 


aZ Qim is Q oun Qinr J Q ouer 4 S, , (11) 
: Tint T Tyr T, 


outH outT 


While eliminating by equation (1) Qoar from the equation (11), the sum of entropy fluxes in 
the discussed processes can also be expressed by the equation: 


S = Qin 4 Qa Qt TE Qian F Que Q © Qu P; 


a Tan Tdi Tyr T 


outT 


+8, (12) 


From the equation (12) it results that the total entropy flux generated in the analysed system of 
heat recovery depends not only on the value of heat fluxes but on average entropy temperatures 
within the process of heat flow in the individual system elements. This property should be 
considered during the system designing, inter alia while adopting the adequate values of working 
media parameters in the individual equipment. However it should be noted that the minimising of 
entropy flux increments results directly in the reduction of the system operating costs, nevertheless 
it is related with the increase of the investment outlays. Although it is already a clichć to state that 
the decision-making processes while designing marine power plants should be supported by a 
complex multicriteria analysis, it is still worthwhile to emphasise this necessity. 

By using the equation (12) the power of water heat turbocharger possible to achieve can be 
determined with the assumed circulation parameters: 


A m > > T a È > TE Q inl is 
P,= Qn T 1 | Qu Q, Jk 6... Ś,) Qui | T a | T us (13) 
inT outH inH 


The equation (13) allows to determine the actual achievable power output of waste heat 
turbocharger. Its first part shows maximum achievable value of this power output under 
assumption that the entire heat flux would be used for the production of mechanical energy and 
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that no heat losses to the environment occur, whereas the processes taking place are not 
accompanied by entropy increase. By use of the above models it is relatively simple to determine 
the optimum parameters of working media in the individual points of transformations. 

3. Entropic Analysis of Selected Equipment of Steam and Mechanical Energy Generation 
System 


While analysing the processes taking place in marine waste heat boilers there should be 
distinguished three functional types of heat exchange section. These are, in case of the most 
complex and extended system, supply water exhaust gas heater, separating sections (high and low- 
pressure) and steam heater. The processes of heat exchange between heating exhaust gases and the 
media receiving heat are accompanied by working media entropy increase. It originates both from 
the heat flow process itself as well as from the pressure losses due to friction. The latter however 
are much smaller in comparison with the former ones. Thus they will be omitted here. 

The entropy increase in waste heat boiler is equal to the sum of entropy increments of all 
bodies participating in the transformation: 


FLAS = AS,,, + AS, (14) 
where: 


AS,,„ - entropy increase of the heated medium (water and steam), 
AS., - exhaust gas entropy increase. 


The increase of exhaust gas entropy flux is calculated under assumption of its constant specific 
heat capacity: 


T, 


. d . dT T, 
AS, = | MC | E MC n—, (15) 
8 | T 8 rs T 8 Peg T, 


where: 


T,» T, - exhaust gas temperature in the beginning and end of the process, 
Ceg ~ average specific heat capacity of exhaust gases, 


m,, - exhaust gas mass flux. 


The increase of entropy flux of water (steam-water mixture or overheated steam) should be 
calculated as the difference of specific entropy values at the outlet and inlet from the section under 
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consideration “1” of the waste heat boiler: 


AŚ = ty: (si; Sow), (16) 


where: 


m, - water mass flux, 
Siwi» Sow; - Specific entropy value of the working medium at the outlet and inlet of the section under 


consideration. 


The study [9] presents the results of exergetic analysis of the processes occurring in the marine 
waste heat boiler. Its completion required the determination of entropy increments in its individual 
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sections. By use of the results of this analysis it can be stated that the biggest losses do take place 
in the evaporating section. These are related inter alia with the water properties, and more 
specifically with the relatively large value of evaporation specific enthalpy with the moderate 
pressure values of steam. Its reduction takes place along with the increase of stem in boiler, 
however, it leads to the reduction of the degree of the heat utilisation of the engine exhaust gases. 
Therefore it is of major importance while designing marine engine exhaust gas heat recovery 
systems to adopt the appropriate steam pressure. 

By use of denominations in figure 2, the increase of entropy flux in steam separator can be 
determined as [9]: 


Osw m 1 


Fig 2. General diagram of steam separator 


a" ; ; ; Sos 1 
AS = m, Biss + (1 -Xi iw ]+ MiwSiw 7 MowSow 7 m. że + sali = val > (17) 
where: 
MM, - steam and water mass fluxes at the separator outlet, 
Mop,» Mos - mass fluxes of feed water and steam-water mixture at the separator inlet, 


Osw 


81,81, - Specific entropy of steam and water at the separator outlet, 
Sos Sow >Soq, - Specific entropy of steam, boiling water and feed water at separator inlet, 
X, - steam dryness degree at separator outlet, 


k. - number of circulations in boiler. 


The adoption of the appropriate value of the number of circulations in boiler should be 
particularly emphasised. Adoption of too big value of this parameter results in the straight manner 
in the increase of the losses. On the other hand, too small value may lead to intensive 
sedimentation of boiler scale in boiler evaporating section. 

The entropy increments in the surface heat exchangers can be calculated in the similar mode as 
in case of waste heat boiler. The simple method is reading of the relevant entropy values from the 
available tables or charts. It is also possible to assume constant specific heat capacity of a medium 
and calculating the entropy increase as: 


3 T 
AS = me, In—, (18) 
T 
where: 
m - medium mass flux, 
c, - medium specific heat capacity of medium, 


T,, T, - medium temperature at the exchanger outlet and inlet, respectively. 
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Mixing-type heater, summing T-pipe and hotwell are the elements where, similar as in steam 
separator, the entropy increments are chiefly caused by the process of mixing of working media. 
Entropy increase in these elements can be determined on the basis of the readings taken from the 
appropriate entropy tables or charts of the individual media “i” at the inlet and outlet of the 
recovery system element under consideration. The entropy increase flux is determined by: 


AS =s, Dr, -D. thy Sy > (19) 
1 1 


where: 
m, - mass flux of i-medium at the inlet, 


So, - Specific entropy of i-medium at the inlet, 
s, - specific entropy of the mixture at the outlet, 


n - number of media fluxes at the inlet. 
The increase of entropy flux in throttle valve can be determined as: 


AS = th(s, —s,). (20) 
where: 
m - medium mass flux, 
S,,S, - medium specific entropy values at the outlet and inlet. 


In the effect of friction in the steam turbine the expansion process is not isentropic. In the result 
the actual expansion work is less than isentropic expansion work and the process itself is a source 
of irrecoverable exergy losses. 

The entropy value can be taken from the diagram of i-s for steam or calculate them basing on 
the appropriate relations as functions of eg pressure and temperature of steam in the relevant points 
of process taking place in the turbine. 

The liquid specific entropy increase in pump can be expressed by formula [9]: 


as=vo =p) >-1} (21) 


i 


where: 


v - liquid specific volume, 
Pi»Po - pressure at delivery and suction side, 


T - averaged liquid temperature, 
n; - pump internal efficiency. 


Conclusions 


Although the enthalpy method allows to determine the effectiveness of energy transformation 
in the designed installation, it does not allow to determine clearly the places where energy losses 
occur or size of these losses. 

The restriction of energy losses is of particular importance while designing the waste energy 
recovery systems in view of small exergy of this energy sources. 

The entropy method allows to determine the energy losses in the individual processes and 
locations of the performance of the thermodynamic transformations which in consequence allows 
to take actions aiming at their reduction. 
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The entropy method can consist the supplement to enthalpy method which is generally applied 
in designing marine power plants, including the designing of waste energy recovery in marine 
Diesel power plants. 

While designing the systems of recovery of waste energy in marine Diesel power plants it 
should be taken into consideration that the reduction of losses is accompanied in general by the 
increase of investment outlays related with the execution of energy-saving arrangement. Thus it is 
important to conduct a complex, multicriteria evaluation of the solutions possible to implement, 
including not only thermodynamic and technical aspects but also inter alia ecological and widely 
understood economic aspects. 
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Abstract 


In paper device for measuring of torque and rotational speed at the ship main drive propeller shaft transferred 
will be described.. This device connected to the power plant monitoring systemt. It may be used for continuous control 
of operating costs of the ship and for the choosing of safe and best engine operational parameters. It is possible to 
observe the actual point of engine's operational state marked against the mechanical characteristic chart. The limited, 
safe area of acceptable working points is displayed on the diagram as well. The presented torque meter has possibility 
to measure instantaneous fluctuations of the torque and engine rotational speed as a function of the shaft rotation 
angle. This is a piece of important diagnostic information for fast preliminary assessment of the load on a particular 
cylinder and of the quality of engine work. 


Keywords: ship’s propulsion unit, main engine, torque measurement, torque meter 
1. Introduction 


Safe and efficient operation of ship propulsion system requires accurate data on current load of 
the main propulsion engine. A device providing those data is a torque meter. Its basic function is 
continuous assessment of load of the main propulsion engine, and signaling when permissible 
parameters have been exceeded. The torque meter can be used also to economisation of ship’s 
motion bases on selecting an optimal relation between the cost of consumed fuel and the obtained 
transporting result. Its measure can be the speed reached by the ship, for instance. Evaluating 
instantaneous fuel consumption requires accurate data on current load of the main propulsion 
engine. 

Torque meters employed in practice on watercraft are used for measuring: 

a) shear stresses — strain gauge torque meters; 

b) shaft torsion angle - string torque meters and torque meters with toothed rings. 
In described below a torque meter is applied photo optic method for measure a torsion angle of 
propeller shaft. 
2. Characteristic of torque meter 


The principle of operation of torque meters bases on photo-optical measurement of a torsion 
angle of the propeller shaft section. Two rings with machined teeth (usually 2.2+2.4 teeth for 1 cm 
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diameter of shaft) are fixed on the shaft at a distance of abt. 400 mm. They are designed in such a 

way that their teeth are in the same plane, and are covered by one measuring head. 

The average torque and power values, measured and calculated for the period of 4 to 240 sec 
(adjustment selected by operator), are displayed as percents of nominal parameters, with 0.1% 
resolution. At the same time the rotational speed is determined at the 0.1 rev/min resolution. 

Fig. | shows a schematic block diagram of the entire configuration of the mentioned torque 
meter. The configuration includes the following basic components: 

e photo optical measuring head — signals source; 

e microprocessors main module — one-chip micro-computer; 

e Programmable Logic Controller — basic signals acquisition, arithmetical calculations, 
communication with outer units, archive measured and calculated output data, tooth rings 
technical condition self test; 

e Operating Panel (LCD colour touch panel) — displaying current results of the engine's load 
measurements (rotation speed, torque, power), data archived in memory, efficiency indexes; 
used by operator for communication and giving commands; 

e LED display — displaying current results of the ,,rotational speed-torque-power” measurements 
on the navigating bridge. 

The device has a modular structure. Its basic component is a programmable PLC controller that 
executes controlling and calculating functions for measurement and calculation data, as well as 
storing of part of them. By using proper interfaces, the controller communicates with an 
, intelligent” module used for preliminary processing of measurement signals, and with other 
terminals. It is connected with the operator's console, displays, PC computer and the system that 
monitors operation of the marine power plant. 


1. Photo 2. Microprocessor 
optical head Main Module 


USB 


R5435 - Modbus 


3. Programmable 
Logic Controller 


3.1 Display 
Unit (LED) 


5.2 Display 
Unit (LED) 


Options: 
shaft generator power 


4. Operating P anel (LCD 
colour touch panel) 


depth below keel 
— S ETHERNET 


7. digital, analog, serial 
inmit.s 


GPS NMEA 0183 


ny 
Fig. 1. Schematic block diagram of the entire configuration of torque meter 


The module used for preliminary processing of measurement signals from a photo-optical head 
(laser measurement head) was designed using a fast micro-processor system that processes signals 
in a real time of an order of several milliseconds. The optional components of torque meter are: 

e PC or Laptop unit equipped with special software for observe torque and shaft speed 
fluctuations which can help in diagnosing the engine. If PLC damage the measured torque, 
revolutions and engine power are show on PC/Laptop screen; 

e analog and digital serial data inputs for acquisition additional data; 


106 


e PC unit for archiving the measured and calculated data and visualisation of actual point of 
engine’s operation state, diagrams and trends (in case of lack of Monitoring System); 
e Monitoring System — as above and communication with the owner’s technical department. 

In the presented torque meter, new technical solutions were applied that base on experience 
gained during design, production and use of previous versions of the device. So far, several dozens 
of torque meters have been produced and installed on merchant and training vessel, as well as in 
university laboratories. The novelty of the present design is the use of a programmable logic 
controller PLC along with own codes allowing easy integration of the torque meter with any 
computer system, or a special system for monitoring operation of marine power plant devices. 


3. The torque meter as a device for efficient exploitation of ship propulsion system 


The torque, power and rotational speed of propeller shaft are the parameters that carries 
important information about technical condition of the engine, and are sensitive to they changes. 
Therefore they can be considered as a diagnostic parameters, especially torque and rotational 
speed of shaft. 
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Fig. 2. The main engine efficiency diagram in power plant monitoring system on the training ship m/s "Horyzont II". 
Marked point of actual working state of main engine 


The knowledge about the torque (or power output) and rotational speed unmistakably identifies 
a point in the propulsion engine operation area. On the Fig. 2 is presented the exemplary diagram 
with the actual point of engine's operational state marked against the mechanical characteristic 
chart. This example of efficiency engine diagram is taken from the monitoring system, which 
together with torque meter are installed on the training ship m/s “Horyzont IT”. 

The area of safe and acceptable working points is displayed on this diagram. This field is 
marked darker colour. The borders of safe and acceptable area of engine working mark: power 
characteristic for nominal or exploitation fuel supply set, propeller characteristic in captive test of 
ship, control characteristics for minimal and maximal rotational speed and characteristic of 
minimal torque (in that case zero thrust curve) [5]. The nominal propeller characteristic is marked 
as a fat line. The dashed lines are the lines of solid torque. The power and torque actual values are 
expressed as percents of nominal parameters, with 0.1% resolution. 

Along with the fuel consumption basing upon the engine characteristic, stored in the device’s 
memory, it is possible to calculate theoretical fuel consumption in each measuring cycle. The 
volume of this fuel consumption is also summed up in a “Last 24hours fuel consumption” counter. 
The result showed by the counter represents sum of power produced by the propulsion engine. 
Pulses delivered from GPS satellite navigation system or from a log, allow assessing a valuable, 
from the point of view of fuel economization, factor referred to as "normalized fuel consumption 
per 1 nautical mile”. This factor immediately responses to any change in parameters of ship’s 
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motion, as well as and to changes in trim, wind direction with respect to ship’s course, propeller 
pitch (if CPP applicable), etc. 
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Fig. 3. The examples of the diagrams on the LCD operating panel 


Mentioned above exploitation parameters are displayed not only on the efficiency diagram but 
also on the LCD handy operating panel (Fig.3). It is possible through the operating touch panel to 
put to the torque meter necessary parameters of engine, propeller shaft, number of teeth on ring 
and threshold values of signalling alarms. Reading the current measurements other exploitation 
and diagnostic parameters and their history, trends from the last 24 hours it is possible also. 

The torque meter has a memory for automatic storing and reading the results of measurements 
and calculations covering several hundreds of consecutive measurement cycles. This function is of 
essential applicability on ships on which a computer system for data storing has not been installed. 

Storing and analysing data from long time intervals, together with other parameters normally 
recorded during marine power plant operation, brings information of great value. Observation of 
trends in those data allows preparing long-term reports, as well as assessing changes in technical 
condition of propulsion system components and the risk of their possible failure. Their correlation 
with fuel consumption at given changes of operating conditions, like state of the sea, wind 
direction and force, draught and trim, propeller pitch, etc., help to make a sensible decision on the 
main engine operation. The analysis of costs resulting from changes in technical condition of the 
propulsion system facilitates making important decisions on repair, cleaning of underwater hull 
parts without docking, or modifying its components. In case of malfunction or worsened technical 
condition of any component of the propulsion system, it is possible to gain information of the 
possible location of its origin. As a result, all this simplifies and shortens servicing and diagnostic 
actions during ship operation. These characteristics of the torque meter allow considering it a 
valuable diagnostic tool. 

Essentially helpful is the torque meter in propulsion systems equipped with the controllable 
pitch propeller, as this provides opportunities for evaluating the most favourable parameters of 
engine operation for a given propeller pitch. A new design of the torque meter makes it possible to 
take into account additional engine load caused by the operation of a shaft generator. In this case, 
along with operating parameters of the main engine, the engine-propeller co-operation point is 
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evaluated as well. 

The torque meter is equipped with procedures that detect and signal failures, or dirt on the 
photo-optical head and toothed rings. 

Optional codes are being prepared for a PC computer that will allow evaluating time-histories 
of shaft torque and angular velocity in a single turn (for a slow and medium rotational speed 
engines). This is a piece of important diagnostic information for fast preliminary assessment of the 
load on a particular cylinder and of the quality of engine work [1]. On the basis of these data one 
can conclude about the power transmitted by particular cylinders to the propeller shaft. Also, a 
“piston-cylinder” system can be identified which reveals visible differences from standard 
performance [2], [6]. 
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Fig. 4. Instantaneous torque fluctuation curve, reconstructed from the frequency spectrum. (the bulk carrier m/s 
“Powstaniec Listopadowy”) 
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Fig. 5. Propeller shaft rotational speed computed and filtered at RPM=118.8 rev/min, torque M; =66.8% Mnom. ( the 
bulk carrier m/s “Powstaniec Listopadowy ”) 


It should be stressed here that the number of points corresponding to one shaft revolution in 
which the torque is measured is equal to a double number of teeth in one ring, while for the 
rotational speed the number of points is four times bigger than the number of teeth. The torque and 
speed rotations fluctuation curves computed from the measured data are significantly deformed. It 
is affected by deformations caused by, among other things, an inaccuracy in manufacturing and 
installing the toothed discs on the shaft, resonances of free torsion vibrations of the shaft, 
vibrations of the shaft deflection, propeller load fluctuations and etc. 

The deformations of measurements of low-frequency can be eliminated using a method of 
spectral analysis which was presented in [3] and [4]. A high-frequency deformation was 
eliminated using a forward-backward type filter. As example, on the Fig. 4 and 5 are presented 
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reconstructed and filtered time-histories of torque and rotation speed fluctuations in consecutive 
shaft revolutions. The presented measurements was made on the bulk carrier m/s “Powstaniec 
Listopadowy” (33767DWT, L=185m, main engine SULZER 6RL66, 8160 kW, 124 rpm). 


5. Conclusion 


Concluding, it should be stressed that described the torque meter can bring measurable profits 
resulted from safe and efficient exploitation of ship propulsion system. This opinion is backed up 
by the following arguments: 
= continuous torque control helps operate the propulsion system without the effect of increased 
current repair cost; 
= current control of mutual relation between real and calculated fuel consumption, the latter being 
indicated by the torque meter, allows evaluating changes in internal efficiency of the engine, 
thus avoiding the effect of increased repair costs and fuel investments; 

= the knowledge of calculated fuel consumption used for covering | nautical mile make possible 
to select the most profitable parameters of the propulsion system (most significant — 
controllable pitch propeller), for instance with respect to the fuel consumption. 

More important: immediate feedback of any action taken to improve vessel operation — this is 

a potential source of fuel economisation; 
= the combined knowledge of the torque and ship’s speed allows controlling changes in technical 

condition of hull and screw surfaces. This, in turn, allows making decisions on possible 

cleaning of those surfaces; 

= the information of shaft torque and rotational speed fluctuation in a single turn makes a form of 
preliminary diagnostic of the propulsion engine operation; 

= storing results of measurements and calculations, along with other parameters kept in the 
device’s memory, allows the ship owner to evaluate, and possibly correct, the way in which the 
vessel is operated. 
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Abstract 


In this paper, a proposal of a new method for assessment of changes occurring in geometric structure of elements 
of turning bearing surfaces has been presented. The bearing curve parameter value changes serve as the measurement 
of the observed changes in the function of stress amplitude, for difjerent use times. The research was carried out 
within BS 06 on the subject , Selected problems of a product manufacture and life engineering". The results of the 
performed empirical tests reveal that description by means of the proposed method of changes occurring in the 
structure of the tested surfaces, throughout the wear process, is possible. 


Keywords: surface geometric structure, bearing curve 
1. Introduction 


The surface layer of the machine elements is a factor which largely determines their properties 
and combination of surface layers of two cooperating elements determines properties of the 
kinematic pair formed by these elements. One of the quantities characteristic for the surface layer 
is the surface of geometric structure which determines the process of wear. 

Extending knowledge on profiling the surface geometric structure (SGP), elaboration and 
application of new ways for its assessment (qualitative and quantitative) has been made possible, 
also thanks to the precise and computer aided measuring equipment. New programming 
possibilities of the surface analysis in a spatial system (3D) which , in combination with improved 
precision of measuring tools, allows for observation and measurement of SGP elements in a 
nanometric scale as well as for its description by means of numerous parameters, not only areal, 
but also volumetric, spatial, hybrid and functional ones, and with the use of characteristics in the 
form of curves. The most popular example is the Abbot-Firestone’s chart, called a bearing curve. 
The parameters of the curve served for the evaluation of the investigated surface changes. 


2. Object and method of investigations 
An oblique ball bearing was accepted for the investigations. This type of bearing was selected 
due to its structural form and the resultant kinematics of its elements, big intensity of phenomena 


accompanying transformation of the surface geometric structure which improves the observation 
conditions. 
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For the investigations, stress values o, = 200MPa, 621MPa, 887MPa, 1381MPa, 1520MPa 
were accepted. Such values demonstrate changes on the whole perimeter of the ring. 

The tests results were registered for three time values t: 

- at the beginning of the investigations (beginning of period I), for ri=0 s, 

- in the period of a fixed intensity of changes ( more or less in the middle of period II), 

1=2,1-10° s, 

- at the end of the bearing life (period III ), 13=3,9-10" s. 

Measuring points with the above coordinates were accepted on the basis of results of overall 
research on oblique turning bearings, contained in works[7]. 

In figure 1, a typical wear process has been presented. Three periods can be easily seen in 
them. 

I -— fast increase in wear intensity, 

II — fixed level or slight wear changes , 

III — repeated fast increase in wear intensity. 


Fig 1. Graphic record of a typical wear process 


According to a traditional approach it is accepted that the life is a sum of period I and II. On 
the basis of investigations [7], it was found that with the life criterion in the form of the motion 
resistance level, a part of period III can also be considered as the period determining the usability 
boundary, that is the bearing life. Therefore, in this period the third point was also accepted, in 
which changes were registered. Coordinates of measuring points: A,B,C , used for measurements 
and observations, are marked in Figure 1. 

Parameters of the bearing curve were accepted as parameters defining properties of the 
examined surfaces. The curve called Abbot —Firestone’s chart describes the material distribution in 
the profile. Since the bearing curve provides information on the profile course, in a precise form, it 
is possible to read from it the profile properties, significant for surface function [4]. Parameters 
characterizing the bearing curve of roughness profile are:[3,6]: 

e Sk — height of the core roughness, um, 
e Spk — reduced height of roughness profile elevation, um, 
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e Svk reduced depth of roughness profile hollow, um, 
e Srl — bearing share of peaks ,%, 


e Sr2 — bearing share of hollows, %, 

e A — the core area, 

e Al — area of elevations filled with the material, 
e A2 — area of hollows free from the material 


Graphic interpretation of these parameters are presented in figure 2. 
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Fig.2. Characteristics of bearing curve [5] 


Parameter Spk can serve as a measurement of effective roughness depth. The value of 
parameter Spk reflects the surface abrasion resistance — the smaller it is the bigger the resistance is. 
Parameter Svk is a measurement of the cooperating surfaces capability to maintain the fluid, 
therefore, the aim of finish machining is to obtain its possibly high value. 


3. Investigation results 


Commonly used amplitude parameters do not fully define properties of the examined surfaces, 
so for this purpose the surface characteristics in the form of bearing curves, were used to extend 
the evaluation. 

In Tables 1, 2 and 3, determined values of the bearing curve parameters are compared for three 
times, in which the effects of changes in SGP were observed, depending on the stress amplitude. 


Tab. 1. Parameters characteristic for bearing curves for t,=0 s 


Parameters 
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In the first transformation period (from T; to t2), parameters Sk, Spk decreased by about 40%, 
whereas, parameter Svk underwent minor changes — it decreased in the range from 11.6% to 
24.6%. This proves that in result of the balls rolling over the path there occurred lowering of peaks 
in the direction to the roughness core. Further operation — until t3 time, caused significant changes 
(increase) of parameters: Spk and Svk by 10 times, and Sk even by 20 times. The cause of such a 
situation is an increase in the wear process intensity reflected in the form of all kinds of damages 
to the surface ( craters, cracks, grooves, etc.). 


Tab. 2. Characteristic parameters of bearing curves for 1)=2,1-10° s 


Parameters 


Comparing parameters of the bearing curve in the function of stresses, presented in Table2 and 
Figure 3, it was found that along with the stresses increase, values of all parameters increased, as 


well. For minor changes o, (200;621), the increase was a few times higher, and then slight. The 
BES Oy 


only parameters whose value in the whole stress changeability range o changed to a small degree, 
were material shares — Srl and Sr2. A slight increase in the value of parameter Srl along with time 
proves grinding in of the observed surfaces of working bearings under the influence of turning 
elements rolling on them, whereas slight changes of Sr2 parameters met the expectations as the 
action of external loads at the micro-unevenness base should be insignificant. 
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Fig. 3. Chart of bearing curve for time t, (621 MPa) 
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In Table 3 and Figure 4, there have been presented dependencies of parameters accepted for the 
surface geometric structure on stresses, for the third period of use. 

The character of changes is similar as for time tz for the smallest stresses, a large rise of the 
examined parameters values was found. 


Tab. 3. Characteristic parameters of bearing curves for 1;=3,9- 10° s 


Parameters 
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Fig. 4. Chart of bearing curve for time t; (621 MPa) 


On the basis of the results contained in tables 1,2,3 it can be found that the paths of turning 
bearings in period II are characterized by the best features of bearing surfaces. A big difference 
between values Srl and Sr2 and a small value of Sk prove that the roughness profile is of the 
plateau surface type character. Surface of this type is featured by high bearing capacity, small 
friction and abrasion resistance, that is features highly desired for working surfaces of turning 
bearings. These observations are confirmed by the expected stabilized intensity of the surface wear 
process in this period. 
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4. Summary 


Dependencies of parameters describing SGP, that is, characteristics of the bearing curves and 
stresses in the points of turning elements contact with the turning path, on the internal elements, 
produced by the bearing external loads, facilitate the choice of the surface layer properties, for 
which it will be characterized by the expected usability features. As these characteristics depend 
on the kind of finish machining and its parameters, being familiar with these dependencies should 
make easier its choice, which is of great importance for the final effect of machining. 
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Abstract 


Rapid development of unit cargo transport caused the contemporary container vessels to be adjusted to carrying 
a considerable number of reefer containers. They need to be plugged in to the marine electric power system right after 
being loaded providing they contain any cargo. They are assumed to be unplugged for no longer than twelve hours 
without any damage to the cargo. At present container vessels are equipped with even a few hundred of sockets for 
plugging in reefer containers. In case of loading the maximum number of reefer containers, they become the most 
powerful energy receiver on board the ship. The regulations of the classification institutions [4] define the marine 
electric power to be consumed by the reefer containers. The paper has been an attempt of estimation of a real value of 
power consumed by reefer containers illustrated by a container vessel 2200 TEU adjusted to carrying 350 reefer 
containers and a comparison of the achieved results with the directives of the selected classification institutions [4]. 
Moreover, a critical approach to the shipyard offices’ methods of establishing the number of reefer containers to be 
plugged in for the needs of ship owning companies ordering new vessels has been presented. 


Keywords: marine electric power system, marine generating set, reefer containers 


1. Introduction 


Within the identification tests concerning real loads of transport vessels’ electric power 
systems in operating conditions, that have been carried out by the author of the paper for years, 
observations were led on various types of vessels including container ships. The adequate 
operational data related to the highest demand for the electric power at fixed time intervals were 
collected by an observer, the engine crew member making use of the measuring instruments, the 
standard power plant equipment. The tests’ methodology and, partially, their results have been 
broadly discussed in the following works [5, 6, 7, 8, 10, 11, 12, 18]. The paper has been focused 
on the results concerning the relation between the number of the carried reefer containers and the 
load of the marine electric power system of the container vessels. 

Six various types of container vessels owned by foreign ship owning companies were the 
subject of the tests. What seems interesting, the number of their carried reefer containers appeared 
considerably smaller than the possible maximum number assumed by the ship manufacturers. The 
analysis of documents [13, 14, 15, 16, 17] of one of the companies showed that such a situation 
has been maintained for several years. At the same time, when viewing the regulations of the 
international classification institutions [4], or the container ships energy balance [1], reefer 
containers with their maximum number assumed by the shipyard designing office appear the most 
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powerful receiver of electric power on board the ship. Thus, the aim was to estimate the 
relationship between the number of the carried reefer containers and the load of the marine electric 
power systems. 

Because of the paper space limits there have been presented results achieved on a 2200 TEU 
container vessel built in Taiwan in 2003. The observations were carried out during a typical 
operational voyage from Europe to West Africa lasting for approximately a little more than thirty 
days. The number of the carried reefer containers usually appears considerably smaller than the 
possible maximum number of reefer containers to be loaded on board a ship (350 FEU). Power 
supply is provided by the marine power station consisting of four identical generating sets of rated 
active power 1200 KW each. The peak load recorded at the consecutive hours of the voyage was 
considered the measure of the active power of the marine electric power system (MEPS). The 
automatic measuring instruments of the power plant enabled the reading of the active power 
produced by the generating sets with the accuracy of 0,01 kW. When carrying the analyses the 
achieved values got rounded to 1 kW. 


2. Results of the field tests 


Due to the record of the peak loads of the 2200 TEU container vessel’s electric power system, 
performed at the consecutive hours of the voyage, the scatter plot of peak loads of the marine 
electric power system depending upon the number of reefer containers was created by means of 
STATISTICA 8.0 with the correlation between the variables determined. The scatter plot [19, 20] 
with the correlation line (full line) and the limits of 95% of confidence interval have been 
presented in fig. 1. 
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Fig. 1. The scatter plot between load of the container ship’s electric power system (MEPS) and the number of carried 
reefer containers 


Correlation between the active power load of the marine electric power system (N mers [kW]) 
and the number of carried reefer containers (k < 350) has been presented by means of the formula: 
N vers = 893,004 + 2,5111k [kW] (1) 
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with the 95% confidence interval assumed. 

In the process of testing there occurred a controlled black-out (break in the continuity of 
electric power supply) caused by the damage of the vessel sea water system and for a few hours 
power was provided by the emergency generating set. Thus, the data on peak load during the 
emergency generating set supply was rejected. In this way modified plot has been shown in fig. 2. 
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Fig. 2. The scatter plot between the load of the marine electric power system (MEPS) and the number of carried reefer 
containers without the data on peak load of the emergency generating set 


The correlation between the active power load of the marine electric power system (N mers 

[kW]) and the number of carried reefer containers (k < 350) is then presented as a formula: 
N mers = 899,8145 + 2,4158k [kW] (2) 
at the 95% confidence interval assumed. 

The rejection of the data on the peak load of the emergency generating set had negligible 
influence upon the scatter plot modification (fig. 1 and 2) as well as the function of linear 
correlation between the MEPS load and the number of reefer containers on board the ship — there 
is only few percent difference. 

According to the data provided by the ship manufacturer [1] (classification supervision of the 
vessel during the process of her construction and operation was carried out by Bureau Veritas), it 
was calculated that the power demand for a single reefer container was 6,5257 kW. Thus, for the 
basic operational state, which is a sea voyage, the marine electric power system load (with the 
linear dependence upon the number of reefer containers) is defined by the formula: 

N mers = 767 + 6,5257k [kW] (8) 
where k < 350. 

The power demand of a single reefer container, achieved on the basis of the field data, appears 
2,7 times lower than the values assumed by the manufacturer. However, it should be kept in mind 
that during the tests the number of reefer containers on board the ship did not exceed 70, which 
makes only 20% of their maximum number. 
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3. Regulations of classification institutions concerning the reefer containers power 
consumption 


On the basis of the analyses presented in the previous chapter it seems to be worth referring to 
the regulations of selected classification institutions dealing with the ways of defining the reefer 
containers power demand. Classification Societies possess their own directives for estimation of 
the reefer containers power demand. The principles of determining the reefer containers power 
consumption have been presented in fig. 3 on the basis of Germanischer Lloyd regulations [4]. In 
case of the lack of field data it is recommended that the following principles are applied when 
balancing the power demand [4]: 

e for a single 20’ container — 8,6 kW; 

e for a single 40’ container — 12,6 kW. 

The level of the assumed coincidence factor is 0,9 [4]. 
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Fig. 3. GL directives for defining the power demand of the cargo hold ventilation system and reefer containers: a) and 
c) in case of 40’ reefer containers; b) and d) in case of 20’reefer containers, according to [4] 


The above presented GL regulations appear even more protective than the ones applied during 
the construction of the vessel, the subject of the analyses presented in the paper. Such approach of 
the classification institutions may cause a considerable overestimation of the reefer containers 
power demand and lead to long lasting operation of generating sets at low load, especially in case 
of a number of carried containers smaller than the one recommended by the ship manufacturer for 
power balancing. This may endanger the shipping company to considerable financial losses due to 
the deterioration of the technical condition of auxiliary engines fed mostly with residual fuels [2, 3, 
9,10]. 
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4. Final remarks 


During the process of 2200 TEU container vessel construction, when balancing the electric 
power demand, the shipyard designing office assumed that at the total number of 350 reefer 
containers to be loaded, the load per one container equaled 6,5257 kW, that is almost three times 
higher than the result obtained due to the calculations carried out by the author on the basis of the 
field data. This low value estimated by the author may result from the relatively small number of 
reefer containers carried at the time of tests. It should be noted that not all ship owners of container 
vessels transport vast numbers of reefer containers, what is proved by the data collected on the six 
container ships, the subject of tests. Thus, the shipyard designing departments should make use of 
such documents as [13, 14, 15, 16, 17] possessed by the ship owner technical services while 
designing the number of socket connections for reefer containers to be plugged in and estimating 
their power demand. This type of exaggerated equivocation of classification companies and ship 
manufacturers may lead to considerable overestimating concerning the assessment of electric 
power demand and adjustment of the number of power generating sets, which puts the shipping 
companies at risk of financial losses. In case of 1100 TEU container vessels, the subject of the 
tests, the ship owner technical services decided to have the residual fuel oil for auxiliary engines to 
be changed for the distilled fuel due to the low loads of generating sets. This way of long lasting 
(several years) operation has significantly increased its costs. Moreover, a part of power plant fuel 
system turned out useless (e.g. residual fuel oil heating system). At the same time the auxiliary 
engines operation at low loads caused accelerated deterioration of their technical condition, 
especially the fuel injection equipment. In order to avoid this type of mistakes, it seems necessary 
to make use of field data at the stage of designing vessels. The knowledge of real electric power 
demand, or the data on the number of carried reefer containers allows for adjusting the product — a 
newly constructed vessel — to the requirements of a recipient (a shipping company). 
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Abstract 


This paper presents action of shock wave resulting from an underwater non-contact explosion, exerted on ship 
hull plating. The impulse load was considered in the range of wave regular reflection and refraction at the boundary 
of two media: water and steel. In most cases the impulse action leads to failures or damages of elements of ship power 
plant as well as shipboard equipment, however without endangering ship's floatability. Typical kinds of failures 
which recurred on ships of various tonnage, are presented, a.o., on the example of ships sailing in Red Sea waters 
during Iraq-Iran war. 


Keywords: explosion, pressure, shock wave, destruction 
1. Introductory remarks 


Underwater non-contact explosion does not cause usually ship’s sinking but only complete 
loss of its maneouvrability due to many failures in the ship’s power plant and shipboard 
equipment. And, ship’s hull plating which takes up the first impulse of load may sustain local 
plastic deformations; the loss of ship’s maneouvrability does not constitute itself a danger in the 
case of ships not engaged in warfare as it took place e.g. during Iraq-Iran war where a dozen 
or so ships flying various flags sustained failures. A good example describing consequences 
of explosion is a general cargo ship, m/s „Józef Wybicki*”[1]. In the case of ships taking part 
in warfare ( e.g. convoys to Murmansk ) such situation became extremely dangerous and usually 
led to ship’s loss. 

On all the ships such failures due to impact load were similar. The most endangered were 
elements made of fragile materials, e.g.: cracks in shaftline casings, lugs of foundations 
of combustion engines and electric motors, tearing off electric driving motors of compressors, 
pumps and hoisting winches (fig. 1, 2, 3). 

Another type of failures are plastic deformations of screw joints, bending deformations of main 
engine crankshafts and propeller shaft segments. Radiocommunication and navigation equipment 
sustains failures of another type [1, 5]. 

This way, i.e. as a result of non-contact explosion of mines, 17 ships in total have been 
damaged during one month in Red Sea waters. 

Depending on magnitude of experienced impulse load the ships had to be subjected to various 
repairs: beginning from minute repairs performed by crew personnel itself to serious repairs 
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in shipyards, lasting many months. 


rv; | 
4 


Fig. 3. Electric driving motor torn off the hoisting winch body 
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2. Reflection and refraction of two-dimensional acoustic wave 


Value of impulse load resulting from shock wave action to ship’s hull plating decides 
on acceleration to which ship equipment elements are subjected. 

In [2] is presented the load resulting from shock wave reflection from non-deformable plane, 
both in the regular and irregular range. 

In this paper the load applied to a flat wall is considered in regular range, with taking into 
account wave refraction at its passing into the other medium. 


Fig. 4. Acoustic wave reflection from and refraction at a flat wall 


Fronts of incident wave and reflected one propagate in the medium I (water). The refracted 
wave penetrates the medium II (steel). In the zone 1 limited by the front of the waves OC and OE 
the both media, I and II, are undisturbed. In the zone 2 between the incident wave and reflected 
one, OC and OD, the medium is disturbed at the parameters of the incident wave. In the zone 3 
between the dividing boundary AB and the reflected wave front OD, the medium I is disturbed 
by the reflected wave parameters. In the zone contained between the refracted wave front OE 
and the dividing boundary AB, the medium II is disturbed by the refracted wave parameters 


(fig. 4). 


Where: 

AB — boundary between two media: I and II, 

OC — incident wave front, 

OD -reflected wave front, 

OE — refracted wave front, 

a; — incidence angle, 

a2 — reflection angle, 

as — refraction angle, 

C1, C2, C3 — Wave propagation velocities in the media: I and II, respectively, 
pı, p2 — density of the media: I and II, 

P1, P2, p3 — pressure of incident, reflected and refracted wave, respectively. 
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By making use of: the regular reflection condition, the crossing point of wave fronts on the 
boundary of media, Snellius principle, as well as the continuity conditions [4] it yields: 


from the reflection condition: a, =a, i C, = cy 


and, from the Snellius principle: 
Śl = sina, ¢, 


| | (1) 
sina, sina, sina, C, 
From the continuity condition of velocity of normal displacements and pressures on the dividing 
boundary of the media it results that: 
C, COSA, — C, COSA, = C, COSA, 


(2) 
Pı + P2 = P3 
The formulas for pressure of reflected wave and refracted one have the form: 
cz 
p6; COSA, — p,c,, 1-4 sin’ a, 
Pa + 3) 
Pı c? 
PC; Cosa, + pc, 1- cz sin? a, 
1 
P _ 2pc; COS Œ; 
Py cz a 


x sD 
PC; COSA, + P,c, w ie ay 
1 


60 


p1=50 MPa 


p2 [MP a] 


Fig. 5. The pressure p, on the reflected wave front in function of the incidence angle a, and the shock wave 
pressure p, 
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p1=50 MPa 


p3 [MPa] 


Fig.6. The pressure p; on the refracted wave front in function of the incident angle a, and the shock wave 
pressure pı 


The shock wave pressure at cross-section of a perfectly stiff wall was considered in [2]. 
Formulating the equation of mass conservation, equation of momentum, and equation of the state 
Teta, one has calculated the pressure acting on ship’s hull plating ( i.e. the pressure due to reflected 
wave ), hence also the impact load. 


The pressure p2 applied to the wall in function of the incidence angle a; and the shock wave 
pressure p, is presented in fig. 7: 


——— LT ZNAC 
65 70 75 80 88 24/77 


Fig. 7. The pressure p applied to the wall in function of the incidence angle a, and the shock wave pressure p, 
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3. Remarks and conclusions 


As results from the above given formulas, the ratio of the reflected wave pressure 
and the refracted wave pressure to the incident wave pressure depends on the acoustic wave 
resistance of the media and the incidence angle œ. In the case when the medium II is more stiff 


than the medium I, i.e. c,>c, , real values of reflection coefficient are obtained at the angle 


=. f m m 
values a, śa,. =arcsin—-. In the considered case: c, ~1500— (water), c, x 5000— (steel), 
Cz s s 


a,, = arcsin 0,3 =17°30'. Ata high value of wall stiffness the pressure p; (fig. 6) and pressure p» 


(fig. 5) differs only a little to each other in this range, hence the assumption on perfectly stiff wall 
results in loads greater than real ones. 

As a result of explosion a part energy is transferred into stiff hull structure and propagated 
inside the ship through particular structural elements which serve as a kind of waveguides. 
Difraction, refraction and interference of waves takes place. In consequence apart from damages 
of elements made of fragile materials, also failures of joints of steel elements occur. 
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Abstract 


Against the background of mentioned semantic, methodical and essentials flaws and errors in articulation and 
generating of the maintenance and operating knowledge, transferred in education, scientific and popularizing 
literature as well as in realization of the maintenance activities, a proposal concerning systemic configuration of this 
knowledge was presented in this paper. Taking into account that the semantic aspects of essence articulation in 
presented proposal make the substance of these considerations, the names of scientific disciplines and the terms of 
specific notions are given in Polish language. 


Keywords: action, exploitation, knowledge, maintenance, methodology, operation, system, science, theory 


1. Introduction 


For effective realization of activities ensuring the technical systems to be set in motion — 
including the main systems of ship driver — apart from material, energy and time resources, the 
adequate knowledge on operating and maintenance is necessary. However, the effectiveness of 
such activities depends not only on essential correctness of this knowledge but also on its practical 
usability. First of all the knowledge ought to be configured and drafted correctly and clearly, in 
essence, methodology and semantics. 

From the systematic view point the operation and maintenance knowledge is a market product 
and the demand for it depends not only on its essential correctness but also on its semantic and 
methodological “wrapping”. Systemic configuration of the maintenance and operation knowledge 
increases the communicativeness of its articulation and perception. 


2. Problem 


Within the limits of this pronouncement it is possible only to relate to some subjectively 
selected principal questions. 

Quite large group of flaws to be noticed in texts fixing the knowledge — not only on the 
operation and maintenance includes the semantic — methodological flaws which consist in 
occurring inadequacy of the terms and ideas a well as some methodological deficiencies. As an 
example may be mistaking the notions such as “thesis” and “hypothesis” or using “thesis in 
situations where it should not be used. It deals also with the following relations: dzialanie — 
funkcjonowanie, działanie — proces, zależność — związek, równanie — teoria, wielkość — wartość, 
liczba — cyfra, wyznaczanie — określanie, pojęcie — termin, masa — materia ( action — functioning, 
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action — process, relationship — connection, equation — theory, magnitude — value, number — figure, 
determining — qualifying, notion — term, mass — matter ). 

To frequent cases belongs also unperceiving the defferences between a physical quantity and its 
symbol and — consequently — mistaking “equation” as symbolic structure with the relationship 
between physical quantities, or the errors consisted in personification of the things though 
attributing them the causative power. 


3. Basic notions 


Postulated usability of systemic formulation ought to be revealed beginning from the 
remembrance of necessary basic notions, starting with the key notion of “system” important for 
further considerations. In accordance with Wintgen s [6] set theory we assume that system (S) is a 
functional whole, or the set of determined elements (E), provided to of determined relations (R), 
what may be written as: 


S=<E,R> 
Consequently, it may be assumed that generator of subject knowledge — in this case the knowledge 


on maintenance and operation — is adequate subjective “scientific discipline”, treated according to 
Greniewski [1] definition, as “relatively individual informatic system” (fig.1). 


THE ELEMENT OF REAL TRANSFORMATION 


THE ELEMENT OF REAL 
DP 


IDENTIFICATION 


DP R 

OBTAINING E 
OF INFORMATION 

A 


TRANSMITING 


SPECIAL FIELD OF SCIENCE 


DN OF KNOWLEDGE 


CREATION OF SPECIAL KNOWLEDGE 


Fig. 1. Model of special field of science as a systems 


DN =< DP, J, ZP, R; W> 


where: DN — scientific discipline, 
DP — subjective domain, J — language, ZP — laws, R — rules, W — knowledge. 


During gaining applicable knowledge numerous circumstances may occur which threaten its 
credibility, quality and finally — its practicability. The reason of such a situation consists in 
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ambiguous identifications included into given system of elements. 
From the systemic view point the language of scientific discipline creates a set of elementary 
systems TPD (fig. 2) called the Ogden — Richards triangle [6] 


J=<T,P,D:RI> 
where: J — language, T, P, D — term, notion, designation, RI — relations among: T, P, D. 


The source of many diversities in meaning at interpretation of the maintenance and operations 
knowledge is lighting the univocal character requirements of relations occuring among the TPD 
system elements, i.e. among the term, notion and designation. 


symbolizing about 


REFERENCE 


describing confirming 


REFERENT | 


| SYMBOL substitute 


relationship 


Fig. 2. TPD System [5] 


Resigning the wider descripion of mentioned flaws, it is still necessery on selected exaples, to 


explain - as: hipoteza, teza, teoria, równanie, zależność, związek korelacyjny, czy 
wspolwystepowanie ( hyphotesis, thesis, theory, equation, dependence, correlation or joint 
appearance ). 


Thus, „udowodnienie hipotezy” ( proof of a hyphothesis ) — as often may be read — is an 
methodological error; the hyphothesis is not a theorem, but only an assumption, the rightness of 
which may be doubtful and its configuration would be needful. 

A mistake is also to make use of the “theasis” in empirical sciences as it a theorem which needs 
to be confirmed, what is possible in the formal sciences only. Commonly used term “równanie 
teoretyczne” ( theoretical equation ) is also a redundancy what may be easily proved on the basis 
of systemic formulation. Both, the “theory” and “equation” are the systems. Equation, as an 
element of the theory, is used to its formal articulation. 

The “theory” is an abstract structure ( R ) created to describing the relations appearing among 
concrete and abstract objects of reality; it forms a functional whole or the system (TE) 
of the laws ( P ), rules ( Z ), theorems ( T ) and hypotheses ( H ): 


TE=<P,Z,T,H;R> 


Whereas the “equation” is an abstract structure being a notion of relations among the symbols 
which, together with the relations, create a functional while, i.e the “system”. 

The fault is also to misinterpret the “eqation” — a formal structure built of the symbols, and the 
“dependence” on “correlation connections”. “Dependence is such an eqation where concrete 
physical magnitudes are the symbols and it expresses the causative — “correlation connection” 
means their coexistence described by an eqation in which concrete physical magnitudes are the 
symbols. An example of wrong interpreting the terms of “process” and “activity” ( dziatanie ) are 
controversions around the definition of “maintenance and operation” ( eksploatacja ). “Activity is 
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a sequence of elementary changes in objective state of reality fragment caused by human 
influence, whereas the “process” is a sequence of elementary changes in objective state of reality 
fragment affected by the impact of various factors, where the human participation is not necessary. 
Thus, the “activity” ( działanie ) consists in creating the process in most possible and expected 
way. In accordance with Polish language convention the name's ending for “activities” is *- anie”, 
while for the “process” — cja”[5]. 


PROCESES ACTIONS 
KNOWLEDGE KNOWLEDGE 
ABOUT SOMETHING TO SOMETHING 


FORMAL 
SCIENCES 


EMPIRICS 
SCIENCES 


THEORY 


THESES HYPOTHESES 


EQUATIONS 


Y=AX+B | CORRELATIONS 
CONECTIONS 


Fig. 3. The structural model of science 
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REASONS 
DEPENDANSE 


In the case of systems such as objectively oriented scientific disciplines, their functions — in 
general formulation (fig. 3) — may consist in shaping and configuration of knowledge: 

- about something, i.e. on the processes, phenomena and relations, or on informatively available, 
existing concrete and abstract objects of the reality; 

- to something, i.e. to realize activities consisted in creating the new existences, not appearing in 
such a shape. 

Sciences of the first group, conventionally and according to custom, are determined as the basic 
sciences, sciences of the second group — as the applied ones. The objective sphere, not fully 
covering the essential division, is division into “formal” ( formalne ) and “empirial” sciences, 
based on methodological premises. In semantic convention, the names of basic disciplines — 
sciences on the processes, phenomena and relations, have in Polish language the ending *- logia”. 
Thus, they may be conventionally named *- logie”, as opposed to appleied disciplines, supplying 
the kowedge necessary for activities, which may be called as “techniques” ( techniki ) or 
“engineerings” ( inżynierie ). From methodological view point both mentioned groups are 
distinguished in this, that “- logies” generate knowedge on the basic of information, whereas the 
“techniques” ( engineerings ) — on the basic of knowledge generated by “logies”. Techniques ( 
engineerings ), generating the knowledge objectively directed, are functioning as the tools useable 
in practical activities. What are them the sources of interesting us “maintenance and operation” 
knowledge? 


4. The maintenance and operation knowledge 


The starting point to articulate the range of knowledge useful for keeping the systems moving 
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is an identification of this action during exiting of the system ( fig. 5 ) as a phase of resource 
circulation process (fig. 4). 


REAL 


SELECTS STORE | WASTES STORE 


Y 


PUTTING STORE 


USEFUL STORE 


MATERIAL FUNCTIONAL 
STORE STORE 


RECYCLING RECOVER USEFUL 
STORE STORE 


ZASÓB ODPADOWY Z RECYKLINGU 


Fig. 4. Real circulation of the resources 


Perceiving the difference between process and activity, and using terminology in convention of the 
Polish language, the process of exhausting usable resources contained in the object will be named — 
“eksploatacja? ( exploitation ) while the activity consisted in extraction of these resources — 
“eksploataowanie” ( exploiting ). 


KONSTYTUOWANIE EKSPLOATOWANIE 
TWORZYWA OBIEKTU ZASOBU 

FORMACYJNE FUNKCJONALNEGO SUBSTANCJALNEGO 

TERIALNE | BIERNE CZYNNE 

FAZA T 
PRZEDOBIEKTOWA TRZYMANIE OBIEKTU RECYKLING 
W RUCHU FAZA 

FAZA OBIEKTOWA POOBIEKTOWA 

OKRES ISTNIENIA OBIEKTU 


Fig. 5. Phases of the objects ` life 


Fig. 5 shows that the process of object exploitation proceeds since the moment of its arising 
until the management of substantial resources contained in it ( object and postobject phases ), 
being in continuous relations with its infrastructure (fig. 6). 


OBJECT 


|| ENVIRONMENT | | 


Fig. 6. Object and environment 
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After decomposition of object infrastructure into: technological infrastructure ( everything 
connected with realization of functions attributed to the object ), exploitation ( maintenance and 
operation ) structure ensuring realization of these functions by the object, as well the 
environmented structure ( reminder of ecological and sozological environment ), the maintenance 
and operation ( exploitation ) situation of the object , in most generalized formulation, will be look 
on the fig. 7. 

During whole exploitation period the object may find itself in the state of passive (fig. 8a) or 
active exploitation. 

At the state of active exploitation, when the object is being exploited it may be in the 
conditioning (fig. 8b) or functional (fig. 7) stage. 


TECHNOLOGICAL ECOSOZOLOGICAL 
ENVIRONMENT EXPLOITATION ENVIRONMENT 
OBJECT 


EXPLOITATION ENVIRONMENT 
MAINTENENANCES SYSTEMS 


Fig.7. General configuration of object exploitation. Object in of functional exploitation state [3] 


When the object is not subjected to the effect of technological and exploitation infrastructures, 
being affected by the environmental infrastructure only, then it is in the state of passive 


exploitation (fig. 8a). 
ECOSOZOLOGICAL 
ENVIRONMENT 


ECOSOZOLOGICAL 


EXPLOITATION 
OBJECT 


EXPLOITATION |__| 


ENVIRONMENT 
OBJECT 


EXPLOITATION ENVIRONMENT 
MAINTENANCES SYSTEMS 


Fig. 8. Exploitation configurations of the object: a — passive exploitation, b — active exploitation, configuration of 
conditioning ( treatment ) exploitation 


Diversity of the phenomena processes and activities taking place during exploitation, as 
resulted from a number of various relations among the object elements and its infrastructure, 
generates demand for the knowledge necessary to realize exploitation activities, i.e. to keep the 
technical systems in motion. This knowledge is being obtained from many different scientific 
disciplines including these phenomena, processes and activities as the subject of scientific 
cognition. 

According to earlier identification of scientific disciplines, the science on exploitation as a 
process of exhausting the resources of used object, will be called “exploitology” 

( eksploatologia ). Exploitology supplies the knowledge on exploitation process, through the 
compilation of knowledge acquired in many specialistic scientific disciplines engaged in cognition 
of the phenomena affecting the course of process. 

The course of exploitation process is affected, apart from the environmental process, also by 
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exploitation activities or exploiting and therefore, getting out the useful resources contained in the 
object. The science supplying knowledge on exploiting will be named “exploitation engineering”. 
The knowledge on exploiting is formed on the basis of knowledge taken from other, kindred 
branches of engineering, as well as from the knowledge on “exploitology” and “exploitics”, 
engaged in exploitation techniques and implements. 

The whole of knowledge on exploitation, i.e. exploitation process and accompanied connected 
process, as well as on exploiting and participating in it activities, may be named “exploitation 
cognitivistics” ( kognitiwistyka eksploatacyjna ) which may be treated as some scientific 
“metadiscipline” fed by many subject — directed disciplines [3]. 

Specification of these disciplines, presented in fig. 9, may be treated as an identification of the 
state of exploitation knowledge sources, necessary to designing, generating and exploiting the 
technical objects, as well as to keeping them in motion. 


Procesy Narzedzia Działania 
LOGIE TECHNIKI INŻYNIERIE 

INŻYNIERIA EKONOMICZNA 
INŻYNIERIA 

EKOLOGIA ERGONOMICZNA 

KWALITOLOGIA CYBERNETYKA INŻYNIERIA 

MATERIAŁOLOGIA EKONOMIKA INFORMATYCZNA 

METODOLOGIA ERGONOMIKA INŻYNIERIA JAKOŚCI 

METROLOGIA INFORMATYKA INŻYNIERIA LOGISTYCZNA 

PSYCHOLOGIA MATEMATYKA INŻYNIERIA 

PRAKSEOLOGIA MECHATRONIKA MATEMATYCZNA 

REOLOGIA MECHANIKA INŻYNIERIA MECHANICZNA 

SOCJOLOGIA METRONIKA INŻYNIERIA 

SOZOLOGIA PROBABILISTYKA MECHATRONICZNA 

SYSTEMOLOGIA STATYSTYKA INŻYNIERIA 

TECHNOLOGIA TRIBONIKA METROLOGICZNA 

TRIBOLOGIA INZYNIERIA 
NIEZAWODNOŚCI 
INŻYNIERIA 
PRAKSEOLOGICZNA 
INŻYNIERIA REOLOGICZNA 
INŻYNIERIA SYSTEMÓW 
INŻYNIERIA SOZOLOGICZNA 
INŻYNIERIA 
TRIBOLOGICZNA 
INŻYNIERIA WIEDZY 

EKSPLOATOLOGIA EKSPLOATYKA INŻYNIERIA 
EKSPLOATACJI 


KOGNITYWISTYKA EKSPLOATACYJNA 


INSTRUMENTARIUM 


PROCESY 


DZIAŁANIA 


TECHNIKI 


Fig. 9. Model of exploitation knowledge resources 
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In particular it deal with the combustion engines, as well as the machines and devices necessary 
to keep these engines in motion within the technological structure such as a ship, and also within 
environmental structure where the ship occurs. 

Structural model of exploitation as a process, considering the activities and exploitation 
states, may be presented in form as given in fig. 10. 

There are presented in the model (fig. 10) the states and activities occurring in objective stage 
of exploitation, thus refered to that, what in exploiting means keeping the system in motion. The 
feedbecks in particular stages of exploitation were neglected in the figure. It was also necessary to 
resign the description of post — object exploiting stage, or the recycling phase, including 
management of usable object resource which left in substancial form after using up the resources 
for functioning. 


Fig. 10. Model of objects exploiting. 1 — classification, 2 — operators service, 3 — supplying service, 4 — awaiting, 
5 — treatment service, 6 — storage, 7 — removal ( liquidation )[4] 


5. Termination 


Presented systemic identification of wide spectrum of the scientific disciplines as the 
knowledge source, showed very strong objective differentiation and wide thematic extension of 
exploitation knowledge [2]. This multiweft specificity inclines to assumption that the exploitation 
knowledge is generated in the sphere of some “meta — science”, which may be defined as 
“exploitation cognitivistics” ( kognitywistyka eksploatacyjna ). 

It is rather difficult to expect that on the basis of knowledge from one discipline will be possible 
to realize the tasks defined in subject matter of this conference: they include designing, production 
and exploitation with the aspect of keeping in motion the ship driving systems, considering the 
safety, diagnostics, ecological abilities, energy consumption, functionality, logistics as well as the 
reliability, serviceability, optimiziation, rationality, controllability and usability. 
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Abstract 


The following article presents the metod of determining the ship's propulsion engine operation basing on the 
engines work parameters by means of dimensional analysis .The ships propulsion engine activity, as noticed by J. 
Girtler in his works [4,5], can be used for its diagnostics. of ship propulsion. Diagnostics engines increases safety of 
ships movement and at the same time protects the sea environment from pollution in case of its sinking. According to 
Girtler engine operation can be considered as, a new physical quantity of dimension Joul multiplied by second 


[J i s]. This quantity can be determined on the basis of algebraical diagram of dimensional analysis constructed by S. 


Drobot. This diagram allows to control the correctness of conclusion rules, in respect of mathematics, used in 
numerical functions of ship propulsion engine operation. 


Keywords: ship propulsion engine operation dimensional analysis, diagnostics of ship engine 
1. Introduction 


The Basic operating problem of ship propulsion systems is diagnostics of ship propulsion 
engines. Loss of operational capability of ship propulsion engine endangers the ships movement 
and in case of ship's sinking may cause dangerous pollution of sea environment. Diagnostics of 
ship propulsion engines, as noticed by J. Girtler in the works [1,2,3], can be carried out by means 
of engine operation. Engines operation is interpreted as energy transmission in form of work or 
heat to the surroundings and expressed by the product of Joul and second [1,2]. 

Identification of engine operation by means of dimensional analysis in aspect of its usability to 
diagnostics has been carried out in the present article. The engine operation has been treated 
according to J. Girtler Works [1,2,3] as a new physical quantity to be used in diagnostics of ship 
propulsion engines. 


2. Forms of dimensional functions of ship propulsion engine performance 
Engine operation as dimensional quantity together with other quantities of his type 
characterizing the movement of ship propulsion system belongs to dimensional space elements. 


Products of dimensional space elements create abelian group together with involution of real 
exponent. It allows to describe dimensional space by means of positive real numbers. These 
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numbers create the underspace of dimensional space. It means that both dimensional and 
nondimensional quantities belong to dimensional space. Out of dimensional space elements we 
can select, a determined by space dimension, amount of dimensionally independent quantities 
called the space base [4,5,6]. 

Elements of the same dimensional space can be arguments of the function defined as not an 
ordinary function and called a dimensional one. A dimensional function must equally well 
describe engines operation in each configuration of units and automatically fulfill the condition of 
invariance. Apart from this condition it must fulfill the condition of dimensional homogeneity. If 
in dimensionally invariant and homogeneous function there are independent and dependent 
arguments so an the basis of Buckingham theorem we can express the last by means of numerical 
function. Such information can be obtained only by means of experiment. 

Analysis of ship engine activity in conditions of its operating creates great difficulties as it 
requires the knowledge of defining functions and also the knowledge of dynamic features of ship 
propulsion system. 

The activity of ship propulsion engine in the course of its operation is defined by the following 
parameters: 

-the effective Power of the engine N, 

- torque of the engine M, 

- rotational speed of the engine n, 


- fuel volume consumed by the engine V 

- supercharging air compression p, 

- time of engine activity t. 

The function form of ship engine propulsion operation D can be determined on the basis of 
functional dependence among the above mentioned parameters. They have adequate dimensions in 
an accepted system of measure units. Besides on the basis of measurement we can attribute to 
them defined numerical values. Taking into account the above mentioned premises quantities we 
can write the following dimensional functions of ship propulsion engine operation D as below: 


D -o| Man, b, r) (1) 


B= of, h, it) (2) 


where: 
Q - symbol of dimensional function, 


2 


2 
M- torque of the engine in k = | ; 
s 


n — rotational speed of the engine in El A 
s 


. 3 
v- fuel volume consumed by the engine in Hi 
s 


N — effective power of the ship propulsion engine in l £ m | ; 


p — supercharging air pressure in | ——— |, 
m:s 


t- time of engine activity in [s] ; 
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D — activity of ship propulsion engine in k = | i 
s 

Functions of ship propulsion engine operation (1) and (2) are described in dimensional space 
of the third grade, which means that among function arguments there are three dimensionally 
independent argument creating so called dimensional base. All possibilities of their choice in given 
functions (1) and (2) are presented in tables 1 and 2. Forms of numerical function of ship 
propulsion engine operation can be determined on the basis of measurements carried out in the 
course of its operation. 


Tab. 1. Choice possibilities of arguments dimensionally independent, so called dimensional bases, in function of ship 
propulsion engine operation D= of Mny p „ 


Ordinal | Form of dimensional function | Dimensional base | comments 
number 
1 M 
D= fpg) 
. M,n,p 
v:p 
= SE ,=n-t 
ASA 
2 M 
D=g(6,9,)— 
n 
i M,n,v 
pnv Wn, 
?, = a7 ZE $, =n-t 
3 
v 
D= fløup) = 
M-n n, vp 
dy = Sorte $, =n-t 
Pp:yV 
4 W 2 
D= f(b,.9,)-¥ pt 
M š 
dy = e: $, =n-t v pt 
v pt 
5 M? 
D= f(6,,6,)-— 
v: p . 
: v, pM 
n:-M v- p:t 
d, yee es % = = 
v: p 
6 
D=f 6,9. "Mt 
p, M 
p:t:v 
=n t, PT A= 
d, emery, 
7 D=f(8,,6,)-M-t 
ptv Mtv 
p, =n: ty, = M 
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Tab. 2. Choice possibilities of arguments dimensionally independent so called dimensional bases in function of ship 
propulsion engine operation D= af ny r) 


Ordinal | The form of dimensional function | Dimensional base | comments 
number 
1 N 

D=f(Y,)-> : 

n Nav 

go, =n-t 
2] =f) 

$,=n't N,v,t 


3. Determination of numerical function form of ship propulsion engine activity on the basis 
of dimensional argument measurements 


To measure work parameters of propulsion engine of general cargo vessel with displacement of 
5500 DWT all devices and measurement apparatus installed as ships standard equipment have 
been used. All measurements the results used in this work were taken during normal 17 days long 
voyage. Those measurements were taken at the time of engines steady work (excluding 
maneouevres) four times, a day at 8,11,14 and 20 o'clock of ship's time. The general cargo vessel 
was propelled by 5RD68 engine made by H. Cegielski — Sulzer. The measurement results 
concerning the arguments of dimensional function of ship propulsion engine performance in the 
course of its operation in nondimensional form have been presented in drawing 1. Dependence of 
variable dependent on independent ones allows, in the domain of real numbers, to determine the 
numerical form of operational function by means of multiple regression and define its constant 
coefficients. 

Drawing | presents measurement caordinates of numerical function of ship propulsion engine 
activity with sharply outlined linear dependence. These coordinates have been fitted to multiple 
regression which helped to obtain the equation of the following form: 


D ei 
a wera + 6,27-10°n-t - 0,001 (3) 
Or 
D=107%p.vi? + 6,27n-t?-M - 107 .Mt Ga) 
where:- 


symbols like in formula (1). 


Correlation coefficient of the above fitting to the straight line amounts to: 
r = 0,99999572 => r° =0,99999144 and after correction r° = 0,99999092 


which gives a standard estimation error equal to 0,00952. 
Calculations that fit the straight line to measurement data included in table 4 were carried out by 
means of the programme STATISTICA 
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Tab. 3. Parameters measurements results of ship propulsion engine opetration which are arguments of dimensional 
function of its operation 


Number | time of | rotational | torque fuel effective | engine Super 
operation | speed kgm? volume | power operation | charging 
t-10° „|>| uj consump. | N-10 D-10” air 

à 3 
[s] $ | k -m° k -m? | Po 
s s? s kg 
Pa 

1 0 2,0833 22094 0,0209 2,8471 0 39226,8 

2 0,0648 2,0900 22511 0,0194 2,9104 1240,65 41188,14 

3 0,1512 2,0733 22094 0,0215 2,8339 6576,57 41188,14 

4 0,1620 2,0650 22511 0,0224 2,8751 1661,34 43149,48 

5 0,1728 2,0733 22511 0,0229 2,8868 8751,93 41188,14 

6 0,1944 2,0833 22094 0,0224 2,8471 10923,92 |41188,14 

7 0,2376 2,0633 22094 0,0215 2,8199 16161,78 | 41188,14 

8 0,2484 1,9900 21052 0,0192 2,5919 16241,63 | 34323,45 

9 0,2592 1,9867 21052 0,0194 2,5875 17655,32 | 35304,12 

10 0,2808 2,0567 22084 0,0231 2,8111 22502,09 | 41188,14 

11 0,3240 2,0633 22094 0,0203 2,8199 30068,14 | 39226,80 

12 0,3348 2,0633 22094 0,0203 2,8199 32106,10 | 39226,80 

13 0,3456 2,0633 22094 0,0202 2,8199 34210,87 | 38246,14 

14 0,3672 2,0683 22094 0,0203 2,8265 38714,45 | 39226,80 

15 0,4104 2,1017 22719 0,0222 2,9538 50530,61 | 41188,14 

16 0,4320 2,1200 22511 0,0208 2,9523 55960,04 | 43149,48 

17 0,4536 2,1217 23136 0,0215 3,0361 63459,73 | 41188,14 

18 0,4968 2,1067 22511 0,0213 2,9332 73542,86 | 41188,14 

19 0,5076 2,1017 22511 0,0214 2,9266 76592,92 | 41188,14 

20 0,5184 2,0967 22511 0,0212 2,9192 79696,81 | 41188,14 

21 0,5400 2,1083 22302 0,0228 2,9089 86147,68 | 42168,81 

22 0,5616 2,0667 22511 0,0244 2,8780 92194,77 | 42168,81 

23 0,5616 2,0950 22511 0,0219 2,9170 93457,23 | 41188,14 

24 0,5832 2,1200 22302 0,0238 2,9251 101040,29 | 47072,16 

25 0,5832 2,1233 22928 0,0227 3,0111 104038,10 | 46091,49 

26 0,6264 2,1283 22928 0,0227 3,0185 120304,64 | 45110,82 

27 0,6372 2,1317 22928 0,0226 3,0229 124102,79 | 45110,82 

28 0,6480 2,1300 22928 0,0222 3,0207 128847,40 | 45110,82 

29 0,6696 2,1317 23343 0,0223 3,0413 140182,41 | 45110,82 

30 0,7128 2,1283 23553 0,0221 3,1001 160027,39 | 45110,82 

31 0,7236 2,1267 23553 0,0232 3,0979 164789,47 | 47072,16 

32 0,7344 3,1733 24595 0,0247 3,3053 181138,80 | 52956,18 

33 0,7560 2,0917 22928 0,0229 2,9663 172222,16 | 46091,49 

34 0,8208 2,0967 23344 0,0219 3,0273 207188,73 | 43149,48 

35 0,8424 2,1283 22928 0,0224 3,0185 217578,31 | 44130,15 

36 0,9288 2,1300 23553 0,0222 3,1023 271925,61 | 45110,82 
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Fig. 1. Measurement coordinates of numerical function of ship propulsion engine activity in an established 
dimensional base p,t,M of the form C=f (A, B). (tablel- position 6). Explanations: ,__D_ +qe-nondimensional index 


of engine operation, A=n-t-10° - similarity invariant of engine rotational speed, p _ P :¥:t -similarity invariant of fuel 
M 
volume consumed by the engine, the rest of symbols like in formula (1) 


In dimensional function of ship propulsion engine defined by the formula (2) one can select 
four different dimensional bases , but only two are correct, in respect of dimensional structure. 
Structures of numerical functions obtained from formula (2), were given in table 2, position 2, 
similarity invariable of numerical function concerning engine operation is independent of its 
rotational speed as shown in drawing 2. On the other hand its value depends on engine work 
parameters in established conditions. 
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Fig. 2. Probability invariant measurements of numerical function concerning ship engine operation in an established 


dimensional base N,y,t of the form H = f(X) (look table 2 position 2). Explanations: H=_D_-nondimensional index 
Nee 
of ship engine operation, X =n-t - similarity invariant of engine rotational speed, 
the rest of symbols like in formula (2) 


Peak occurence taking place in similarity invariant in numerical function of engine operation 
shown in drawing, is caused by the average of propulsion engine. After twelve days of the ship's 
voyage broke the liner of engine cylinder head on third unit. First it caused, a decrease of 
operation value brought about by a faulty action of the third system and next an increase of 
operation value as a result of augmentem work of the remaining systems. Numerical structure 
given in table 2 position 2 on the basis of invariant changes of the engine operation similarity 
(numerical function dependent variable) can be recorded in the following way, drawing 2. 


jj 501015 0S 0-15 0577510" 


3 = 


5 (4) 
5 =0,1016 > 0,5775-10° <n-t 
or 
D=0,1015-N-t? >0<n-t<0,5775-10° (43) 
a 
D=0,1016- N -t° >0,5775-10' <n-t 
where: 


symbols like in formula (2). 
4. Conclusions 
Two different forms of numerical functions concerning propulsion engine operation (1) and 


(2) are possible to obtain, by means of algebraic diagram of dimensional analysis given by 
S. Drobot. Different numerical structures are given in tables 1 and 2. 
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These structures allow us to define superficially, dependence among dimensional quantities 
describing ship operation. Obtained on their basis numerical functions estimators of engine 
operation are their models can be defined accuratelly according to the established coefficients 
determined on the basis of engine work parameters. 

Formulas (3) and (4) which define the ship engine operation can be treated only as correct 
proposals is respect of dimensional aspect. Forms (1) and (2) of dimensional function of engine 
operation differ one from another due to number of explanatory variables. 

Reducing one explanatory variable [function form (2)] make it possibile to obtain numerical 
function structure of one variable. In case of five explanatory variables (form 1) one can obtain 
numerical function structures of two variables. 

The best numerical function model of ship propulsion engine is the one of a simple form and 
easy to interpret physically. Numerical function models of operation distinquish themselves by 
the fact that they take into account the essential quantities which describe the work of an engine 
depending on the time of its operation. So they are of dynamic character and due to this can be 
used for diagnostic and prognostic purpose. 

The form of numerical function models of ship propulsion engine operation can be defined on 
the basis of engine operation parameters. They can be true only for the engine at which the 
measurements were carried out. 
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Abstract 


Objectivity and rationality in making decision, assumed optimal in given conditions, forces to apply an 
evaluating (quantitative) approach to the problem, hence to search for such their parameters (indices) which, in a 
given decision situation, can be deemed most adequate. 

And, to precisely determine the task it is necessary to specify also its duration time, apart from conditions in 
which it will be realized.. When considering propulsion engine, i.e. the main element of ship propulsion system, 
especially important becomes not only the problem which amount of energy could be at one’s disposal but also within 
which time interval it could be delivered. Therefore apart from applying the commonly used reliability indices, it 
seems sensible to consider the operation in such evaluating approach as it could be determined by energy and time 
simultaneously. The presented method may be deemed a valuable supplement to the ways have been applied so far of 
description of reliability features of the driving system, considered crucial for ship power plant and ship itself. 


Keywords: operation, ship power plant, Poisson process, Markov process 


1. Introduction 


As described elsewhere [4, 5], in operational practice of ship devices ( like in the case of other 
complex functional systems of mechanical devices) working life of the same elements installed in 
different power plants is not an adequate criterion for unambiguous determination of an operation 
strategy because it is not an unambiguous measure of their wear and tear. 

In decision making situation, reality of operation can be better represented by introduction of at 
least one random variable as a model parameter though it leads solely to more or less probable 
conclusions , that results from the fact that to particular values of decision variables not only one 
value of criterion function but many values occurring with different probabilities, are attributed. 
To the most often applied probabilistic decision models belong those in which the expected value 
of consequence of taking the decision is selected to aid in choosing optimum value of decision 
variable [10]. Hence in choosing an optimum value of a decision variable one should be aided by 
the expected value of criterion function [4]. 
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In such situation it is easier - from formal point of view - to present decision-making procedure 
in one of the structural forms most commonly met, namely : decision tree or decision table. In 
general case the decision tree takes the form shown in Fig. 1. 

For the below presented tree, the criterion function is constituted by maximization of the 
expected value of the consequence c(d;, s;), which , for particular decision tree nodes which 
symbolize the fact of taking a given decision dj , can be determined as follows [8, 9]: 


E(e/d,)=>[p(s,)/4, -eld,,s,) i=1,2 k  j=1,2. n (1) 
i=l 


And, it should be observed that the decision situation is deterministic as it consists in choosing 
only one decision out of n possible ones, in spite of that occurrence probabilities of the state s; (i = 
1,2, ..., k) under assumption of taking the decision dj (i = 1,2, ..., n) appear. 


DECISION-MAKING 
SITUATION 


+nanazaazzzana 


Fig.1. General case of decision tree. n — number of considered decisions; p(s;)/d; — conditional probability of 
occurrence of the state s; in the case of taking the decision d,; c(d;, si) — consequence of occurrence of the state s; 
in the case of taking the decision d, 


Application of the decision procedure shown in Fig. 1. requires, apart from determination of a 
repertoire of possible decisions, to know how to determine values of the conditional probabilities 
p(si)/dj. 

In practice, predictions concerning tasks to be realized are usually based on the wide-understood 
notion of service reliability of an object or system. And, when considering the notion of service 
reliability of power devices ( e.g. ship diesel engines ) attention should be paid to that from the 
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user’s point of view the most important problem is quality of execution of a given task ( in 
extreme case — its not fulfillment). This way the notion of reliability is inseparably associated with 
unambiguous determination of the task in question. 

Moreover to precisely determine a task , apart from assuming conditions in which it will be 
performed, it is necessary to specify is duration time. The problem is especially important in such 
domains as sea shipping where specificity of tasks is as a rule associated with necessity of the 
functioning of crucial mechanisms and devices (e.g. installed on ships) for long periods. 

Therefore it is especially important not only which quantity of energy would be available 
during operation of a given power device but also for which duration time it could be delivered. 

The presented approach is made realizable by considering engine’s operation (further referred 
only to ship main propulsion engine) in such evaluating way as to make simultaneous determining 
it by energy and time possible. 

In this case the operation (D) within the time interval [0, t] can be interpreted as a physical 
quantity determined by the product of the time-variable energy E = f(t) and time t , which can be 
generally expressed by the following relation [7]: 


D=fE(e\dr=20[M, ned (2) 


2. Estimation of occurrence probability values of ship motion limitations resulting from the 
lowering of total efficiency of ship main propulsion engine 


One of the possible ways of estimation of occurrence probability values of limitations in ship 
motion with a given speed is the above mentioned quantitative evaluation of operation. The way is 
as much versatile that it is possible to apply it in the case when results of operational investigations 
are lacking, and adopted assumptions and model parameters are determined exclusively on the 
basis of engine’s technical and operational documentation and determined task realization 
conditions. 

In the case of ship main engine, with a view of taking into account the so-called design sea 
margin as well as service power margin [10] for the engine operating under partial loads , the 
process of the decreasing of available power output ( hence also of the possible operation Dy) will 
be realized in two phases: 

e In the first phase only an increase of hourly fuel oil consumption will take place ( at a 
relatively constant value of developed engine torque ), hence operational cost will also 
increase; 

e inthe second phase a limitation of effective power developed by the engine will appear due 
to design limitations and lack of possible increasing fuel charge. 

If partial engine load is assumed constant the phenomenon can be interpreted as follows: 

e In the first phase the time-variable drop of total engine efficiency results first of all in 
increasing its hourly fuel oil consumption (increase of specific fuel oil consumption). It can 
be described as a series of the recordable events F consisting in increasing the fuel charge 
gp” by the increment Ag, at a relatively constant value of the torque M, ( appropriate to a 
given engine load state) . This way an increase of engine’s operational cost is generated, 
however without any limitations imposed on ship motion parameters, in principle; 

e Gradual degradation processes during further engine operation result in that the recordable 
events U which consist in decreasing the engine torque Mo at the constant fuel 
consumption gy ( 1.8. gp = Gp max), to occur. Further long-lasting operation of the engine 
results in significant worsening its characteristics which impose serious limitations on ship 
motion with a given speed or course. In heavy weather conditions such situation will 
obviously lead to producing a hazard to ship safety. 
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In the context of the above mentioned quality of task realization to know the following data 
becomes important in making decisions with the use of a probabilistic decision process : 
a) expected value of increased cost of task realization, resulting from increased fuel oil 
consumption, 
b) value of occurrence probability of such number of F events which cause the fuel oil charge 
gp” to increase up to the value Gp max , and ship motion limitations to occur subsequently. 
In the considered case of ship main diesel engine the problem defined in a) can be solved by 
making use of the assumption that the number of repetitions N4gp of the event F within the time 
interval (0, t) is a random variable of non-negative integer values. 
The dependence of the random variable on time constitutes the stochastic process {N(t) t> o}. 


3. An example of application of stochastic process theory to estimation of value of the 
probabilities p(s;)/d; 


Under assumptions on stationarity [3, 11], lack of consequences and flow singularity, the 
Poisson’s homogeneous process can be applied to the process of increasing the fuel charge gp” as 
a result of decreasing the engine's total efficiency q. (in steady load conditions of the engine), and 
the random variable N 4gp is characterized by the distribution [1]: 


k 
PIN, =k)= b i) exp(-A,1); k=12,...,n (3) 
where: 
àf- a constant interpreted as occurrence intensity of the event F (increasing gp” by the 
value Agp). 


Main particulars of the engine, on which the carried out calculations are based, deal with 
ORT-flex60C-B Wartsila engine (according to [12]); hence for the selected contract parameters : 
the contract power output Nx = 80% Nari = 17370 kW, engine speed nx = 90% nn; = 102,6 rpm it 
can be assumed that the state of engine operation under 85% load is typical one as in this point of 
engine's working area specific fuel oil consumption value reaches its minimum. Hence the 
parameters of the typical state of operation are determined as follows: 


e n=97,2 rpm, 
e N.=14 764,5 kW, 
= Ss 
e g.=169,4 wn: 
Taking into account the above mentioned , one can determine value of the fuel charge gi” as 
follows: 
m 14764,5kW - aan : > 
85% e Õe S 
= fe = =473 4 
p S d 21 2 9 
60 s 
G pmax = gp” can be simultaneously assumed as precise data concerning parameters of 
injection apparatuses are lacking, and: 
19107kW - an : > 
pag WI (5) 
3 Ew, 9 
60 s 


therefore: 
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AG ymax =Bp gi” =9,8g x10g (6) 

Next, value of TGpmax at which value of Gpmax resulting from degradation processes is reached, 
should be estimated. In the case of lacking appropriate results of operational investigations the 
only one, practically useful way of its determination is to analyze engine’s technical and 
maintenance documentation on the basis of which the most unreliable engine’s elements and 
systems can be selected in line with indications of engine’s manufacturer. 

In the considered example such analysis revealed that the most unreliable functional subsystem 
of the engine is its fuel system (which complies with many results of other independent research 
[6]) . Therefore when analyzing the service time-scale of the above mentioned system as well as 
other maintenance operations, to assume T gpmax = 1000 h has been deemed justified. 

If accuracy class of the applied flowmeters is assumed equal to 0,5 and value of recorded fuel 
flow rates is also assumed, then in the considered case the value of Agp can be determined on the 
level of Ay, = 0,004 kg/s. Therefore on the basis of the above mentioned assumptions it can be 
stated that during the time TGpmax the following number of F events which would result in 
imposing limitations in realizing the assumed ship sailing speeds: 

x% 
Ni, == Cente ŚRI 95a (7) 
ý Ag 4 


p 


Hence the determined value of Ar amounts to: 


U 
Ag 


„AG. 3 1 
"Tye, 1000 h 


A (8) 


Finally, for the considered example the occurrence probability of k’ number of F events 
amounts to: 


-expl 3-107 -t) (9) 


k 
P(N exp( dt) = 


Agp 


(3-10 -t7 
! 


and, the function P(N = k')= f(t) is presented in Fig.2. 


Moreover , to make the figure more clear, the occurrence probability of a lower number (k = 2) 
as well as higher ones (k = 4 and 5) of F events is also shown in it. Interpretation of the 
probabilities presented in function of time confirms expectations as it can be observed that 
occurrence probability of a lower number k of F events decreases along with time in favour of 
their greater values. In decision-making process such analysis of the probabilities in question 
makes it possible to estimate expected consequence values and to select a task realization variant. 
The assumptions concerning T Gpmax aS well as Agp may be deemed somewhat doubtful , however 
in any instant their determined values can be replaced by different, more realistic ones, but the 
procedure itself does not undergo any modifications. 


Agp 
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Fig. 2. P(N =k)= SO 
In the next phase the above presented problem was solved by using a model in the form of the 
Markov process 4X '(t) ; t = 0}. The state transition graph of the considered process with taking 
into account the values obtained from the relation (8) is presented in Fig. (3) 


ot A12 A23 
A30 


Fig. 3. State transition graph of the process {X’(t) ; t > 0} 


and the initial distribution of the proces: 
Po =P{X (0) = xo} = 1, 


pi = P{X (0)=x;} =0 dlai=1,2, 3. (10) 


Therefore the set of Kolmogorov — Smirnov equations takes the form: 
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AP(t) _ ho P(t) + Ay BO 

dt 
dRG) _ „| : 

wie du HO) + Agi: AO 
dF,(t) _ . ; 
e e PO +n AO (1) 
AO _ 4 l 

= Axo BG) +A: P E) 
P(H)+PO+PRO+RO=1 


The set can be transformed by applying Laplace transform [2] as well as the assumed initial 
distribution of the process to the set of linear equations in the domain of transforms having the 
following form: 


s-Fy(s)-1=-hy, Py (8) + Aso PB (5) 
s+ P(s)=—Ay P (8) + ho Pi (5) 

s: P (8) =—Ays Pi (8) +A F (5) (12) 
s: P; (8) =—Agy P (0) +A P, (8) 
AO+EO+RO+RO=l 


Attempting to solving the above mentioned set one should determine the transition intensity 
Aij. To estimate values of the parameters is possible provided expected values of the random 
variables Tij are known. In practice the mean time of remaining the process in the state x; 
provided the next will be the state x; , can be considered to be the estimator E(T;). In such case the 
searched for value A is expressed as follows: 


Ay =a Z hy = (13) 


In the case in question, as 3 transitions of the process, which result from the lowering of total 
engine efficiency by the same value, are considered, it is justified to assume the following 
assumptions associated with determining Aij values : 

e inthe most common case of lacking results of operational investigations it is possible to 
assume that in a rational system of operation to distinguish duration time of any of the 
states except the state xk , i.e. in this case - x3 is unjustified (as there is no reason to 
claim that any of the states should last for a longer or shorter time), 

e itcan be therefore assumed that E(To1) ~ E(T 12) ~ E(T23), hence Ag, %12 AB =A, 

e the mean duration times of the distinguished states xo, x; and x2 should be then 
assumed the same, it is therefore justified to assume its value equal to: 

T G 


aa 7 
pa (14) 


On substitution of the data for the considered case the value : A = 0,003 was obtained. By 
taking into account the above specified assumptions and the notation : A39 = U , the set of 
equations (12) can be presented in the following form: 
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s- P} (s)-1=-A- Pi(s) +H Pi (s) 
s-P'(s)=-A- B'(s) +A- P(S) 

s- P} (s)=-A- Pi (s) +A-P'(s) (15) 
s-P;(s) =-A-P/(1) +h- PIs) 
P(t) + PA) +P) + P0) =1 


The solving of the above mentioned set of equations in the domain of transforms and the 
subsequent executing of inverse Laplace transform makes it possible to find the searched for 
distribution P3(t) graphically presented in Fig. 4. _ 
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Fig.4. The probability distribution P3(t) 
4. Summary 


The presented method may be deemed a valuable supplement to the ways have been applied so 
far of description of reliability features of the driving system, considered crucial for ship power 
plant and ship itself . Its basic advantage consists in connecting energy assessment with duration 
of time in which a task is realized. The time is very important in the case of sea shipping tasks 
usually long lasting. 

Making use of it one is able to determine, for a given instant , useful work ( useful energy ) 
which can be produced by the whole driving system, as well as to determine value of occurrence 
probability of such number of F events which would cause additional limitations to form during 
realization of a task ( due to not possible propelling the ship with a given speed) or to make its 
realization impossible at all. Value of the probability can be hence taken as that of reliability 
index and implemented in making operational decisions. Its additional advantage is versatility 
which makes it possible to apply it to reliability analysis of every ship power device or subsystem 
including those not being machines, e.g. heat exchangers. 
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Abstract 


The World is strongly dependent on crude oil for its transport needs. In order to diminish this dependence, 
we need to introduce clean, CO>-efficient, secure and affordable transportation fuels. The current production of 
liquid biofuels in the EU25 is less than 1% of the market. Recent assessments have concluded that the 2010 tar- 
gets, 18 Mtoe used in the transport sector, are unlikely to be achieved. There can be three basic possibilities of 
accomplishing this target: i) the use of alcohols (first of all ethanol) and their mixing with petrol; ii) the use of 
fatty acids esters (methyl or ethyl) of vegetable oils and their mixing with diesel fuel, iii) the use of synthetic 
hydrocarbons of the synthetic gas coming from biomass resources and eventually their mixing with other “clas- 
sical” hydrocarbons. 

This paper presents a new way of utilizing alcohols as fuels for a diesel engine. It is proposed to use heavy 
alcohols as a mix with vegetable oils and conventional diesel fuel. It is presented another way to use alcohols. 
Namely the use of heavy alcohols as a solvent for vegetable oil (called the biomix or BM) and after the obtain- 
ment of the density which would be similar to diesel fuel, mixing the biomix with diesel fuel to obtain biomixdie- 
sel (BMD). This solution will be shown for example with butanol as heavy alcohol, rape oil as vegetable oil and 
conventional diesel fuel. The investigations are carried out with a simple diesel engine on the engine test bed. 
Main parameters of engine (power output, torque, specific fuel consumption) and the main exhaust gas compo- 
nents (in this case CO, NO,, PM) were investigated. There were better results achieved than one expected. Con- 
trary to existing experiences, the maximum of power output and the torque of engine is higher in the whole range 
of the rotatory speed of the engine crankshaft when the engine biomixdiesel (BMD) is reinforced. The addition of 
the biomix component to fuel influences the specific fuel consumption. Generally with the larger part of the bio- 
mix component the specific fuel consumption grows. Because the power of engine also grows up one should 
expect that in exploitation the specific fuel consumption should not increase. It is very important that this fuel 
could be used to reinforce old, existing now and the future diesel engines. 

It’s worth paying attention that the presented solution in which a virgin vegetable oil (contrary to today’s 
situation in which as a fuel ingredient we have only fatty esters) is an ingredient for fuel. 

The production of butanol is known (from biomass and in other way with electrolysis of ethanol). The possi- 
bility to get butanol from ethanol gives a very good perspective for the use of ethanol from today’s overproduc- 
tion and moreover without the essential change of infrastructure. 

All this leads to the conclusion that fulfilling the expected requirements of European Union regarding the 
biofuels is fully possible.The introduction of new fuel needs carrying out of a lot of complicated investigations, 
but chosen direction may be interesting. 


Keywords: diesel engine, butanol, biofuel, ecofuel 


Introduction 


The World is strongly dependent on crude oil for its transport needs. In order to diminish 
this dependence, we need to introduce clean, CO»-efficient, secure and affordable transporta- 
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tion fuels. The current production of liquid biofuels in the EU25 is less than 1% of the market. 
Recent assessments have concluded that the 2010 targets, 18 Mtonnes used in the transport 
sector, are unlikely to be achieved. 

There can be three basic possibilities of accomplishing this target: i) use of alcohols (first 
of all ethanol) and their mixing with petrol; ii) the use of fatty acids esters (methyl or ethyl) of 
vegetable oils and their mixing with diesel fuel, iii) the use of synthetic hydrocarbons of the 
synthetic gas coming from biomass resources and eventually their mixing with other “classic- 
al” hydrocarbons. 

It is known that today in EU we have an overproduction of ethanol. Together we have an 
overproduction of petrol. We need more and more fuels for diesel engines. Also in the nearest 
future, thanks to increasing of engine efficiency, the situation will be not change. 

This paper presents a novel way of using alcohols as fuels for a diesel engine. Namely it 
shows the use of heavy alcohols as a solvent for vegetable oil (this mixture is called the bio- 
mix or BM) and after the obtainment the density by biomix which would be similar to the 
density of diesel fuel, mixing the biomix (BM) with diesel fuel (D) to obtain biomixdiesel 
(BMD). 

This solution will be shown for example with butanol as heavy alcohol, rape oil as vegeta- 
ble oil and conventional diesel fuel. 


Investigations 


The technology of butanol production is known. There are three separated ways to use. 
=" Directly 
o Conventional 
Traditionally, the production of fermentation products such as Bioethanol has relied 
on feedstock’s that are rich in either sugars (cane or beet) or starch which is easily 
broken down into sugars (wheat, corn or rice). The future is in new plants the work- 
ing with another art of microorganisms like thermophiles. 
o TMO technology process 
TMO's technology platform is based on a select group of thermophilic micro- 
organisms (Thermophiles). The optimal feedstock for bioconversions would be waste 
biomass (e.g. straw, wood chips and paper pulp effluent) and crops specially grown 
for their high biomass production rate (kenaf, miscanthus and short rotation woody 
crops). Such sources can be described as “cellulosic biomass” for their high cellulose 
and hemicellulose content. 
= Indirectly — from ethanol 
Electrolysis. The electrolysis is an old process, but with develops a fuel cell and new 
art of solar cell the electrolysis process may by develop to the new process of butanol 
production from ethanol. The elektrolyse is going from chemical eqation; 


2C2H;OH + 2C;HsOH — 2C4HoOH + 2H> + O2 


The properties of butanol and its comparison to the light alcohols and conventional engine 
fuels are as following. 
Tab.1. The properties of different fluids 


>3 cSt 
0.40.8 cSt 
29.2 MJ/1 0.43 MJ/kg 3.64 cSt 
19.6 MI/I 0.92 MJ/kg 1.52 cSt 
Methanol 16.0 MJ/1 1.20 MJ/kg 0.64 cSt 
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It is interesting that the butanol is similar in the energetically properties to the petrol. The 
solvent properties of butanol for solving heavy hydrocarbons (such diesel fuels) are very 
good. The mix is homogenous (both fluids don’t separate after several months). This is in 
contrary to the ethanol, which doesn’t solve the diesel fuel. It is important that the butanol 
practically doesn’t solve water (contrary to the ethanol which solves water in any proportion). 

There were prepared three mixtures for investigation. 

In the first step the rape oil (vegetable oil) is mixed with butanol from which we have the 
mixture, the density of which would be similar to the density of diesel fuel. This mixture is 
called as a BM (Bio Mix). In the second step this fuel (BM) was mixed with conventional 
diesel fuel (EN 590) to obtain biomixdiesel (BMD). These fluids were mixed in the following 
proportions. 

e biomix (BM) 10 % v/v, diesel fuel (D) 90 % v/v, called as 10BMD, 

e biomix (BM) 20 % v/v, diesel fuel (D) 80 % v/v, called as 20BMD, 

e biomix (BM) 50 % v/v, diesel fuel (D) 50 % v/v, called as 50BMD. 

All above mixtures are very homogeneous. It is impossible to see any phases (contrary to, 
for example, the mix of rape methyl ester with conventional diesel fuel). 

The engine investigations are carried out with a simple diesel engine on the engine test 
bed. 

Main parameters of engine (power output, torque, specific fuel consumption) and the main 
exhaust gas component (in this case CO, NOx, PM) will be evaluated. 

The investigation results are given in the following figures. 
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Fig. 1. The change of engine output parameters as a dependant factor of crankshaft rotation speed of the PER- 
KINS AD.3.152.UR engine supplied with 10BMD. 


The additive of 10% ,,bio” phase to the mineral diesel fuel leads to the remarkable in- 
crease of engine torque by unchangeable specific fuel consumption. Consequently the emis- 
sion of NO, decreases by 30% to 60%. At the same time emissions of CO and PM decrease, 
too. 

Likewise the parameters change by supplying the engine with 20BMD. 
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Fig. 2. The change of engine output parameters as a dependant factor of crankshaft rotation speed of the PER- 
KINS AD.3.152.UR engine supplied with 20BMD. 


Undoubtedly the biggest changes of engine parameters occur with supplying the engine 
with 50BMD. In this case the maximum of engine power and engine torque increases with the 
simultaneous increase of specific fuel consumption. Fuel consumption increases only with the 
higher engine crankshaft speed. The fuel consumption increases by the percentage as the en- 
gine torque, but increasing of torque was registered with all engine crankshaft rotation speeds. 
All the emissions of toxic parts of engine exhaust gases were decreased considerably. 
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Fig. 3. The change of engine output parameters as a dependant factor of crankshafi rotation speed of the PER- 
KINS AD.3.152.UR engine supplied with 50BMD. 


The above results were obtained without any change of engine control parameters (the 
engine control parameters were the same as with supplying the engine with conventional di- 
esel fuel), so without any optimization of control parameters. It appears that after optimization 
of engine control parameters; the results will be much better. 


Conclusions 


There have been given the first results of investigations of application of heavy alcohols as 
an ingredient of diesel fuel. In this case there have been presented the first results of investiga- 
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tions of the mixture of butanol (as heavy alcohol), rape oil (as vegetable oil) and conventional 
diesel fuel. This mixture was called the biomixdiesel (BMD). 

There were better results achieved than one expected. Contrary to existing experiences, 
the maximum of power output and the torque of engine is higher in the whole range of the 
rotatory speed of the engine crankshaft when the engine biomixdiesel (BMD) is reinforced. 
The addition of the component biomix to fuel influences the specific fuel consumption. Gen- 
erally with the larger part of the biomix component the specific fuel consumption grows. Be- 
cause the power of engine also grows up one should expect that in exploitation the specific 
fuel consumption should not increase. It is very important that this fuel could be used to rein- 
force old, existing now and the future diesel engines. 

The production of butanol is known (from biomass and in another way from electrolysis 
of ethanol). The possibility to get butanol from ethanol gives a very good perspective for the 
use of ethanol from today’s overproduction and moreover without the essential change of in- 
frastructure. 

There has also been presented another way to use alcohols (and vegetable oils — without 
transestrification) as diesel engine fuel. 

AII this leads to the conclusion that the use of ethanol overproduction and European pro- 
duction of vegetable oils will contribute to fulfilling of the expected requirements of European 
Union regarding the biofuels. 

It is very important that this fuel could be used to reinforce old, existing now and the fu- 
ture diesel engines. 

The future works should refer to, first of all, the better adjusting of the engine to fuel (es- 
pecially engine control parameters) and also fuel to the engine for specific exploitation needs. 
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Abstract 


The article presents the results of the experimental research of the heating cylindrical surface location influence 
on the local heat transfer coefficients in the bubbling fluidised bed in the physical model of the marine fluidised bed 
boiler. A particular feature of this boiler model is positioning it in a kind of cradle simulating ship’s motion on 
sinusoidal type regular wave. The heating surface is located in the traditional ship’s centre line or in her midship 
section. The testing has been conducted in the conditions where the fluidising column remained motionless while 
being at the same time deflected from the vertical by a constant angle and in the conditions of continuous oscillating 
movement. In the effect there has been obtained the distribution of the heat transfer local coefficient values in the 
fluidised bed in three points located at the diameter of the column at a certain distance from the separating grid with 
various locations of the heating surface. The research has indicated that column oscillations contribute to the increase 
of the bubbling fluidised bed local heat transfer coefficients, at the parallel positioning towards ship’s centre line. The 
further results of research may provide the grounds to formulate guidelines for designing the shipboard fluidised bed 
boilers. 


Keywords: heat exchange, fluidised bed boiler, design, ship 
1. Introduction 


While designing the marine steam boilers, certain rules due to ship’s rolling motion on waves 
should be observed in order to make the boiler operation safer and the risk if its damage gets less. 
In case of water-tube boilers with natural circulation such rule is for instance positioning of boiler 
drums in such a manner that their axis of symmetry is parallel to ship’s centre line. On the other 
hand, their positioning with the axis in perpendicular to ship’s centre line would cause excessive 
fluctuations of water level while rolling thus threatening the proper and safe boiler operation. 

Whereas in case of forced circulation boiler heating surfaces, e.g. evaporator or superheater 
consisting of pipe bunches, their location towards the ship’s centre line becomes less significant. 
There are arrangements encountered involving pipe positioning both rectangular and in parallel to 
this line. The ship’s side rolling has no such significant impact in this case onto the pump-forced 
working medium flow. Also in connection with exhaust gas flow in view of its strong turbulent 
motion and its minor inertia, the ship’s rolling motion practically does not affect the performance 
of the said heat exchangers. 
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On the other hand the situation changes completely if we consider the heating surface 
immersed in the bubbling fluidised bed formed by grains of fuel and inert material suspended in 
the stream of flowing air. Such is the case with fluidised bed boilers. Under the influence of ship’s 
rolling such bubbling bed behaves similar to liquids. Also the movement of grains in the bed itself 
changes which influences the local coefficients of heat transfer in the bed [5]. In this situation it is 
interesting to investigate the influence of the positioning of heating surface on the values of the 
local coefficients of heat transfer within the bed. 


2. The Course of the Research 


In order to find the answer to the question asked in the introduction several experiments have 
been conducted on the test stand specified in detail in [2]. The same method of determination of 
the local coefficients of heat transfer has been applied as referred to in the studies [5, 2]. The 
equivalent to a pipe in boiler bunch is here a cylindrical heating element (copper pipe) with 
electric spiral placed inside. The scope of research has required conducting of series of 
measurements for two various pipe positions, i.e. while the pipe axis is situated in parallel to 
midship section, and then when the pipe axis is situated in parallel to ship’s centre line. In the 
present construction of the stand the pipe can only be positioned in such a manner that its axis 
remains perpendicular to the column axis. This restriction results from the circular column cross- 
section. In the future it is planned to apply new version of the stand with the column of square 
cross-section which will allow to locate the pipe along the side walls of the column. The diagram 
of both positions of the pipe, distribution of thermocouples and the directions of movement of 
fluidising column are shown in Fig. 1, and the Fig. 2 demonstrates the view of a stand part. 


i 1 A. 
I I 
| 4 l 4 
3 2 1 I 
T T 
4 20 20 l 
+ =a i Gy 
a d A a a 
= uj — 3 | 
l | 
I a I al 
I a i a 
| l 
EM | © SEA 
| | 
sb 3 e wae 
=== ship's centre line a 
f midship section | 
| Ś 
J L 
+ — | 
of mm > 
J f i 
I 1 — fluidising column Å 
| 2 — heating element 
I ; f 
i 3 — separating grid 


i 4 — thermocouples 


[l 

l 

-| 

| 

—> directions of deflections | 


Fig. 1. Diagrams of the positioning of heating element and distribution of thermocouples 
3. The Results of the Research and their Analysis 
At the stand in question there have been previously tested the local heat transfer coefficients in 


the bed involving one permanent positioning of heating pipe, i.e. in parallel to midship section, 
whereas the bed height in rest condition and number of fluidisations have been changed. The 
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influence of column deflection from vertical has been also investigated with regards to the values 
of the local heat transfer coefficients. A major conclusion from this investigation has been the 
statement that the column deflection from vertical results in reduction of mean values of heat 
transfer coefficients, while the reduction tended to grow at the side towards which the deflection 
occurred. The results of these investigations have been presented inter alia in [5]. 


fluidising column 
thermocouples 
heating element 


fragment of grain return layout 


fluidised bed in rest 


Fig. 2. View of the stand 


The first test upon the change of heater positioning to that in parallel to ship’s centre line has 
been the checking of the influence of the number of fluidisations on the values of heat transfer 
coefficients in the bubbling fluidised bed. The examinations have been conducted for the 
permanent bed height in rest, i.e. 0.12 m. The bed material, similar like in the former testing, have 
been the poppy seeds. The heater power has been kept at the constant level of 45 W. According to 
the suppositions, in view of full symmetry of the system, in the column vertical position, the 
nature of the changes of heat transfer coefficients has remained the same as in the positioning 
referred to above — heater parallel to midship section. As the fluidising air velocities grow, also the 
heat transfer coefficients grow as well. Their distribution along the radius R of the column for the 
velocities 2.6 m/s and 4.5 m/s at the distance of 0.245 m over the separating grid in points 
separated by the value r from the column axis is presented in Fig. 3. 

As can be observed from the course of the curves, the values of the local coefficients of heat 
transfer are somewhat smaller closer to the walls on account of the circulation occurring inside the 
bed caused by well entry loss. 

The subsequent measurement series have been conducted for the column inclined by 
permanent angle of 22 ° to the left, and then inclined permanently by the angle of 28 ° to the right. 
The investigation has been conducted for the fluidising air velocity of 4.5 m/s. The results are 
presented in Fig. 4. They indicate the symmetry of the distribution of the coefficient value of heat 
transfer regardless from the direction of column inclination (V — vertical, L — left, R — right). Thus 
in this case the characteristic reduction of coefficient values at the inclination direction does not 
occur. On the other hand there can be noticed general reduction of all coefficient values at the 
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column inclination. The bigger decrease of the values can be observed at a bigger inclination , in 
this case for the inclination towards the right (alpha 4.5-R). But a slight reduction of the coefficient 
values in each case at the right side of the column results from the constant difference in the 
indications by thermocouples 1 and 3 (Fig. 1). 
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Fig. 3. The values of the local coefficients of heat transfer along the line at a distance of 0.245 m over the separating 
grid with the column in vertical position 
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Fig. 4. The values of the local coefficients of heat transfer along the line at a distance of 0.245 m over the separating 
grid with the column vertical position (V), inclinations to the right (R) and to the left (L) 


Fig. 5 presents the diagrams showing the distribution of the values of the local heat transfer 
coefficients for the same fluidising air velocity of 2.4 m/s with two different heater positions, i.e. 
in parallel to ship’s centre line (alfa 2,4-R longitudinally) and rectangular to it (alfa 2,4-R 
transversely) in another manner, in parallel to midship section, with column permanently inclined 
to the right. To make the comparison better there is also presented a diagram showing the 
distribution of the values of local heat transfer coefficients with the column positioned vertically 
(alfa 2,4 V). As already mentioned with the vertical column positioning the heater orientation is of 
no consequence. 

The same measurements have been repeated for the column inclined to the left with the same 
fluidising air velocity, i.e. 2.4 m/s. Fig. 6 presents the results in the form of diagrams. The line 
marked as “alfa2,4-L longitudinally” shows the distribution of the values of heat transfer 
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coefficients with heater positioned in parallel to ship’s centre line, whereas the line marked as 
* alfa2,4 — L transversely” with the positioning rectangular to this centre line. Similar as in Fig. 5 
also here a diagram is presented to show the distribution of the local values of heat transfer 
coefficients with column vertical positioning (alfa 2.4 V). 
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Fig. 5. The values of the local coefficients of heat transfer along the line at a distance of 0.245 m over the separating 
grid with the column vertical position (V) and column permanently inclined to the right (R) 
and heater transverse and longitudinal positioning in relation midship section 
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Fig. 6. The values of the local coefficients of heat transfer along the line at a distance of 0.245 m over the separating 


grid with the column vertical position (V) and column permanently inclined to the left (L) 
and heater positioned longitudinally or transversely in relation to ship 's centre line 


As visible in figures 5 and 6, in case the heater is positioned transversely to the ship’s centre 
line, there can be observed previously mentioned typical reduction of the value of the local heat 
transfer coefficient at the side to which the inclination has occurred. However, the most 
meaningful conclusion is that with the heater positioned transversely the mean values of the local 
heat transfer coefficient, regardless of the column inclination direction, along through the line of 
the heater are less than the values of the coefficients along the heater while it is positioned 
longitudinally to ship’s centre line. 

In the maritime practice the conditions when a vessel sails with a permanent major list to one 
side are quite rare. More complex ship’s movements on the waves are however general. In many 
shipbuilding issues, however, it is sufficient to take into account only so called simple rolling, i.e. 
of one degree of freedom, occurring only in ship’s centre line or midship section. Such approach is 
most frequently related to ship’s rolling movements which on the their intensity due to direct 
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relation with the stability safety and also the effectiveness of the performance of some systems and 
operations are of special significance [3]. Simple rolling without being combined with any other 
movements is however a little too far fetched approximation of actual conditions, but still 
acceptable in the research of the behaviour of the fluidised bed in these conditions which is proven 
also in other research studies, e.g. [4]. 

The stand where the aforesaid testing has been made allows to conduct the examination of the 
average values of the local coefficients of heat transfer during the constant oscillating movement 
approximating just the ship’s rolling on a regular wave of sinusoidal type. It is also assumed that 
the vessel is on even keel. 

The examinations have been conducted for the transverse and longitudinal heater positioning 
with constant oscillating movement of 28 s period and inclinations of 22 ° to each side. The figure 
7 shows the results in the form of diagrams. The line marked as “alfa L-R longitudinally” shows 
the distribution of the values of heat transfer with heater positioned in parallel to ship’s centre line 
during pendular motion, whereas the line marked as “alfa L-R transversely” with its positioning 
rectangular to this line also during pendular motion. 
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Fig. 7. The values of the local heat transfer coefficients along the line at a distance of 0.245 m over the separating 
grid with the column positioned vertically (V) and during column pendular movements with transverse or longitudinal 
heater positioning in relation to ship’s centre line 


For the sake of comparison also diagram is presented that shows the distribution of the values 
of the local heat transfer coefficients with permanent motionless vertical column positioning and 
air velocity of 2.4 m/s (alfa 2,4-V). During the column oscillating motion the fluidising air 
velocity has changed by itself within 2.4 up to 3.2 m/s. The velocity fluctuations are caused by the 
varying resistance characteristics of the fluidised bed during inclinations. The fluidised bed 
location out of parallel in relation to the separating grid during inclinations causes the drop in flow 
resistance and faster transfer into fluidised state, as well as increase in air velocity. These issues 
have been discussed more in detail in the study [1]. 

As can be observed from the courses in the diagrams the smallest average values of the local 
heat transfer coefficients along the line at a distance of 0.245 m above the separating grid occur 
during oscillation with the transverse heater positioning towards ship’s centre line. They are 
however bigger than with the constant list to port or starboard, in such heater positioning with the 
fluidising air velocity of 2.4 m/s (figures 5 and 6). The biggest values of the local heat transfer 
coefficients are reached with the longitudinal heater positioning and during column oscillation 
(alfa L-R longitudinally). Therefore the conclusion can be drawn that the column oscillations 
contribute to the increase of the local heat transfer coefficients, with parallel heater positioning in 
relation to ship’s centre line, in comparison to the values of the coefficients obtained in motionless 
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column positioned vertically (alfa 2.4-V). Another reason of such result are the local changes in 
the material concentration within heater area. 


4. Conclusions 


On the basis of the presented research results it can be concluded that the manner of heater 
positioning in relation to ship’s centre line does influence the average values of the local 
coefficients of heat transfer with permanent list or in continuous oscillating movement. In the 
examined positions the arrangement with the heater in ship’s centre line has proven better. 

On the basis of the results obtained, however, it cannot be unambiguously stated that in the 
construction of marine fluidised bed boilers — on account of heat exchange conditions — there 
should be preferred the longitudinal positioning of the pipes of heat exchangers immersed in bed 
in relation to ship’s centre line. It is necessary to complete the examination of the heat transfer 
coefficients with the heater also positioned in parallel planes to ship’s centre line and the examined 
midship section with the application of fluidising column of square cross-section. 

The conclusions from the above investigations, conducted with the use of certain 
simplifications, give the grounds to continue experimenting with the application of more advanced 
simulators of ship’s full movement on irregular wave. It also seems worthwhile to extend the 
scope of research by the inclusion of circulating fluidised bed. 
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Abstract 
Technical systems are more complex every day as their electronics and mechanics. Technological advances 
tend to be autonomous in its performance and perform an auto-diagnosis that allows determining an abnormality 
existence in a component or subsystem and deciding if the system has to be stopped or not. The conventional 
maintenance does not allow an integrated diagnosis analysis of a system. Among the factors that generated condition 
can be found a lack of communication between units: bad information management, ignoring relevant information, a 
lack of a clear monitoring policy and variables tendency. 


Keywords: diagnostics, monitoring, system of exploitation, reliability, damages 


1. Introduction 


The energy processors theory is based on a main energy flow analysis, where a system balance arises between input 
energy N, dissipated energy N, and useful energy N,. 

Through a residual process series as vibration, noise and heat the input energy is dissipated in one of these 
phenomena, that reflect the technical system wear (accumulated dissipated energy), therefore the dissipated energy study 
through the system provides inference about the artefacts wear and therefore is a interest determining the system technical 
condition. Figure 1, shows the methodology followed by the mini-central technical diagnosis. This process involved the 
next diagnosis stages: 

e Register and acquisition of residual energy (reception and selection points of vibration signals and operation 
variables). 

e Signal processing. 

e Statehood monitoring of Francis turbine. 

e Multidimensional Monitoring Condition. 
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e Maintenance ratios establishment. 
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Fig.l. Technical diagnosis methodology for a Francis turbine in slow time 


2. Study Case 

As study case for presented methodology is La Herradura’s Mini-central Hydroelectric property of Empresas Públicas 
de Medellin, which is located in the municipality of Frontino, near of Medellin, Antioquia. The mini-central has two 
Francis type turbines of horizontal axis, each one with a rated power of 10.4MW and a 5m'/s flow, with a rotation speed of 
900rpm and a design net jump of 230.6r.p.m. Fig.2 shows up a general scheme of parts forming the turbine. 
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Fig.2. Francis turbine 


The Mini-central has two systems able to obtain information of Mini-central technical condition, the first one is the 
vibration monitoring system and the second one is the Mini-central monitoring and control system. 
The permanent vibration monitoring system for the generator is based on an instrument with the serial number “VDR-24” 
(Vibro Diagnostics Recorder — 24 channels), in the VDM data module and in the “ATLANT” diagnosis program. The 
vibrations measurement chain is showed on Fig.3. 


170 


Fig.3. Measurement chain, vibration analysis system [1) La Herradura’s Mini-central. 2) Sensors — Variable Measurement. 3) VDM- 
Data Module. 4) VDR — Processing and signal acquisition.. 5) ATLANT Software. 6) ATLANT signal and analysis transformation] 


Through this system the r.m.s signals vibration value is monitored (speed and displacement) in points (1-8) presented 
on Fig.4. 


Fig.4. Acquisition point scheme of variables related to vibration (1-8 points) 


The system used to central monitoring and control from operation station is the V7 Monitor Pro from Schneider 
Electric (Fig.5), this allows the data acquisition, monitoring and real time control and has a setting Server-Client and an 
unlimited number of TAG’s (variables). 


Fig.5. Measurement chain, monitoring and control system 
[1) La Herradura Mini-central. 2) Variables measurement sensors (Pressure, temperature, current, voltage). 3) Control panel. 4) V7 Monitor Pro 
Software. 5) General deployments of generation and monitoring units. 6) Historical board and reports] 


3. Measurement points selection 
Based on measurement points algorithm the mechanical vibration signal is analyzed, received from the three hydro- 
dynamical bearings based on independence criterion and information quantity. Fort the first case the information 


independence will be given by inverse area under the curve of coherence oa ( f ) between two signals measured at the 


same place, for the current study the signals are taken from vertical and horizaontal speed in every bearing therefore, there 
will be a greater information independence when the maximum area, according to the next expression [5]: 
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[uF 
(1) 


To determine the information quantity the coherence values are taken between signals depending on certain frequencies 
(system characteristic frequencies) and a criterion under the following expression [5]: 


l 
1 
In,,(6;)= Log 


AC xy and In,,, values are registered in a data base, given the data volume to analyze, an optimization problem is set 


out, which aims to determine the generation bearing in which the independence and information quantity are the highest. 
According to Fig.6, the reception points of diagnosis signal are located almost at the same distance from the optimal point. 
This means that has reliable information for the technical diagnosis of the generation unit. 
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Fig.6. Reception point analysis from diagnosis signal 
4. Symptoms Calculation 


During the diagnosis model implementation, a series of new symptoms were calculated to make a registered signals 
by the vibration monitoring system, Fig. 7 shows an example of the spectrum amplitude (225Hz):2VV symptom: this refers 
to the vibration speed amplitude of impeller blades frequency flow (225Hz), in vertical direction of 2 bearing. A data 
tendency was observed during monitoring time, which indicates the evidence of a system abnormality, thus the evidence 
the real system statehood condition, allows detecting, locating and evaluating failures in the system. On [6] is showed the 


symptom definition and some examples are suggested. 
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Fig.7. Vibration speed frequency from impeller blades flow (225Hz), in a vertical direction of 2 bearing 
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5. Observation Matrix Elaboration 


To elaborate the symptom matrix 25 generation and control variables were considered, corresponding to specific 
average values registered during the monitoring day. This selection was made along with the systems operators, seeking to 
include the variables that in a certain case can provide abnormality evidence in the system. Regarding to vibration 
monitoring signals, 210 symptoms were calculated each monitoring day, also based on operators experience and 
appropriate literature. Among them highlights scalar estimators as: average, RMS value, peak value, peak to peak value, 
shape factor, standard deviation, bias, among others. In this way a new observation matrix of symptoms will contain 235 
variables and 157 observations. The symptoms observation matrix from a system is represented as follows [7]: 


Sintomas 
Varaibles 
GOA Si; Si | >80 
Sa Saa 0 Sag Sor | 6, 
5 i ) 3 : : Observacione: 
pr i Sa Sip ce Sy co Su | >O, Mediciones 
Sra Spa Spo Spr] 70, (3) 


Where columns j =1,2,---,7 are different measured symptoms and rows i=1,2,-:-,p are observations or 
measurements made for every symptom in different life cycle times of technical system. 


6. Limit value establishment for estimators 


Symptoms limit value of diagnosis systems were calculated according to the following relation: [5]: 


Shim =S+ Op 
2A (4) 
Where: § is the symptom average value during machine operation time 0 and the symptom standard deviation op, A - is 


the tolerable level of unnecessary established repairs and G is the machine availability. 

Fig.8 shows the relation between the manufacture established limit and the previous calculated method. The set limit 
observed can be found way above from normal data behavior therefore do not show variable changes evidence. On the 
other hand, the calculated limit can identify subtle variable changes being in an historical behavior range. 
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Fig.8. New limit calculation 
7. Singular Value Decomposition 


Following the proposed methodology, symptoms observation matrix, its dimensions are reduced through PCA. Then 

the Singular Vale Decomposition (SVD) is applied with the purpose of extracting different failure modes that evolve in a 

system, assessing the wear advance used en new indexes and ratios. The SVD an application for sizing the symptoms 
observation matrix can be expressed as follows [8]: 

= r 

0, =U, *#Di n V] 8 


U P> dimension orthogonal matrix. p , are the left singular vectors. V,„: is an Fr dimension orthogonal matrix of right 


singular vectors. Lie : is a diagonal matrix of singular values. 
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The failures profiles are determined using singular values and vectors O,,U,,V, found with SVD, obtaining a 

conditon evolution interpretation of technical system. These failures are given by [9]: 

SD, = O,, XV, =O; "Uy © 

Where SD, is the left singular vector amplified by a respective singular valueo,. Hence this value leads as bug and 
information about intensity of failures due to the inclusion of o, [9]. 


The total generalized failure profile P(A) or SumSD , which represents the general evolution of condition of 
technical system is determined through [10,11]: 


P(0)= SumSD = yD, (o) 
= (7) 


8. Implementation Methodology 


During the implementation methodology at La Herradura’s Mini-central an evolving failure on Francis turbine was 
detected, which is showed on Fig.9. 
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Fig.9. First Francis turbine failure evolution 


The technical diagnosis showed up that the evolution profile of machine was strongly correlated with variables that 
describe the technical condition of first Francis turbine bearing. The temperature, the relative axis displacement with 
respect to combined bearing and the 225 Hz spectral component are part of diagnosis parameters that dominate in this first 
evaluation unit. A continuous component increase was observed during the machine operation time. This occurs as a 
consequence of interaction between impeller blades and distributor moving blades, a pulse is generated due to the 
frequency flow pressure of impeller blades (225 Hz, this pulse is labirynths transported from turbine seals causing an axis 
push in axial sense, generating a vibration at the same frequency level. With the turbine seals wearing increase, the pulse 
effect increases, hence, the axial push increases generating vibrations increase. 

On Fig.10 the diagnosis parameters tendencies are showed for Francis turbine operating time related to identify 
failure. 
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Fig.10. Vibrations evolution related to failure profile 
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The probabilistic decision model (Fig.11) and reliability symptoms function (Fig.12) of Francis turbine were 
determined with important information in any strategy of critical operating systems maintenance. These machine’s 
behavior patterns allow making correct decisions just in time and reducing risk. It is important to remember the 
implemented methodology during this project, which is based on real data of Francis turbine condition during utility time. 


- Probabilidad de acierto en la 
evaluación de buen estado 88%. 


- Probabilidad de acierto en la 
evaluación de la averla 84% 


- Probabilidad de no detección de la 
averla 16% 


- Probabilidad de falsa alarma 12%. 


Fig.11. Probabilistic diagnosis model 
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Fig.12. Reliability symptoms function of Francis turbine 
9. Conclusions 


Techniques and algorithms used in different branches of science can be applied to technical condition monitoring, as 
the case of main components analysis and the singular values decomposition. 

The multidimensional symptoms monitoring allows identifying changes in the system technical condition and 
establish possible causes from that condition. 

This kind of monitoring in the specific analyzed case generate a maintenance decision-making support, which impacts 
in cost reduction related to maintenance and optimal personal use besides, it generates an increase in the system’s 
availability and reliability. 

The study was made in real exploitation conditions, considering dynamic variables of generators with the purpose of 
obtaining information about the general technical condition of system. 
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Abstract 

The multidimensional monitoring of symptoms applied to railway systems allow to detect and 
locate the sections of the track (straight and curved) that generate the decrement of the safety and 
comfort of the passengers. It also evaluates the technical state of the rail-vehicle interface. In 
addition, it allows observing, evaluating, and controlling the reliability and availability of the system. 
The objective of the study is to propose an alternative to evaluate the condition of the technical state 
of railway systems from a dynamical point of view which guarantees safety and comfort of the 
passengers. One looks for diminishing operative costs of maintenance, improving the use of 
equipment for maintenance tasks of the track, vehicle and auxiliary equipment, optimizing the time of 
the maintenance personal, the maintenance frequencies (corrective, preventive, etc.). It also aims to 
identify the variables related to maintenance actions that have a high influence on the technical state 
of the system. This paper presents the results obtained when applying a modeling of this type to a 
railway system, being centered mainly in the application of SVD theory to the technical diagnosis of 
systems. 


Keywords: multidimensional monitoring of condition (MMC), Singular Value Decomposition (SVD), 
railway technical diagnosis 


1. Introduction 


Technical systems are becoming more complex when talking about their mechanics and 
electronics. Technologic advances tend to be more auto-sufficient, and able to auto-diagnose 
themselves, which allows to determine if any anomaly is present in any subsystem or component, to 
finally decide whether the system must or must not be stopped. 

The conventional maintenance with some factors, such as the lack of communication between 
dependencies, a not proper management of information, not having a clear monitoring policy or 
variables trends among other, do not allow the performance of an integrated diagnosis of the system. 

Due to the previous factors, it is necessary to implement new methodologies of technical 
diagnostics, in order to satisfy all of the company’s needs for obtaining results, an integrated diagnosis 
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through computing simulation, tools, analysis methods, and information evaluation from the 
machine’s technical state. 

The multidimensional monitoring of symptoms — MMC — has been developed as an alternative 
that fits the new technology demands; it assures obtaining performance indicators such as the 
reliability before any damage within the system, giving as a result less critical failures, unexpected 
pauses of production, dead times, maintenance costs, and the mistreatment of human resources, 
among others. 
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Fig. 1. Diagnosis Model with a holistic approximation for the railways systems 


The holistic approximation presented in the previous figure, obtains information about the dynamic 
system operation, and the maintaining activity history through: 
e The portable diagnose system (PDS) which registers acceleration and force signals in order to 
evaluate the security, comfort, and railways fatigue according to the International Law 
UICS5 18". 
e The managing strategies of the company, taking into account the permanent way and vehicle 
variables. 
e A geometric modeling of the wheel-rail interface, through which it is possible to obtain the 
equivalent iconicity and the contact efforts, which are the factors describing the relation 
mentioned above. 


1 UIC518, Tests and approval of rail vehicles from the point of view of their dynamic behavior, safety, railway 
wear and quality of the travel (comfort) 
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e = The holistic approximation, based on the application of MMC methodology which is based on 
the singular values decomposition method (SVD) as well. SVD determines the general failure 
of the system, which is the indicator of the integral technical state of the system. Its evaluation 
allows the improvement of the actual management actions and the rise in the availability, 
reliability, and maintainability indicators of the technical system. 


2. Multidimensional monitoring of condition 


Figure 2 shows the implemented methodology for this study case, 4 main groups are warned; in 
the first one, the energy processors theory is reflected, where through the exit energy study, the 
technical state of the system and its deterioration can be inferred. First, signal processing (second 
group) must be done due to the fact that phenomena like vibration, noise and heat can be measured by 
data acquisition equipment, and instrumentation (symptom monitoring). 
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Fig. 2. Methodology for holistic approximation 


At this point, it is important to establish what kind of measuring points of the physical variables 
are appropriate in order to obtain relevant information, that is why a selection of measuring points will 
be performed to establish the point generating a larger amount of information, and a smaller amount 
of useless information about the technical state of the system. 
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Once the measuring points are established, the monitoring of the physical variables can be made, 
assuring relevant information about the system. 

Acquired signals are then transformed into a symptom series. The calculation parameters of these 
symptoms might be obtained by any norm alignment, or by a temporal series statistic analysis, in 
which basic statistic scalars can be calculated, among many others, as the half, standard deviation, 
bias, form factor, clearance factor. Way more complex analysis like the spectral, the envolvent, the 
cepstral, frequence time, ciclostationarie and superior order statistic analysis might also be developed 
in order to obtain symptoms which reflect the measuring signal. For the present study, the signal 
processing is being done through SPD. 

A third group within the methodology is composed of the technical system’s proper information, 
like system operation variables, system parameters, and a technical system geometric modelation. 
Among those sources of information, a series of periodically measured variables can be found 

To perform a multidimensional monitoring is necessary as the first step to determine the symptom 
observation matrix which is composed of the symptoms in its columns and the measures of those 
symptoms in its rows. This matrix contains information about the system’s technical state and the 
proposed methodology seeks to extract the relevant information. 

The observation matrix is reduced through a Principal Component Analysis PCA which, along 
with a series of reduction criteria, allows identifying the symptoms containing relevant information 
about the system’s technical state. 

Based on the reduced matrix (space of reduced symptoms), and by using the Singular Value 
Decompositions SVD, the general failure that presents the evolution of the system’s technical state is 
determined. It is a variable which consolidates the information about the system’s technical state, 
contained in principle in the symptom observation matrix. 

Once the general failure and its limit value have been determined (through statistical methods), it 
is possible to perform an integral diagnosis of the system by estimating the conditions under which the 
general failure exceeds its limit value. Forecasting tasks can be made, beginning with the general 
failure, in order to predict the future technical state of the system. It is also possible to create 
indicators of the system’s technical state from the general failure, such as the functions of risk and 
reliability which support the system diagnosis. And finally, it is possible to establish general models 
(through multiple linear regression) which describe the behavior of the general failure in order to enter 
current, future, real or assumed values to the model and so determine their influence over the general 
failure. 


3. SVD General Theory 


The symptoms observation matrix defined by O,, € R™” is an array of symptoms at certain 
moments of the life of the technical devicef It can be said that the symptoms observation matrix is a 


discrete way of observing symptoms (Cempel, 2000, p.2). 
The aim of obtaining a symptoms observation matrix O,, is to obtain all the information 


related to the condition of state of the technical system and to distinguish different forms of failure 
that evolve during its operation (Cempel, 2000, p.2). The symptoms observation matrix of a system is 
represented as follows (Natke-Cempel, 2001, p.599): 
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Variable 


Symptoms 
EM NE =~ NY Sis >b, 
Sy) NE ate NY cn NE => 0, 
: : k : : Observations i=1,2,---, 
O,, ={8,}= p 
ý | Sa Sip ct Si, Si | >6,| Measurements j =1,2,---,r 
[Spa Sa Spg ooo Si. > 8, 


Where columns j=1,2,:--,r are the measured symptoms and rows i=1,2,:::,p are the 
measurements made to each symptom at certain moments of the technical system’s life. 
In some cases the information of several columns can be redundant and varies among symptoms 
contained in the matrix O,, This can be avoided by removing the correspondent (redundant) columns 


using tools such as the PCA. 

The Singular Value Decomposition SVD is used in order to obtain different modes of failure 
which evolve in a system, evaluating the wear progress through indexes and indicators. 

The application of the SVD to dimension the symptoms observation matrix is expressed as follows 
(Bongers, 2004, p.42), (Cempel, 2000, p.4), (Cempel, 2004, p.2), (Cempel, 1999, p.179), (Cempel, 
2003, p.1293), (Cempel, 2003, p.2), (Cempel, 2000, p.5), (Cempel-Tabaszewski, 2006. p.3), (Cempel- 
Tabaszewski, 2003. p.216), (Cempel-Tabaszewski, 2005. p.3), (Wall-Rechtsteiner-Rocha, 2003, p.2): 

Op = U 3 Że M (1) 
U „: is an orthogonal matrix of dimension p of the left singular vectors 


V „: is an orthogonal matrix of dimensionr of the right singular vectors 
È „ : is a diagonal matrix of the singular values: 


in diag(0,,0,,-*:,0,) (2) 
The singular valueso, are expressed by numbers and u,,v, are the singular (orthogonal) vectors 


U,V... At the same time, these 


representing the columns of their corresponding matrices > pe es 


pr? 


vectors are creating independent sub-matrices (o, ) which totally describe the modifications of the 


operational system, that is its evolution, wear or its failure modes at time ¢ =1,2,---,z (Cempel, 
2003, p.1293). 
The singular values different from zero o, >0 indicate failure or damage. These values are used 


to detect changes in the system and to evaluate their intensities. Failure intensities are organized by 
their magnitudes, as shown in the main diagonal of the singular values of O, (Natke-Cempel, 2001, 


p.610). 

General failures are determined by using the singular values and vectors o,,u,,v, found 
through the SVD, obtaining an interpretation of the evolution of the state condition of the technical 
system. These technical failures are given by (Cempel, 2000, p.4), (Cempel, 2004, p.2), (Cempel, 
1999, p.179), (Cempel, 2003, p.1293), (Cempel, 2003, p.2), (Cempel, 2000, p.5), (Cempel- 
Tabaszewski, 2006. p.3), (Cempel-Tabaszewski, 2003. p.216), (Cempel-Tabaszewski, 2005. p.3), 
(Wall-rechtsteiner-Rocha, 2003, p.2): 
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SD, =O,,XV,=0,*u, (3) 
where SD, is the left singular value amplified by the corresponding singular valueo,. Hence, this 


value has independent information in the form of failure as well as information in the intensity of 
these failures because of the inclusion ofo,. 


Therefore, for a lifetime 6 of the system, quantities SD, (8) and o, (8) are independent from one 


another and act as indexes of the system changes during the operation. 
From the viewpoint of energy equivalence of a norm, the following can be said: o, (0) can be 


treated as an index of the wear progress (intensity), while the SD, (8) is the momentary form during 
the evolution of the condition of state of the system (Cempel, 2003, p.1294). 

One of the indexes which can be obtained through the SVD is the profile of total general failure 
P(0) or SumSD which represents the evolution of the general state of the technical system and is 
estimated as (Cempel, 2003, p.1296), (Cempel-Tabaszewski, 2005. p.4): 


P(0)= SumSD = > SD,(0) (4) 


The profile of total general failure SumSD, will be used to describe the technical state of the system 
under research and from now on will be called general failure. 

Through the general failure, it is possible to perform diagnostic and forecast tasks, as well as 
determine some additional indicators of the condition of the technical system, and define state models. 


4. Implementation of the multidimensional condition monitoring 


We implement multidimensional monitoring of symptoms in a railway system type metro and 
present the information sources of the system’s technical state as well as the conditions considered for 
the implementation and the results obtained through a computational tool based on the proposed 
methodology. 


4.1 Description of the experiment 

The analysis focuses on a commercial railway and it is developed for a sample of approximately 
16 straight sections and 37 curved sections, as well as a representative sample of passenger vehicles 
(17 out of a 42 cars fleet). 

The related symptoms are the estimators considered by the UIC518 standard for the evaluation of 
security, comfort and wear of the railroad; the geometric properties of rail and vehicles as well as 
estimators related to the wheel-rail interface. Two groups of matrices are formed: 


i~1,2,---,16 
Estimadores Estimadores Estimadores aD 5 
J = 325 ... š 
O,, =O, = da — riel UIC f 
pr 16x32 “a p a k- 23,48,91 0,13,- 18 
2 > k 1=1,3,4,:--,15 


The first group is made of the estimators calculated from the UIC518 standard W, ,, the properties 
of the railroad Y,,, and the estimators related to the wheel-rail interface Z,, (for the lot of 31 


estimators). There is a sample of straight sections i as observations of these symptoms. So we have 
17 matrices of this type, one for each passenger vehicle analyzed. 
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The second group is similar to the first one, but here a sample of curved sections i is 
considered for a total of 35 estimators. 17 matrices of this type are obtained as well. 


i+ 1,2,---,37 
Estimadores Estimadores Estimadores 1-12... 5 
J — sł, ... > 
O,, =O. = d jel UIC fi 
pr 37x35 rue p o i j ty pz 1,4,6,8,9,10.1 113, o 20 
y ij i,k i, 1=12,:--15 


To contract the matrices mentioned before, we used information obtained from maintenance 
employees of the company, as well as from experimental measurements performed in normal 
operation conditions of the system. 


4.2 Experiment Results 

Each of the formed matrices is implemented into the multidimensional monitoring of condition 
MMC. As the first step, they are transformed considering the mean value and the standard deviation 
of each of the estimators; then, they are implemented into the dimensionality reduction module 
through the analysis in main components. Once the estimators providing little information about the 
technical state have been established, they are removed from the original matrix (with no 
transformations). 

The reduced matrix is again transformed, but this time with respect to the initial value of each 
estimator. It is included in the singular value decomposition module SVD through which the general 
failure is determined. This is an indicator of the general state of the technical system whose limit 
value is calculated by the Neyman-Pearson method, considering the availability of G=0.9 and the 
level of unnecessary repairs of A =0.1 

Once the general failure and its correspondent limit value have been calculated for each group of 
matrices, diagnosis is performed, the indicators of state are generated and its general models are 
proposed. We now present the obtained results for each group of matrices formed. 

Figure 3 shows the general failure obtained from an observation matrix where symptoms are the 
variables measured on the railroad, the estimators obtained after applying the SPD and the ones 
calculated for the wheel-rail interface. The observation of these symptoms was made in vehicle 05 
running through different curved sections of the railroad. 

Axis x represents the observations obtained from the general failure (curved sections) while axis 
y represents the amplitude of the failure. The diagnosis is, hence, developed when determining the 
railroad section in which the general failure exceeds its limit value; section where the company should 
consider performing maintenance. 


Figure 4 shows the co-variances (axis y ) between the general failure and the symptoms considered 


at the analysis (axis x ). Symptoms having a high relative covariance when compared to other 
symptoms are assumed to be responsible for the exceeded value of the general failure (with respect to 
its limit value). Therefore, once we identify the road section in which the general failure exceeds its 
limit value, we find through covariances the symptom generating this condition and, in this way, 
maintenance actions can be focused on particular road sections and on possible phenomena happening 
in these sections. 
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Fig. 3. Profile of the general failure 
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Fig. 4. Covariance of symptoms with the general failure 


4.3. Results for UIC estimators, railroad properties and estimators of the wheel-rail interface for 
straight sections. 

The procedure described before is applied to the 17 matrices in which the evolution of symptoms 
is given by the straight sections of the railroad. In 64.7% of the cases (analyzed units), the variable 
representing a greater importance is the railroad corrugation, 11.8% is due to the rail corrugation and 
its alignment, 11.8% is due to the alignments only, a 5.9% is due to the equivalent conicity and 5.9% 
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to the lateral forces. 
In 76.5% of the cases, the system’s reliability is reduced in the straight section ENV-86. In 70.6% of 
the cases, reliability is reduced in the straight section 77-81. This happens because the general failure 


exceeds its limit value. 


4.4 Results for UIC estimators, railroad properties and estimators of the wheel-rail interface for 
curved sections. 

In 41.2% of the cases (analyzed vehicles) the most relevant variable was the railroad corrugation, 
17.6% is due to the railroad corrugation and the horizontal alignment, 23.5% due to the lateral forces 
and 17.6% due to the quasi-static lateral accelerations. In 76.5% of the cases, the system’s reliability 
decreases in curve 2 and in 70.6% of the cases, reliability decreases in curve 5. 


4.5 Models of state for the general failure 

Once the general failure and its relationship with the measured symptoms has been determined, an 
analysis of multiple regression is developed considering the matrix formed by the estimators of the 
UICS18 standard, the geometric properties of the railroad and the estimators related to the wheel-road 
interface, for straight and curved sections. 

The general failure found after the application of the SVD is considered as a dependent variable. 
The most representative symptoms in each case (symptoms having a great participation on the general 
failure) are considered as independent variables. 

Particular models for each of the 17 analyzed vehicles are determined, and then each of these 
models is applied to real data from the other 16 vehicles in order to obtain a general model adaptable 
to most of the cases. 

We now present the general model obtained for the general failure SumSD from the straight 
sections of the railroad: 

SumSD = 18.74 + 24.80Y, , + 14.28Y,,; —11.36Y,,, — 7.917, „, —5.02W, „s 


— 4.62Y, ; + 3.09W, , +3.00Y, , + 2.49Y, , —1.26W,, — 0.97Y,, — 0.44Z,, — 0.207, ,, 


This model adapts with an acceptable global significance to 7 of the 17 vehicles under research (the 
level of significance is assumed a =0.05). Figure 5 shows the three best adaptations of the general 


model. 


(5) 
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Fig. 5. Representation of the three best adaptations of the general model obtained for straight sections 
of the railroad 
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The following is the general model obtained for the general failure SwmSD from the curved 
sections of the railroad: 
SumSD = 44.13 + 27.07Y, , + 8.00Y, |; + 7.78Y, „, + 4.30Z,, +1.80W, ,, 


+1.38Z,, + 0.32W, , — 0.31W,, — 0.13W, , +0.06W, 1s 


This model adapts with an acceptable global significance to 8 of the 17 vehicles under research 
(the level of significance is assumeda =0.05). Figure 6 shows the three best adaptations of the 
general model. 

By determining these models we tend to establish the condition of the general failure of the system 
for current conditions of the geometric parameters of the rail-vehicle interface. This allows to make a 
rapid diagnosis of the system in order to identify the necessary actions of maintenance. 


(6) 
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Fig. 6. A representation of the three best adaptations of the general model obtained for curved 
sections of the railroad 


5. Conclusions 


The multidimensional monitoring of symptoms helps us identify the railroad sections and vehicles 
needing a detailed revision by the personnel of railroad maintenance and movilization. It also gives us 
clues about what happens in the road sections under research. Hence, this type of monitoring in this 
specific case becomes a tool that supports maintenance decisions which translate into a reduction in 
maintenance-related costs, an optimal use of the maintenance personnel, and it additionally increases 
the availability and reliability of the system. 

The research was performed under normal conditions of operation considering dynamic variables 
of the vehicles, estimators related with the wheel-rail interface, and with parameters commonly used 
by the personnel of maintenance and movilization in order to obtain information about the general 
technical state of the system. 

Rails corrugation is the main factor affecting the reliability of the system in straight and curved 
sections since it influences the general failure, which lets us infer about the general technical state of 
the system. It is also known that the railroad corrugation produces vibrations on the suspended masses 
of the vehicle and such vibrations are transmitted to the passengers, causing discomfort. 
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Other influencing factors are the lateral forces measured in the axes’ extrems. These forces affect 
the stability and security of the vehicle when it travels through curved sections, as well as the railroad 
alignments and the quasi-static lateral accelerations measured in the vehicle’s box (which measures 
the influence over the passengers’ comfort). 

Road sections generating a loss of reliability on the system are: straight section ENV-86, straight 
section 77-81, straight section 85-87, curve 2, curve 5, curve 4 and curve 1. Therefore, we could 
suggest a change in the maintenance routines in order to give priority to these road sections. 

Two general models for the evaluation of the general failure in straight and curved sections were 
established. This was made in order to quantify the general failure in current conditions of the rail and 
the vehicle, supporting the maintenance activities. 
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